PROPAGATION  OP  HARMONIC  WAVES  IN 
LAYERED  TRANSVERSELY  ISOTROPIC 
SHELLS 


HENRY  EVANS  KECK 


UNIVERSITY  OF  FLORIDA 


DEDICATION 


From  women'e  eyes  this  doctrine  I derive: 

They  sparkle  still  the  right  Promethean  fire; 

They  ace  the  books,  the  acts,  the  academes, 

That  show,  contain,  and  nourish  all  Che  world; 

Else  none  at  all  in  aught  proves  excellent. 

Love’s  Labour's  Lost,  IV,  ill 
Shakes peace 


These  words,  spoken  for  me  by  Shakespeare,  reveal  nor 
true  feelings  as  1 dedicate  this  work  to  Jeanne,  ay  wife. 


B.E.K. 


ACKI^O^^LEDG^EN^S 


The  author  thanks,  with  sincere  appreciation, 

Dr.  A,  E.  Artnenlkas,  dissertation  chairmen  of  hie  super- 
visory conmltCec,  for  guidance  and  tutelage  In  the  coin- 
mission  of  this  work.  The  debt  here  can  not  be  estimated. 
Additional  thanks,  for  continuing  interest  and  concern  in 
the  author's  progress,  are  given  to  Dr.  I.  K.  Ebeioglu, 
committee  chairman,  Dr.  G.  E.  Neville,  Jr.,  Dr.  S.  Y.  Lu, 
Dr.  K.  Yee,  and  Dr.  L.  Janos,  members  of  the  supervisory 
conmittee. 

A separate  word  of  appreciation  is  proffered  to 
Dr.  K.  A.  Eisenberg  for  the  many  and  thoughtful  discus- 
sions, given  in  the  spontaneity  of  academic  interest. 

The  author  is  further  indebted  to  the  United  States 
Air  Force,  in  a manifold  manner:  for  providing  the  chance 

to  pursue  graduate  studies;  for  providing  the  assistance  of 
the  personnel  at  the  Flight  Dynamics  Laboratory  and  the 
Digital  Computation  Division,  Hright-Fatterson  Air  Force 
Base,  Ohio,  without  Whose  aid  numerical  results  could  not 
have  been  obtained. 


11 


TABLE  OF  COflTENTS 


ACKNOWLEDGMENTS  li 
LIST  OF  FIGURES  vi 
KEY  TO  SYMBOLS  vli 


ABSTRACT  xi 

CHAPTER 

I.  1 

II.  10 

1.  Basic  Equations  10 

2.  Displacement  Equations  of  Motion  16 

111.  20 


1.  On  Che  Use  of  Displacement  Potential 
Functions  to  Solve  the  Equations  of  Motion  20 

2.  A Completeness  and  Representation  Theorem 
for  General  Harmonic  Waves  in  Transversely 


Isotropic  Bodies  27 

3.  Alteration  of  the  Representation  for  Two 

Special  Cases  35 

4.  Specification  of  >f(r.e)  and  f(r,e)  39 

IV.  43 

1.  Analysis  of  the  Radial  Wave  Number  s,  and 

f(r)  43 

2.  Analysis  of  the  Radial  Wave  Numbers  p and 

q,  and  ?(r)  48 

V.  68 


1.  Preliminaries  of  the  Problem  and  Nondinen- 

sionalization  of  Equations  66 

2.  Derivation  of  the  Frequency  Equation  75 

3.  The  Cl-l;  Plane  85 

4.  Turning  Point  Line  Expressions  for  the  C^j  92 

5.  Complex  Zone  Line  Expressions  for  the  C^j  97 


ill 


TABLE  OF  CONTENTS  (Continued) 


the  Frequency  Equation  Representation 
in  Complex  Zones 


Frequency  Equation  for  Three  Layered, 
Traction  Free,  Isotropic  Cylinder 
Frequency  Equation  for  Haves  of 
Infinite  Length 

Frequency  Equation  for  Axisymmetric 
Traveling  Haves 

Frequency  Equation  for  Axisymnetric 
Haves  of  Infinite  Length 
Frequency  Equation  for  One  and  Two 
Layered  Shells 

layered  Cylinders  with  Clamped-Free 
and  Clamped-Clamped  Boundary  Conditions 
Frequency  Equations  for  Cylinders  with 
Clamped-Free  and  Clamped-Clamped  Bound- 
ary Conditions  When  ■?  = n =0 
Frequency  Equation  for  Cylinders  with 
Very  Thin  Outer  Layers 
N-Layered  Rods  and  Cylindrical  Shells 


The  Frequency  Equation  in  Sector  1 
for  y-»«o 

A.  The  Equation  Dg’  =0  (1  = 1,3) 

B.  The  Equation  Dj"’’’  =0  (i  = l,3) 

C.  SuMtiary  of  Sector  1 Limits 
The  Frequency  Equation  in  Sector  2 


A.  Sandwich  Shells  - Sector  2 

B.  Core  Layered  Shells  - Sector  2 

C.  General  Shells  - Sector  2 

The  Frequency  Equation  for  Higher  Sectors 
Suninary  of  the  Asymptotic  Behavior  of  the 
Frequency  Equation  When  ^»1 


The  Frequency  Equation  for  Infinite, 
Three  Layered  Plates 


98 

104 

104 

108 

113 

117 

121 

124 

132 

135 

139 

144 

148 

151 

159 

160 

161 

170 

173 

174 

176 

160 

180 


TABLE  OF  CONTENTS  (Continued) 


Page 

187 


1.  Method  of  the  Numerical  Analysis  187 

2.  Discussion  of  the  Frequency  Spectra  193 

3.  Discussion  of  the  Vibration  Displacement 

Distributions  200 

BIBLIOGRAPHY  218 

BIOGRAPHICAL  SKETCH  226 


LIST  OF  FIGURES 

Figure  Page 

1.  versus  (i  = 1,2,3)  54 

2.  Inequality  (4:24)  57 

3.  Radial  Wave  Numboirs  versus  60 

4.  Radial  Wave  Numbers  versus  62 

5.  Cress-Scction  of  Cylinder  69 

6.  Ci-t  Plane  with  Zone  Linos  91 

7.  The  veraus  b'"^  Plane  157 

6.  Plate  Cross-Section  181 

9.  Frequency  Spectrum:  Homogon.,  H/R  = 1,0  203 

10.  Frequency  Spectrum:  Hoisoscn.  (cont),  H/R  = 1.0  204 

11.  Frequency  Spectrum:  Case  I,  H/R=1.0  205 

12.  Frequency  Spectrum:  Case  1 (cont),  H/R  = 1.0  206 

13.  Frequency  Spectrum:  Case  II,  H/R=1.0  207 

14.  Frequency  Spectrum:  Case  II  (cont),  H/Rei.O  208 

15.  Frequency  Spectrum:  Homogen,,  H/RtO.l  209 

16.  Frequency  Spectrum:  Homogen.  (cont),  H/RcO.l  210 

17.  Frequency  Spectrum:  Case  1,  H/R  = 0.1  211 

18.  Frequency  Spectrum:  Case  I (cont),  H/RvO.l  212 

19.  Frequency  Spectrum:  Case  II,  H/R  = 0.1  213 

20.  Frequency  Spectrum;  Case  11  (cont),  H/R  = 0.1  214 

21.  Longitudinal-Shear  Vibration  Distributions  215 

22.  Radial  Plane-Strain  Vibration  Distributions  216 

23.  Clrcumfercntial-Shoar  Vibration  Distributions  21? 

vl 


KEY  TO  SYMBOLS 


H)OOgh  the  usage  of  symbols  in  this  work  Is,  In  the 
main,  consistent  with  modem  engineering  usage,  and  though 
each  symbol,  when  Introduced,  is  defined,  the  following 
listing  should  prove  of  use  when  Che  reader  is  tempted  to 
begin  his  study  in  the  middle  of  Che  text,  rather  than 
Che  beginning.  All  symbols  listed  below  may  occur  unsu- 
perscripted  (in  Chapters  I through  IV)  and  may  in  subsequent 
chapters  bear  the  superscript  "(1)",  where  1 = 1,  2,  or  3. 
This  marking  denotes  the  material  layer  of  the  three  layered 
shell  to  which  the  symbol  la,  at  that  moment,  ascribed, 
where  “(I)"  is  for  Che  innermost  layer,  and  "(3)"  is  for 
the  outer. 

tij  the  1,J  component  of  stress  in  an  orthog- 

onal, cartesian  coordinate  frame 

e-.  the  i,j  component  of  strain  in  an  orthog- 

onal,  cartesian  coordinate  frame 

Ft  the  I component  of  body  force 

Ciju  fourth  order  tensor  of  elastic  stiffnesses 

for  a general  anisotropy 

Tj  the  i component  of  stress  in  contracted 

notation,  defined  by  (2:4) 

the  i component  of  strain  in  contracted 
notation,  defined  by  (2:4) 

sii  ^ijkl  >^^^B^^hed  corresponding  to  the  con- 

tracted notation  of  (2:4) 

vil 


^ .'f.N'i 


elastic  compliances,  defined  by  equation 
(2:9) 

the  1 ccCTponent  of  displacement 
dens icy 

strain  energy  density,  defined  initially 
on  page  11 

Lame  isotropic  elastic  constants 
Poisson's  ratio,  Voung's  modulus 
cylindrical  coordinates 
cartesian  coordinates 

nondlraensional  radial  coordinate,  defined 
by  (5:22) 

wave  operator,  defined  by  (3:3) 

Laplace  operator 

Laplace  operator  in  r and  e , defined  by 
(3: 11a) 

various  potential  functions  used  generally 
in  Section  1 of  Chapter  III 

number  of  circumferential  waves,  defined 
initially  by  (3:12) 

wave  number,  defined  following  (3:12) 
wavelength,  defined  following  (3:12) 
natural  frequency,  defined  following  (3:12) 
phase  velocity,  defined  following  (3:12) 

wave  speeds  in  infinite  Isotropic  media, 
defined  by  (3:4b) 

group  velocity,  defined  by  (5:39) 
nondimensional  frequency,  defined  by  (5:2) 

frequency  of  first  antisymmetric  thickness 
shear  mode,  defined  by  (5:1) 


nondimensional 


number,  defined  by  (5:3) 


Jn,Y„ 

Ai,&i 


3".>n 
46tU  )1 


wave  velocities,  defined  by  (5:4)  and  (5:6) 

nondlmenslonal  wave  velocity  ratios,  defined 
by  (5:5)  and  (5:7) 

...  ,6  nondlFensional  elastic  con- 
stants, defined  by  (5:10) 

radial  wave  numbers,  defined  by  (3:17b) 
and  (3:l7c) 

factor  of  the  radial  wave  numbers  p,  g,  s, 
respectively,  defined  by  (4:3),  (4:37),  and 
(4:38) 

1 = 1, 2, 3 corresponding  to  the  shell  layers, 
sign  factors  defined  by  (4:2),  (4:34),  and 


nondlmenslonal  radial  wave  numbers,  defined 
by  (5:12)  and  (5:18) 


1 = 1,2,3  corresponding  to  critical  values 
of  the  given  ratio,  defined  In  (4:12)  and 
(4:13),  for  each  layer  of  the  shell 


nondlmenslonal  complex  zone  line  numbers, 
defined  by  (5:15) 


regular  Bessel  functions  of  order  n 
modified  Bessel  functions  of  order  n 


regular  or  modified  Bessel  function  of  the 
first  kind,  of  order  n,  defined  by  (4:40) 

regular  or  modified  Bessel  function  of  the 
second  kind,  of  order  n,  defined  by  (4:40) 

i = 1,2,3  corresponding  to  specific  displace- 
ment potentials , Integration  constants  In 
the  general  solution 

n=0  or  1,  special  representation  of  plate 
solution,  defined  by  (8:7)  and  (8:9) 

used  to  denote  the  determinate  of  ( ) 


lx 


sed  to  denote  special  de- 


1 an^inteser,  US 


R 

( )• 
R.fte 
I,  Im 


elements  of  the  frequency  equation,  defined 
initially  in  (5:37) 

complex  elements  of  the  frequency  equation, 
defined  by  (5:50) 


slope  of  a ray  emanating  from  origin  of  the 
A'Y  plane,  defined  by  (7:1) 

slope  of  ray  on  frequency  spectrum  corre- 
sponding to  Rayleigh  wave  velocity,  defined 
by  (7:12) 

ratio  of  d^^^'s  of  adjacent  layers,  defined 
by  (7:14) 


angular  argument  of  asymptotic  expressions 
for  regular  Bessel  functions  of  large  argu- 
ment, defined  by  (7:27).  Here  "(1)"  refers 
to  n^aterial  layer  of  the  shell 


shell  layer  thickness,  defined  in  Figure  5 


mean  radius  of  the  middle  layer  of  the  shell, 
defined  in  Figure  5 


complex  conjugate  designation 
denote  real  part  of  complex  argument 
denote  imaginary  part  of  complex  argument 
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The  propagation  of  hannonic  waves  in  infinitely  long, 
right  circular  cylindrical  shells  and  rods  is  treated  on 
the  basis  of  the  three-dimensional  theory  of  elasticity. 

The  shell  or  rod  is  assumed  made  of  three  concentric, 
transversely  isotropic  cylinders,  each  of  different  mater- 
ial, bonded  perfectly  at  their  interfaces.  The  displace- 
ment field  in  each  cylinder  is  represented  In  terms  of 
potential  functions  and  these  ere  subsequently  used  to  un- 
couple and  solve  the  Nevier  equations  of  motion.  The 
boundary  and  Interface  conditions  of  the  cylinder  are  used 
to  specify  the  integration  constants  of  the  solution.  The 
condition  of  nontrivial  solutions  for  these  constants  pro- 
duces an  18X18  determinant,  representing  a transcendental 
relationship  between  the  frequency,  wavelength  of  the  tra- 
veling waves,  and  the  niunber  of  circumferential  waves. 


This  frequency  equation  yields  ell  possible  modes  of  props* 
gatlon  In  the  cylinder.  The  cases  of  Infinite  wavelength, 
axlsymmetrlc  waves,  and  axisymmetrle  waves  of  Infinite 
wavelength  are  treated  for  traction  free  and  rigidly  clamped 
cylinders  and  rods.  The  frequency  equation  is  developed  for 
the  case  of  isotropic  layers.  The  case  of  very  thin  outer 
layers,  In  the  framework  of  the  exact  solution,  is  con* 
sidered.  An  asymptotic  analysis  of  the  frequency  equation 
is  conducted  for  the  case  of  very  short  wavelength  travel* 
ing  waves.  The  solution  for  the  cylinder  is  treated  in  the 
limit  of  infinite  radius  and  the  corresponding  solution  and 
concomitant  wave  equations  are  developed  for  a transversely 
isotropic,  three*layered  plate.  Finally,  a numerical  analy* 
sis,  based  on  a computer  program  capable  of  every  aspect  of 
the  above  problem,  is  presented  for  axisyiranetric  waves  in 
isotropic  layered  cylinders. 


CWAPTER  I 


the  iitedunlcs  of  wave  motion  of  all  types  has  long 
been  a topic  of  interest  to  physicists  and  scientists.  At 
the  macroscopic  level  perhaps  tbs  theories  are  furthest 
advanced  in  connection  with  the  clectrornsgnetic  and  elastic 
field  theories.  In  the  elastic  domain  the  desire  for  knowl- 
edge is  far  surpassed  by  the  need  for  practical,  applicable 
results.  Ue  can  not  know  enough  about  wave  motion  to  pro- 
gress beyond  application  of  this  knowledge  in  the  engineer- 
ing world  of  our  age. 

To  obtain  a complete  understanding  of  the  wave  patterns 
and  their  ensuing  stress  and  strain  fields  in  any  elastic 
body,  an  investigation  must  be  phrased  in  terms  of  a general 
elastic  body  - representative  of  the  particular  case  of  in- 
terest. Moreover,  a source  must  be  hypothesized  in  the  body, 
exciting  a generality  of  local  time  varying  strain  fields. 
Prom  this,  the  formulation  of  the  problem  must  be  framed  in 
an  acceptable  mathematical  model.  Appropriate  boxindacy  con- 
ditions and  differential  equations  arise  in  this  process 
which,  with  Che  exception  of  almost  trivial  eases,  resist  the 
most  persistent  efforts  of  solution. 

Humerous  avenues  of  escape  exist  from  this  dilenina.  One 
may  use  a more  crude  sachematical  model,  or  barring  this  (but 
which  is  almost  the  same),  may  linearize  and  otherwise  modify 


the  equations  of  the  problem  to  render  them  soluble.  An- 
other alternative  Is  to  modify  the  elastic  body,  specifying, 
for  instance,  a semi-infinite  medium,  an  infinite  plate,  a 
cylinder  cf  infinite  axial  length,  etc.  For  such  a specifi- 
cation some  of  the  boundary  conditions  are  simplified  and 
the  work  involved  in  manipulating  the  rigorous  equations 
may  be  rendered  more  tractable.  Further  reduction  in  the 
complications  of  the  general  problem  is  obtained  if  initial 
conditions  for  the  time  dependent  part  of  the  solution  are 
in  some  nay  obviated. 

Hltb  the  source  located  at  infinity  in  the  axial  direc- 
tion of  an  infinite  eUstic  cylinder,  it  is  clear  that  all 
those  forms  of  veve  motion  uhose  amplitudes  are  attenuated 
with  distance  will  not  be  observed  unless  the  patterns  of 
reflections  in  some  way  reinforce  each  other.  If  wave  mo- 
tion is  observed,  it  is  then  natural  to  assume  that  the  am- 
plitudes of  passing  waves  do  not  change;  that  is,  it  is  of 
little  consequence  at  whet  point  along  the  elastic  body  one 
picks  to  observe  the  wave  motion  - as  long  as  the  source  is 
an  infinity  away.  It  Is  also  convenient  to  pick  a source 
which  will  generate  infinite  trains  of  sinusoidally  varying 
traveling  waves.  This  Is  hardly  the  compromise  It  seems,  as 
can  be  seen  from  [l,2]]^.  The  method  of  Fourier  analysis,  us- 
ing these  harmonic  wave  results,  makes  it  possible  to  study 
more  complicated  wave  shapes.  Thus,  tdien  the  elastic 


^Numbers  in  square  brackets  refer  to  references  listed 
in  the  Bibliography. 


properties  ere  described  for  Che  body  in  question,  and  uhcn 
the  boundary  conditions  on  the  bounded  surfaces  are  decided, 
there  ceraalns  only  Che  Cask  of  choosing  a machemaclcsl  inodel 
which  will  yield  to  some  method  of  solution. 

Much  work  was  done  In  the  early  1800's  by  Poisson, 
Navler,  Stokes,  Rayleigh,  and  I^agrange  ^3,4,5,  and  6^. 

These  early  efforts  were  based  on  a finite  number  of  modes 
of  wave  motion  in  semi-infinite  bodies,  plates,  and  bars. 
These  first  "approxlinate''  solutions  were  valid  only  in  the 
low  frequency,  long  wavelength  range.  Solutions,  based  on 
linear  elastic  theory,  of  wave  problems  in  infinite  plates 
by  Lamb  [7]  and  infinitely  long  cylinders  by  Pochhammer 
[83  and  Chree  [93  were  presented  at  the  turn  of  the  cen- 
tury, Because  of  Che  complexity  of  Che  resulting  fre- 
quency equations,  they  were  not  evaluated  numerically 
until  recently.  However,  their  numerical  evaluation  is 
difficult  even  using  electronic  computers,  and  for  this 
reason  several  modem  approximate  theories  have  been  es- 
tablished, Mirsl^  and  Herrmann  [l0,ll3  produced  a five 
mode  Tlmoshenko-Cype  theory  containing  the  usual  membrane 
and  bending  effects  plus  the  influence  of  transverse  shear 
deformation  and  rotatory  Inertia,  A three  mode  theory,  de- 
veloped by  KcKlven,  Shah,  and  Sackman  [12,133,  containing 
the  framework  for  inclusion  of  higher  modes,  utilised  cor- 
rection factors  which  matched  the  approximate  mode  curves 
to  the  exact  theory  mode  curves  at  two  places  along  each 
exact  curve.  This  theoty,  which  is  for  hollow  cylinders,  is 


4 


an  extension  of  a rod  theory  developed  by  Klndlln  and  Mc- 
Niven  [l4],  Kaghdi  and  Cooper  [15,16]  provided  a three  mode 
approximate  theory  for  thin  shells , Two  sets  of  governing 
differential  equations  were  presented  based  generally  on  the 
inclusion  and  exclusion,  respectively,  of  small  magnitude 
terms  involving  the  shell  thickness  and  middle  surface  strains. 
The  first  set  of  these  equations  reduces  to  tove's  first  ap- 
proximation [17]  and  Che  second  set  to  Donnell's  equations 
[I6],  when  transverse  sheer  deformation  and  rotatory  Inertia 
are  neglected. 

The  advent  of  computers  has  brought  about  a thorough 
examination  of  the  Pochhammer-Chroe  frequency  equation  (that 
is,  the  exact  solution).  Mindlin,  Onoe,  and  McNlven  [19] 
and  Klndlln  and  Pao  [20]  have  fully  described  the  behavior 
of  the  branches  of  the  frequency  equation  for  real,  Imagi- 
nary, and  complex  propagation  constant  (wave  number),  Gaels 
[21,22]  extended  the  exact  theory  for  traveling  waves  and 
vibrations  to  thick  walled  hollow  cylinders , and  Greenspon 
[23]  independently  presented  the  vibration  theory  for 

The  dispersive  characteristics  of  axially  symmetric 
waves  in  infinitely  long,  two  layered,  elastic  rods  for  the 
ease  of  exlsymmetric  nontorsional  waves  were  investigated  on 
the  basis  of  the  linear  theozy  of  elasticity  by  HcNlven, 
Sackman,  and  Shah  [24]  and  of  axisymmetric  torsional  waves 
by  Armenakas  [25],  Baltrukonis,  et  a^.  [26]  solved  the 
problem  of  an  elastic  cylinder  rigidly  encased,  end  later 


treated  longitudinal  shear  vibrations  in  shells  composed  of 
two  elastic  layers  [27]  and  plana  strain  vibrations  in  com- 
posite shells  with  an  elastic  casing  and  incompresalve  inner 
layer  [28].  Armenakas  [29]  presented  a unified  theory  for 
aymmetrie  and  antisynnietric  harmonic  waves  in  two  layered 
cylinders  with  traction  free  surfaces.  Amonalcas  and  Keck 
[30]  presented  the  asymptotic  behavior,  for  very  large  pro- 
pagation constant,  of  Che  frequency  equation  for  traveling 
harmonic  waves  in  composite  two  layered,  traction  free, 
elastic  rods.  Jones  and  Whittier  [3l]  have  presented  stress 
and  displacement  field  distributions  on  the  basis  of  linear 
elastic  theory  for  the  case  of  exisymmetrie  wave  propagation 
in  two  layered,  traction  free,  elastic  cylinders. 

The  development  of  the  theory  - as  applied  to  anisotrop- 
ic elastic  bodies  - has  received  little  attention  prior  to 
recent  times.  This  was  partly  due  to  the  increased  diffi- 
culty of  Che  mathematical  aspects  of  the  problem,  and  partly 
because  there  was  little  need,  other  Chan  pedagogical,  for 
such  developments. 

The  development  of  the  exact  linear  theoiy  as  applied 
CO  tubes  and  cylinders  of  anisotropic  materials  was  Initi- 
ated by  Chree  [32]  in  1890  with  a discussion  of  the  propa- 
gation of  harmonic  waves  in  transversely  isotropic^  rods  of 
infinite  axial  length.  His  solution,  however,  was  based  on 

^Transverse  isotropy  is  characterized  by  any  material 
which  possesses  a plane  in  which  the  elastic  properties  ace 
the  seme  in  all  directions,  and  Che  elastic  properties  in 
the  normal  direction  to  this  plane  ace  different. 


approximate  methods  and  has  little  value  outside  the  low 
frequency,  long  wavelength  region.  In  1954  Morse  Q33]  de> 
veloped  the  exact  solution  for  compres s iona I waves  travel* 

Ing  In  solid  rods  of  a transversely  isotropic  material.  In 
the  1960's  Mirsky  C34,3S]  produced  an  exact  solution  for  axl- 
synmetrie  vibrationa  in  orthotropic^  cylinders . This  presen- 
tation utilised  Che  method  of  Frobenius  to  solve  the  equations 
of  motion  and,  consequently,  represented  a significant  de- 
parture from  the  usual  method  of  solution  of  problems  of 
this  kind  by  use  of  potentials.  Subsequently,  Mirsky  C36|] 
presented  an  approximate  theory  for  thick  walled  orthotrop- 
ic cylinders  of  the  same  nature  and  degree  of  accuracy  as 
the  Mlrsky-Kercmann  isotropic  cylinder  approximate  theory. 

Two  years  following  the  vibration  solution  by  Mirsky,  Nowln- 
skl  [37J  published  a paper  giving  the  exact  solution  for 
traveling  waves  in  orthotropic  rods.  He  also  used  the  method 
of  Frobenius  to  deal  with  the  equations  of  motion.  Kirsl^ 
published  a paper  [38]  in  1965  dealing  with  the  completeness 
of  solutions  for  anisotropic  elastic  problems  when  the  dis- 
placements are  represented  in  terms  of  potential  functions, 
and  a two  part  paper  [39,40]  providing  Che  solution,  in 
terms  of  displacement  potential  functions,  for  traveling 
waves  in  transversely  isotropic  rods  end  cylinders. 

The  exact  solution  of  Che  traveling  wave  problem  in 


^Orthocropy  is  the  name  given  a irmcerlal  possessing  the 
property  of  three  mutually  orthogonal  planes  of  elastic  sym- 
metry, Such  a material  Is  characterized  by  nine  Independent 
elastic  constants. 


composite  anisotropic  cods  and  cylinders,  that  is,  in  rods 
and  cylinders  with  two  or  more  bonded  layers  each  of  a dif- 
ferent inatcrial,  has  not  been  established,  Vu  [4l])  has  pre- 
sented an  approximate  theory  for  vibrations  in  isotropic  sand- 
wich elastic  ^lindrieal  shells  and  Chu  [42]  for  honeycomb 
core  sandwich  cylinders.  Both  of  these  theories  include  the 
effects  of  rotatory  inertia  and  transverse  shear  deformation. 
The  validity  of  such  theories,  however,  must  be  determined 
experimentally  or  by  comparison  with  solutions  established 
on  the  basis  of  the  three-dimensional  theory  of  elasticity. 

Use  of  monocrystalline  materials  (for  example,  qtiarts 
crystals  in  electronics)  has  created  an  additional  need  for 
wave  theory  in  the  ultrasonic  frequency  range  applicable  to 
single  crystal  plates  of  various  anisotropies.  Furthermore, 
plastics,  fiber  reinforced  shells  with  hexagonal  elastic 
symmetry  or  near  hexagonal  symnetry  [43],  and  recently  in- 
troduced engineering  materials  in  the  aerospace  industry, 
such  as  graphite  end  Berelllum,  all  present  differing  re- 
quirements for  establishing  the  transverse  isotropic  theory. 
Moreover,  although  polycrystalline  materials  exhibit  little 
tendency  towerd  anisotropic  properties  under  ordinary  condi- 
tions [44],  when  rolled  into  thin  sheets  metals  such  as 
steel  exhibit  a mild  transverse  isotropy  [45],  Sheet 
Berelllum  [46]  likewise  exhibits  this  property,  but  to  a more 
distinct  degree. 

The  theory  contained  herein,  which  is  established  for 
three  layered  cylinders  and  plates,  aside  from  being  a 


contribution  to  the  exact  elastic  theory,  will  provide  the 
required  analytical  basis  for  nany  of  the  problem  areas 
cited  above. 

In  this  investigation  the  propagation  of  an  infinite 
train  of  harmonic  vaves  in  infinitely  long,  right  circular, 
cylindrical  shells  and  rods  will  be  treated  on  the  basis  of 
the  throe-dimensional  theory  of  elasticity.  The  shell  or 
rod  will  be  assumed  made  of  three  concentric,  transversely 
isotropic  cylinders,  each  of  different  materials,  bonded 
perfectly  at  their  interfaces.  The  displacement  field  in  each 
cylinder  will  bo  represented  in  torms  of  potential  functions 
and  these  will  subsequently  be  used  to  uncouple  end  solve  the 
Wavier  equations  of  motion  boundary  and  interface 

conditions  of  the  cylinder  will  be  used  to  specify  the  inte- 
gration constants  of  the  solution.  Nontrivial  solutions  for 
these  constants  produce  an  18X18  determinant  which  is  var- 
iously referred  to  as  the  frequency  or  dispersion  equation. 
This  equation,  representing  a transcendental  relationship 
between  the  frequency  and  wavelength  of  the  traveling  waves , 
and  number  of  clxrcumf erential  waves,  yields  all  possible  modes 
of  propagation  in  the  Q’linder.  It  will  be  shown,  as  in  the 
case  of  simple  Isotropic  shells,  that  the  fTrequency  equation 
degenerates  into  simpler  cases  for  axially  syncietric  waves 
(torsional  and  nontorslonol)  and  for  waves  having  infinite 
axial  wavelength  (longitudinal  shear  and  plane  strain). 
Asymptotic  expressions  will  be  presented  for  the  behavior  of 
the  frequency  equation  when  the  wave  number  is  very  large. 


Ihe  limiting  case,  as  the  radius  of  the  cylinder  becomes  in- 
finite, will  be  treated,  end  the  frequency  equation  for 
three  layered  plates  established, 

A computer  program  will  be  written  for  the  general 
problem  outlined  above.  Because  of  limited  availability  of 
machine  time  numerical  results  will  be  presented  only  for 
the  ease  of  axisyirmetric  waves  traveling  in  shells  composed 
of  isotropic  layers. 


CHAPTER  It 


1.  Basle  Eouet Ions 

The  principle  of  consocvatlon  of  linear  monenrum  of  an 
element  dV  contained  In  a body  D bounded  by  the  surface  S 
gives  [47.A8] 

■Cij.l  (2:1) 

vhere  ^ ie  the  mass  density  per  unit  volume  of  D,  ty  ate 
the  components  of  the  stress  tensor  from  the  linear  theory 
of  deformable  continue,  and  Fj  are  the  components  of  the 
body  forces.  Differentiation  with  respect  to  time  and  apace 
have  been  indicated,  respectively,  by  the  “dot"  and  "eonroa" 
conventions.  A repeated  subscript  indicates  summation 
over  all  possible  values  of  the  Indices  (1,J  ° 1,2,3).  The 
components  of  displacement  Uj , are  related  to  the  compo- 
nents of  strain  of  a linearly  deformable  continuum  by  the 
relation 

+%l)  • <2:2) 

tensor  by  the  following  constitutive  relation: 

‘ = (2:3) 

Since  the  stress  and  strain  tensors  are  symmetric.  It  fol- 


11 


This  expression  reduces  the  nunber  of  independent  constants 
in  the  tensor,  from  81  to  36. 

If  the  etrsin  components  in  a unit  element  of  tho  body 
D are  Increased  from  to  tho  «orh  done  is 

M - 'tvj'Jeij  ) 
and,  using  (2:3),  this  gives 

dW  ' , 


Alternative 


, interchanging 


duiraty  subscriptE 


Therefore,  since  the  order  of  differentiation  is  im- 
material Ci73,  it  follows  that 
^IJW  ' Cuij  - 

This  last  relation  reduces  the  number  of  independent  elastic 
constants  from  36  to  21,  which  Is  the  namber  of  constants 


12 


for  a general  anisotropy. 

To  avoid  dealing  ulth  double  suns  and  large  groups  of 
subscripts  the  folloulng  notation  Is  Introduced 


In  teems  of  (2:4),  eqtjatlons  (2:3)  become 


= CijCj-  . (i,j  " 4).  (2:5) 

The  constitutive  relation  (2:5),  uhich  Is  the  general- 
ized Hooke's  Law  of  the  linear  small  deformation  elastic 
theory,  can  only  be  further  elrnplif^od.  If  the  mechanical 
response  of  the  material  considered  is  symsetrle  with  respect 
to  axes  or  planes  In  the  solid  etc.  If  the  elas- 

tic properties  are  symmetric  with  respect  to  three  mutually 
orthogonal  planes,  the  number  of  Independent  elastic  constants 
is  reduced  to  nine  and  (2:5)  yields  the  six  equations  of  or- 
thotropy 


(2:6) 


If  the  elastic  properties  of  the  solid  are  symmetric 
with  regard  to  an  axis  of  symmetry,  z,  the  nine  Independent 
constants  in  (2:6)  are  reduced  to  five  and  the  solid  is 
termed  transversely  Isotropic.  For  this  case  equation  (2:5), 
written  in  expanded  form,  given 


(2:7) 


The  fact  that  the  strain  energy  of  an  elastic  solid  in- 
creases through  positive  values,  from  zero  in  the  strain 
free  state,  Imposes  restrictions  on  the  range  of  values  to 
which  Che  elastic  constants  may  obtain.  These  limitations 
are  brought  out  most  easily  by  noting  that  the  strain  ener- 
gy density  function  is  a positive  definite  quadratic  form. 
Hence, in  the  usual  manner  of  establishing  the  criteria  for 
positive  definiteness  of  a quadratic  form  [49^,  each  of  the 
inequalities  below,  for  p = 0 to  5,  must  hold 


Hence,  for  t 


3 isotropy  in  particular 

1 >0. 


are  the  limitations  liiposed  on  the  elastic  constants  by 
positive  definite  strain  energy  density. 

The  elastic  stiffnesses  in  fomiulas  (2:6)  and  (2:7) 
can  be  related  to  engineering  moduli  and  Poisson's  ratios  in 
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the  followins  fashion.  The  determinant  of  the  elastic  stiff- 
nesses physically  can  not  be  zero.  Thus,  equations  (2:5)  can 
be  inverted  and  the  strains  expressed  as  linear  combinations 
of  the  stresses 


' '.i.  •.'*1,  (a) 

where  (2:9) 

datlS(il  drtlt-,jl  = 1.  (b) 

The  are  termed  elastic  compliances  and  bear  the 
follouing  physical  significance  [44]: 

a)  s^q  (q  = l,2,  or  3)  relates  an  extensional  stress 
to  an  extensional  strain  both  in  the  sane  direction, 

b)  s^^  (qcr;  q,r  = l,2,  or  3)  relates  an  extensional 
strain  to  a perpendicular  extensional  stress,  and 

c)  8^^  (related  to  (q  = l,2,  or  3;  r = 4,5,  or  6) 

relates  an  extensional  strain  to  a shear  stress  and 
vice  versa,  both  in  Che  same  plane, 

6^^  (related  to  s^j^j^j)  (q»l,2,  or  3;  ro4,5,  or  6) 
relates  an  extensional  strain  to  a shear  stress  in  a 
perpendicular  plane  and  vice  versa: 

d)  s^q  (q  = 4,5,  or  6)  relates  a shear  strain  to  a sliear 
stress  in  the  same  plane,  and 

5,1.1  ° ' ^‘*'1  ■’ 

e)  8^^  (q=r;  q,r  = 4,S,  or  6)  relates  a shear  strain  to 
a shear  stress  in  a perpendicular  plane. 
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Because  of  the  analytical  similarity  of  equations  (2:9) 
and  (2:5),  the  requirements  used  to  reduce  (2:5)  to  the  or- 
thotropic or  transverse  Isotropic  cases  will  reduce  (2:9)  to 
the  same  configuration,  except,  taking  note  of  (2:4), 

The  c^j  can  be  related  to  engineering  parameters  E , 
only  through  the  compliances.  From  (2:9b)  the  relation 

Cij  ‘ |Sij  | , (2:10) 

where  is  the  minor  determinant  associated  with  Sj^j,  is 

easily  obtained.  Consequently,  the  following  representation 
for  the  Cj^j  in  terms  of  the  Sj^j  and  the  engineering  constants 

, . . 

. 5,>  -Sii6»s  


Similar  relations  are  easily  obtained  for  orthotropy. 
The  ultimate  reduction  in  the  number  of  independent. 
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Lastlc  constants  Is  obtained  when  the  mechanical  properties 
E the  solid  are  the  same  in  any  direction.  This  is  isotropy 
id  is  characterized  by  tr«>  independent  constants.  The  stress- 
train  relations  are 


r. (k-4,5,6). 

The  elastic  constants  in  (2:8)  may  be  associated  with 
the  Lame  coefficients  ^ and  h,  or  Young's  modulus  E and 
Poisson's  ratio  v through  Che  following  schema: 

C„  e.(i-w)/(u»Xi-lw) 

C.,  A 

(c„-C„)A  I*  E/tdtP) 

(c„*icj(c,„ -c,iV(c„*0  E 


2.  Disolncemant  Eouat ions  of  Motion 

Since  Che  development  required  in  this  work  must  ul- 
timately be  phrased  in  cylindrical  coordinates,  all  deriva- 
tions will  be  carried  forward  in  the  coordinate  frame  r,  e, 
z,  where  z corresponds  to  the  axial  coordinate  and  r and  6 
are  defined  in  the  usual  way. 

In  cylindrical  coordinates  Che  equations  of  motion  may 


be  written  explicitly 


|I- ^l|T'e^|5r.  , 

ti35i.  +?3sb  . , 

Sr  r r ga  "SI  ‘ 

3trt  4.  Lt  . * i2S«  1-  i?M  = - 

tr  r * r 99  9*  'gl*- 

Because  the  stress-strain  equations  for  anisotropy  do 
not  yield  to  s convenient  Censorial  notation  as  in  the  case 
of  isotropy,  the  development  uill  be  continued  in  en  expand 
form  to  arrive  at  working  equations. 

The  geometric  strain-displacement  equations  are 
t„  se,  = , 

= e = Ur  , 


By  substituting  equations  (2:13)  i 
lowing  stress-displacement  equations  at 


* T,  = C„-aUr 


(2:14) 


Substituting  equations  (2:14)  into  the 
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equations 


■ i|s,. . (.) 

•"“(!;.• -tit-  -Ff)* 

•f:g’  -'"’S'  ■ ti?-  . (b)  (2:15) 

iSi-'  M 


'‘■‘•‘“'iS:  + "•“> 

‘"■S.'S**  * & -'"(St”  *?5t-n“*)-S?S''  "■> 


(c)  (2:l6J 
(cont^ 


t,  * c»,  1 ■^‘wi  - - i3av\  , 

”ir-  ' 3e»-  fir) 

Finally,  for  isotropy,  the  displacement  equations  of 
motion  tahe  the  compact  form 

i(o.-oji»,.ii  (‘■'.‘-.’■I.  «:i» 


CHAPTER  III 


!•  ^ the  Uae  of  DisolflceinanC  Potential  Functions  to  Solve 

the  Eauations  of  Morion 

The  classical  general  solution  of  the  isotropic  equa- 
tions of  motion  was  established  by  Poisson  [50]  (1829).  He 
showed  that  every  (sufficiently  regular)  solution  of 

admits  the  representation 

U:(X,1)  =u'.(-x,t) (a) 

(3:2) 

; U;,.  ' 0,  (b) 

where  t^(x,t}  and  u'^(x,t)  satisfy  Che  wave  equations 

□|cf  = 0,  DjU-=o  (3:3) 

°n  = V*  -f  , (n=  l.l),  (a) 

" (3:4) 


and  V Is  l:he  Laplace  operator. 

Poisson  thus  established  Che  fundamental  fact  chat  a 
complete  solution  of  (3:1)  can  be  expressed  as  a superpo- 
sition of  an  irrotational  and  equivolureinal  notion. 

The  first  appearance  of  the  displacement  field  described 
in  teems  of  potentials,  that  is,  the  Scolces-Kelmholtz  reso- 
lution, 

20 


21 


■ (3:5) 

UBS  apparently  due  to  Lama  ^513  (1832),  who  noted  that  every 
vector  of  the  form  (3:5)  satisfies  (3:1)  if 

O^Cf  " o.  '0  . (3:6) 

Notice  that  (3:5)  end  (3:2)  differ  only  in  the  definition  of 
the  rotational  portion  of  the  displacement  field. 

Clebsch  CS2^  (1863)  stated,  and  later  Somlgliana  [53^ 
(1892),  Tedone  [5ft]  (1897),  and  Dohom  [5S]  (1898)  independ- 
ently proved,  that  any  solution  of  the  equations  of  motion 
(3:1)  admits  the  representation  (3:5)  and  (3:6).  These 
facts,  along  with  Duhem's  version  of  the  completeness  theorem, 
are  given  in  a paper  by  Sternberg  [56].  Duhem's  completeness 
theorem  vithout  proof  is  included  here: 

Let  u^(x^,X2,X2,t)  be  a particular  solution  of 

(2:17) 

in  a region  of  space  D and  for  t^<  t<t2.  Then  there 
exists  a scalar  function  <f(x^,X2,X2,t)  and  a vector 
function  V^(Xj^,X2,X3,t)  such  that  Uj^(x^,X2,Xj,t)  is 
represented  by  (3:5),  with 

’t'i.i*®.  (3:7) 

and  and  satisfy 

D^<f  = °,  Dj  ■('i  * o , (3:3) 

where 

U = (a)  (3:ft) 
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(b)  p:4)^ 


The  operator  In  (3:4)  is  referred  to  as  the  wave  opera- 
tor, and  the  velocities  and  C2  in  (3:4)  are,  respectively, 
the  dilatational  and  equlvoluinlnal  wave  speeds  in  an  infinite, 
linearly  elastic,  hosMgeneous  isotropic  solid. 

Equations  (3:S)  and  (3:7),  taken  together,  indicate 
that  u.  nay  be  represented  in  terms  of  three  independent 
potential  functions : cp  and  two  of  the  It  is  to  be  noted 

that  Uj^  given  by  (3:5)  satisfies  (3:1)  provided  (3:3)  hold, 
even  if  i|/j^  is  not  solenoldal.  Consequently,  equation  (3:7) 
may  be  relaxed  to  the  form 


This  representation  has  been  used  in  many  of  the  contempo- 
rary investigations,  for  example,  [21,28,29,30],  etc.,  in 
establishing  the  complete  solution  of  (3:1)  when  u^  represent 
the  displacement  field  for  Infinite  trains  of  harmonic  waves. 
Shifting  attention  now  from  Isotropy  to  anisotropy,  it  is 
early  noted  that  the  wealth  of  formal  solutions,  so  prevalent 
In  the  isotropic  theory.  Is  almost  entirely  lacking.  Ap- 
parently only  the  ease  of  transverse  Isotropy  has  successfully 
been  treated  by  the  use  of  potentials.  Lekhnitsky  [57]  and 
Elliot  [58]  solved  the  displacement  equations  of  tocslonless, 
exisysimetric  equilibrium  for  several  Important  problems  in 
terms  of  a single  stress  function  satisfying  a fourth-order 
partial  differential  equation.  Eubanks  and  Sternberg  [59] 
developed  a completeness  proof  regarding  the  use  of  this 


ti.i  = 


(3:8) 
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potential  solution  . Mirsic}’  [38^  subsequently  gave  a com- 
pletensss  theorej-ti  extendins  tho  applicability  of  the  Le'chait- 
sky  potential  solution  to  include  Che  nonsynnetric  elasto. 
static  problem. 

Mlrsky  [38]  also  provided  a general  theorem  concerning 
the  solution  of  the  elastodynamic  problem  in  transversely 
isotropic  media.  Through  use  of  a sec  of  potential  functions 
defining  the  displacement  field  according  Co  Che  follot^ing 
equations 


* 13^  , 

3r  r 38  ’ 


(3:9) 


u 

^ Sm 

ho  showed  that  tho  equations  of  siotlon  for  a transversely  iso- 
tropic solid  are  satisfied  if,  in  turn,  U,  V,  and  W satisfy 
the  differencial  equations 


I V ‘8  ■ .V4.  3^1  T > 


(a) 

(3:10) 


AV  «AW  • 0 . 


(b) 


Here 


(a) 


- V*  -V  [(c„c„  -c^,  -2.c„c.,)/c„c*.]£,v,^- 

' * If;’  1^,  - 


(3:11) 

(b) 
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{cent)  (b)  (3:11) 

“ (cont) 

As  will  subsequently  be  shown  in  section  2 of  this  chap- 
ter, the  potentials  U,  V,  and  U used  above  by  Mirsky  follow 
closely  from  the  Stokcs-Helmholtz  resolution  of  the  vector 
displacement  field.  However,  the  third  of  (3:9)  represents  a 
minor  (and  nonsimplifying)  departure  from  the  usual  specifi- 
cation of  Che  Helmholtz  potentials. 

Any  attempt  to  factor  the  operator  A of  equation  (3:11) 
into  a product  of  two  second  order  operators  leads  to  a sys- 
tem of  algebraic  constraints  between  the  elastic  constants, 
which  apparently  are  only  satisfied  by  the  condition  of  isot- 
ropy. Hence,  any  solution  u.  of  (2:15),  in  a region  of  space 
D and  for  t^  C < Cj,  may  be  represented  by  (3:9)  provided 
V satisfies  (3:10a)  and  U and  V satisfy  (3:10b).  It  is  read- 
ily apparent  Chat  Che  wave  operators  no  longer  play  a role  in 
the  general  solutions . 

Equations  (3:11)  are  indicative  of  the  mathematical  for- 
nldlty  of  the  elastodynamlc  problem  in  anisotropic  linear 
elasticity.  No  formally  acknowledged  potential  solution 
(Galerkin,  Lekhnitsky,  etc.)  has  been  successfully  applied  to 
even  so  simple  an  anisotropy  ss  Chat  under  discussion  in  this 
paper.  By  successful  application  is  meant,  of  course,  a ren- 
dering of  Che  elastodynamlc  problem  into  a solvable  system  of 
equations,  capable  of  being  shown  to  provide  the  complete, 
general  solution  of  the  problem. 

In  applying  the  Scokes-Helmholcz  resolution  of  Che  vector 
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displacement  field  to  the  isotropic  equations  of  motion  for 
an  Infinite  body,  it  is  natural  to  obtain,  in  uncoupled  form, 
the  effects  of  dllatational  and  rotational  motion.  In  an 
infinite  isotropic  medium  the  displacement  vectors  associated 
vith  dllatational  and  rotational  waves  are  collnear  with,  or 
orthogonal  to,  the  wave  front  vector,  respectively.  However, 
in  a general  anisotropic  infinite  continuum,  the  three  dis> 
placement  vectors  associated  with  the  equivolumlnal  and 
dllatational  waves  have  no  such  orientation  end  the  displace' 
ment  vector,  which  is  collnear  with  the  wave  frost  vector,  is 
generally  comprised  of  both  equivolumlnal  and  dllatational 
displacement  componente  Furthermore,  the  two  equivolumi- 

nal  wave  velocities  are  generally  dissimilar  in  an  aniso- 
tropic body.  The  magnitudes  of  the  three  wave  velocities,  in 
addition,  generally  depend  on  the  direction  of  propagation 
within  the  continuum.  These  facts  imply  that  in  a general 
anisotropic  continuum  the  equivolumlnal  and  dllatational  com- 
ponents of  motion  will  be  coupled. 

When  solutions  of  the  equations  of  motion  are  sought  in 
cylindrical  coordinates  for  a transversely  isotropic  material, 
the  equations  are  as  complex  as  those  for  an  orchotropic  ma- 
terial, if  the  elastic  symmetry  axis  is  not  chosen  to  coin- 
cide with  the  z axis  of  cylindrical  coordinates.  No  substi- 
tution, in  terms  of  potential  functions,  has  been  found  for 
the  displacement  componente  which  renders  the  orthotropic  equa- 
tions soluble.  In  fact,  use  of  the  conventional  resolutions, 
for  example  (3:9),  for  the  displacements  results  in  a s 
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of  three  third  order  partial  differential  equations  with 
variable  coefficients,  each  involving  the  three  potentials. 
Due  to  the  variable  coefficients  in  the  equations  it  is  not 
possible  to  organize  the  three  eqiiations  into  useful  opera- 
tor forinulatlon,  since  the  operators  do  not  convnute.  Futher- 
more,  among  the  three  resulting  equations  there  is  no  degree 
of  similarity  to  aid  in  the  solution.  Consequently,  though 
it  is  conceptually  possible  to  manipulate  the  system  of  equa- 
tions to  obtain  a nev7  aet  involving  each  potential  function 
separately,  it  is  no  longer  possible  to  expoct  each  potential 
to  satisfy  the  same  higher  order  differential  equation  aa  is 
the  ease  if  the  differential  operators  commute  C61.62]- 

Frequently  simplifications  in  the  governing  differential 
equations  can  be  obtained  by  assuming  the  functional  fonn  in 
which  certain  variables  enter  the  formulation.  Additionally, 
^en  potential  function  formulations  do  not  aid  the  solution 
process , it  is  customary  to  seek  particular  solutions  of  the 
equations  of  motion  in  torras  of  the  displacements  directly. 
However,  except  for  the  simplest  of  anisotropies,  these  two 
reductions  in  complexity  are  of  little  aid.  This  is  easily 
llluBtrated  by  the  following  (classical)  particular  solution. 

U,  cosne  crttki  -tut). 

u.  ■ VCri  (3:12) 

U4  ••  Wlrl  tosnS  Sm(Wi  -lolV 

When  introduced  in  equations  (2:16).  Here  k = 2t/^  is  the 
wave  number  and  "h  the  viavelength;  cj  is  the  natural  frequency, 
se  velocity  of  the  wave  is  given  by  v = <J/k. 


The  phas 
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^ (c.i,  t,s>^  «•  (a) 

"“0-  ^ ^ ^ ^ -**'«)  ■?.  V - 

<3 

-{c^+c,^^-*  w = 0,  (b) 


essentially  all  the  complications  p 
1 Mitsky  [34,35]  and  Kowinski  [37]  a 


.'s  [63] 


ited.  If 


2S 


in  some  (not  necessarily  simply  connected)  region  D,  then 
g(r,6)  can  bo  represented  as 

g‘^B)  = g,(r,el  ^ g,(r,B).  (a) 

where  (3:14) 

. (b) 

and  is  given  in  (3:11a). 

Theoreni:  Let  (i  •=  r,®,z)  be  a solution  of  Che  equations 

of  motion  (2:15)  of  a transversely  isotropic  elastic  con- 
tinuum in  some  region  (not  necessarily  simply  connected)  of 
space  D and  for  c<  C2>  Then,  there  exist  scalar  func- 
tions Cf(r,6)exp(lkz  - iut)  and  f (r,s)exp(ikz  - luc)  such 
that  the  displacement  field  is  represented  as 

(3:15) 

Ui  - 7iplr,e)]eipl(V.i-wlV 

where  are  constants  with  respect  to  the  variables  r,6, 

and  z,  is  given  by  (3;lla),  and  ^ and  f satisfy,  respec- 

□^(  = 0.  (3:16) 

The  wave  operators,  D^(i  = 1,2,3),  are  defined  as 

(a)  (3:17) 
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1b  which  the  "+*'  and  Bigna  are  associated  with  and 
q^,  respectively,  and 

. ij,„.  . W (3:17) 

Proof:  According  to  the  Stokes-Hslnholtz  theorem  C64,653> 

every  sufficiently  smooth  vector  field  u,  nay  be  written  as 


U;  * Uji  + (3:18) 

Expanding  the  vector  (3:18)  into  components  yields 

(tm  "E* 


Vk 


+ A^e  ^ iV,.  -O’ 


(b) 


In  equation  (3:19b)  the  term  3V^3z  has  been  eliralnated 
by  taking  into  account  the  fact  that  the  vector  satisfy- 
ing Che  first  of  (3:18),  nay  contain  an  arbitrary  gradient 
quantity.  That  is 

V = Vd  . 

where  q,^^  J-a  chosen  equal  to  -bV^/bz.  Consequently,  the 
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divergence  of  V.  vanishes  if  (3:19b)  is  satisfied.  This 
may  be  accomplished  by  choosing 

Using  this  substitution,  end  subsequently  letting 

FCr.el  = IK'-.e)  + iWt(T,e), 

(3:20) 

G(r,e)  =Vj(r,e), 

equations  (3:19)  reduce  to 

»e'(-f|5  -|f  )e»piUa-ul).  (3:21) 

Uj  =(cUF  •*-  -o>0  I 

with  = W*’  - Vf  • 

Using  equations  (3:21)  in  (2:15)  yields 
F i (i5u.*^-CuU’-)F 

+ J:|e[c«V,^G  + (?■-*  -=«V-K]'  0,  (a) 

^l!fc,.c4F 'S*)]- 

- ] - 0.  (b)  (3:22) 

Appropriate  differentiation  of  the  first  two  of  (3:22) 
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F -(r’- 

[c,4^,^G.  *(tw‘  -c„v'-)&]  = U/'-.e), 


vJViiCr.el  = Vfh^(r,e)  =0. 

Substitution  of  (3:23a  and  b)  in  (3:22) 

h2(r,e)  such  that 

aj?,  ^ _ 1 'Si.,  - Sisi  . 

3r  r 3®  ’ t 36  "s." 

Equations  (3:22c)  and  (3;23a)  form  a set 
llns  P and  Z: 

L,r  - E =b,U,e). 


L,"  c„ 


L^* 

Lj«  Ct,j+2.c„')^f-  + (^u5^ 

L*-  C^vf- 
Cons oquently. 


i]]  = )',h,61.(a) 

(b)  (3:23) 

(c) 

relates  b^(r,e)  and 
(3:24) 

of  relations  coup- 
(e) 

(3:2S) 

(b) 
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(L,L^-i- L^.Li)F  ' L^K,  = (^w>- -CjjWMW,  (a) 

ana  (3:26) 

^ (l,L.  * LtLjz.  ' -L,Vi,  = -CjjV'^)K,  • (b) 


L,L,  + Lj^Uj  = C„C^^V^V,*’  » 

V,^  + (3:27) 


+ ( - cjjk'-)C^w'^  - =44 '■'■')  ■ 

Solutions  of  (3:26)  are 

F = f'"  + F'«  , 

Z = 2.“’  + E ; E « 7?Z./ik  , 

ana 

= h,Cr,e)/i5w»--C44V‘  , 

z“'  • -W,('-,e)/^uj^  -=44'J’  • 


(3:28) 


(3:29) 


(3:30) 


In  vleu  of  (3:27),  equation  (3:29)  n:ay  be  vritten  as  a 
factored  operator  if 

(v‘  * p-)(  vf-4  ciM  =^-^(l,L*  +L,Lj)  , 


Sl*=44p'«(^  =(Su>’-C„V’-K^u.^'C44k’-'). 
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Solving  the  above  for  and  yields : 

1 (3:31) 


Accordingly,  the  general  solutions  for  F*'’and  Z '' may 
be  represented  as 


2-"'* 


. ( v‘^'^q.'){pi,Z-2j  *0.  (c) 


(3:33) 

(b) 

3 find  a 


F -F.  +F‘“  , 

E. 

Equation  (3:33a)  uith  (3;32c)  is  used  in  (3;23a) 
representation  for  Z,  This  yields 

(3:34) 

Using  this  expression  for  Z,  and  (3:33a)  for  F,  in  (3:22c) 

0 = 


which,  using  (3:31),  vanishes  identically. 

Turning  now  to  (3:23b),  a solution  is  given  by 

G -•  6,  + G‘^ 


(a)  (3:35) 
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h,  fr.e^ 


A representation  for  the  displacement  field  is  obtained 
using  (3:33a),  (3:34),  and  (3:35)  in  (3:21)  such  that 

"'■®  * -5  3‘>p «»-»'). 


Using  equations  (3:24),  (3:30),  and  the  second  of 
(3:35b),  the  following  simplifications  are  made: 


Define 


5 

Ky-,0)  * &,(r,e), 


□fDtc?  = 0,  DjU 


U) 

(3:36) 

0>) 


(3:37) 


□I  = V^  + p*, 


(3:38) 
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V^  * s\ 

with  and  q’’  as  given  by  (3:31),  and  by  (3:35), 
Then 

Ue--  -^^]eTcpHwt-wt). 

and 


(3:39) 


Altefarlon  of  the  Reoras entatlon  for  Two  Special  Case 


There  are  two  specific  situations  that  need  to  be  men* 
tloned  which  might  be  thought  to  effect  the  generality  of 
Che  above  representation:  (1)  isotropy,  and  (2)  Cj^j  + c^^^O. 

If  Che  material  body  D is  isotropic,  then 
033-0^1,  044-^=66  and 

ff  . ,C|^  = 5^  (3:40) 

No  essential  difference  exists  between  this  condition 
and  Che  theorem  conditions,  hence 


Ub  = _|P  )etplCvt-<,>i), 


(a)  (3:42) 
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Q^<{,  = 0 


(b) 


( 


(c) 


The  second  case  must  be  alloued  since  the  positive 
definiteness  of  the  strain  oner®’  density  function  requires 
only  that  Cj^j^,  c^g,  c^^,  and  positive.  The 

constants  and  Cj^2  *’°*^*’  nesative  to  extents  which 
are  limited  by 

C„'=tt.>0  , (3:63) 


Consideration  of  the  situ 
treated  by  starting  with 
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illy 


(3:44) 


Equations  (3:64)  lead  to  (as  done  previously) 

□^F'h,(r,e),  vfW,  =0,  (a) 

OjGi  * h^(t,e),  vfhj=0,  (b)  {3;45) 

□ jF  t 0*2  =0.  (c) 

where  the  Dj^(i  = 1,2,3)  operators  are  as  previously  defined 
with 


(3:46) 
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Fu^th€fIl]o>^e , h^(r,e)  and  h2(r,d}  are  related  by  (3:24). 
Operating  on  (3:45c)  with  0^  end  noting  (3:4Sa), 


=-Dih|U,S)  = (3:47) 

since  the  □,(1  = 1, 2, 3}  operators  comnute. 

Consequently,  Z is  represented  by 

(a) 

with  (3:48) 

0^2-,=  atZ-i  “0, 

The  solution  of  (3:45b)  Is  given  as 

■vG"’  (a) 

with  (3:49) 

G'"  - , (b) 


Equations  (3:45a  and  c)  lead  directly  to  the  following 
representation  for  F In  terms  of  Z and  hj^(r,e): 

(3:50) 

Substitution  of  (3:48a),  (3:49a),  and  (3:50)  Into  (3:21) 
yields 


(3:51) 


= ikZj^espaWi-uit). 
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n representation  (3:39) 
tlons  (3:36),  In  the 


limit  as  (Ci3  + c^4)  e 


LCWg-tuil,  (a) 

in  (3:520)  vaniehes  a.  (0,3.0,,)  approaohes  oero.  Althoogh 
(3:46)  gives  the  proper  form  for  p»  when  (0,3  + 0,,)  is  exact- 

(*=13*°44>>  P’ 


' <'13*=44>  " 


eipi(kt-«t)-  (3:54) 
consequence  of  equations  (3:52a  and  b)  and  (3:54), 
r that  the  special  case  (c,3  + c,,)»0  is  Included 
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in  the  formulation  of  the  completeness  theorem,  for  If 

£,  Z-i  = D^tfi  = 0 , 

• “4, 

and 

C.  = ^ □j&'Djf'O, 

equations  (3:51)  become  just  (3:52a  and  b)  and  (3:54). 

As  a consequence  of  the  above  results  for  isotropy  and 
(Cj^j  + = 0,  Che  representation  given  in  the  theorem  of 

the  previous  section  is  con^letely  general, 

*.  Specification  of  <?(r.e)  and  f (r.9) 

Returning  to  the  general  representation  for  the  displace- 
ment field,  (3:39),  the  scalars  (f  and  f are  further  restricted 


(f(r,e)  • C(l(r)lcosT>e  -visinnel, 

(3:55) 

Ur,e)=  Uotsmne -V  tosne] . 

Kahing  Che  usual  specification  Chat  the  solution  be 
the  real  part  of  the  above,  and  f are  chosen  as 
!eCr,e,i,V)  = Cftr)c05ne  cosCVl-uil, 

(a) 

l(r,e.s,t)*  Krtemne  cosCM-wl), 


U,  = (cc.,40^''‘  - ^fW]cosne  tosCki,-«t). 

Ue  - " he  costVi-ut),  (b) 

Uj  - +c*4(?^-l.%Wj«bnesm(ia-wt)- 
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If,  alternatively,  the  Intaginary  pert  of  (3:55)  is  chosen 
for  Che  solution,  Chen 

cp(r,e,»,0  * 4it»r>e  Sif\(kh-Ut), 

^ (a) 

?u,e,e,f)=  ?Wc£>sn6  woUi-wt), 


Ug  =:  - ((t,i+c«)Qc(>Cr)  -V  ^'•''‘yosne  sinCkt.  (b) 

“i  ‘ ■ ^(^'1  '^1^'')  -^):?Cr1^5inn6  cos(ki-wt). 

Comparison  of  (3:56)  and  (3:57)  reveals  that  in  Che 
case  of  axisymmeCric  displacements  (n  = 0),  Che  torsional 
displaccmonc  field  (u  vu  =0)  is  the  imaginary  part  of 
(3:55)  and  depends  on  the  single  potential  f,  while  the 
non~torsional  field  (ug»0)  is  the  real  part  of  the  so- 
lution and  depends  only  on  if. 

Recalling  (3:37)  and  (3:38),  in  conjunction  with  (3;56a), 
it  follows  that 


whore  now 


(3:58) 
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Consequencly,  the  operators  011^(1  = 1)2,3)  are  non 
Bessel  operators  and  solutions  for  Cf  and  f are 

and  (3:59) 

= A3J„(40  Bj-fnCsrl.  (b) 

As  a consequence  of  the  completeness  theorem  given 
prevlousl;^,  equations  (3:55)  with  (3:59)  represent  Che 
general  solution  of  the  equations  of  motion  (2:15)  and 
(3:56b)  and  (3:57b)  are  the  resulting  displacement  fields. 

An  unusual  and  significant  fact  concerning  the  solution 
for  the  equations  of  motion  obtained  in  this  chapter  Is  that 
Che  rotational  and  irrocatlonel  portions  of  the  displacement 
field  both  rel7  on  and  ^2  definition.  In  the  isotrop- 
ic theory  (f2  does  not  appear  In  the  rotational  (equivolumi- 
nal)  field  and  does  not  appear  In  the  dilatatlonal  field. 
This  additional  complexity  of  the  anisotropic  theory  was  al- 
luded to  in  Section  1 of  tills  chapter. 

Writing  the  dilatatlonal  and  rotational  fields  as  the 
divergence  and  curl,  respectively,  of  Uj^  produces  the  fol- 
lowing expressions 

ui;, cosne  c.s(k* -wt)  , (3:60) 

“ij  (a) 

(3:61) 
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SmtiB 


(O  p:61) 


«?  = ' ‘t*.- 

Thus,  c(>j^  and  ^2  each  partially  dllatatlonal  and  equivo- 
Itiinlnal  in  character.  Since,  in  the  limit  to  Isotropy, 
passes  over  to  an  equivoluminal  potential,  it  uill  be  referred 
to  as  fluasl-eouivolunin.-xl.  SimiUrly,  <?2  >’®  referred  to 

as  miasl-dilatationnl.  Only  f(r)  retains  its  singular  rSle 
as  an  equivoluminal  potential  in  Che  transition  from  isotropy 
to  transverse  isotropy. 


CHAPTER  IV 


The  regular  Bessel  functions  and  Yj^Cpr)  are  os* 

dilatory  in  nature  (for  real  p)  and  the  number  of  zeros  of 
J in  some  Interval  of  r,  say  r,^  r<r,,  is  governed  by  the 
magnitude  of  p - in  a uay  similar  to  that  for  cos(pr).  Hence, 
the  nuignitude  of  p is  associated  with  the  number  of  waves  in 
the  interval.  The  numbers  p,  q,  and  s,  defined  in  equations 
(3:31)  and  (3:35),  will  be  referred  to  as  radial  wave  numbers 
as  opposed  to  the  axial  wave  number  k. 

The  general  solution  for  (f  and  f (given  in  equation 
(3:51))  will  be  directly  effected  by  changes  in  the  qualita- 
tive Tiature  of  p,  q,  and  s.  For  instance,  if  s is  imaginary 
f(r)  ceases  to  be  oscillatory  in  nature  and  behaves  in  an  ex- 
ponential fashion.  Similarly,  if  p and  q are  complex,  then 
(f(r)  exhibits  an  oscillatory,  exponentially  damped  behavior 
as  r is  varied.  In  this  eliapter  an  investigation  of  the 
range  of  values  to  be  associated  with  the  radial  wave  numbers 
will  be  conducted.  For  each  range  of  values  found  possible, 
an  appropriate  form  for  the  solution  will  be  specified. 

1.  Analysis  of  the  Radial  Wave  Humber  s.  and  f(r) 

Since  the  expression  for  s*'  is  the  sinqilest  of  the  three 
radial  wave  numbers,  it  will  be  considered  first. 
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(6:1) 


If  > 1,  then  9 > 0,  otherwise  is  zero  or 

negative.  The  magnitude  and  sign  of  s will  be  indicated  by 
absolute  value  signs  on  the  parenthesized  factor  in  (4:1) 
and  a sign  factor  such  that 

’ 1-1  <1. 


(4:2) 


Then 


In  this  way  the  value  of  6 is  always  positive  and  JTj  is 
either  +1  or  i for  the  conditions  listed  in  (4:2).  The  val- 
ue (s)  of  (uVk*  which  causes  a radial  wave  number  to  vanish 
is  of  considerable  interest.  This  number  (c^i,  this  in- 

stance) will  be  called  a turning  point  number  in  view  of  the 
change  in  behavior  of  tho  Bessel  functions  when  passes 

from  greater  than  to  less  than  the  turning  point  value. 

Using  (4:3)  and  (4:2)  in  Che  expression  for  f(r)  yields 

Based  on  the  relations  for  modified  Bessel  functions  of 
the  first  and  second  kind  C^^3i  f(:^)  assumes  the  forms 
f AjJ„(VsO  (a) 

((i)=|  (4:4) 

ifjej.i,  (b) 
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vhere  the  arguaent  of  s is  either  0 or'n/2.  Equations  (4:2), 
(4:3),  and  (4:4)  completely  specify  the  solution  for  f over 
all  permlssabla  values  of  real  u and  real  (positive  or  neg- 
ative) k. 

The  wave  equation  admits  solutions  for  which  the  fre- 
quency aj  may  be  real,  complex,  or  imaginary.  Physically, 
however,  It  Is  only  acceptable  for  Che  frequency  Co  be  real 
and  positive. 

Imaginary  or  complex  domains  for  Che  axial  wave  number 
k are  physically  accessable  and  acceptable.  If  the  axial 
wave  number  is  complex,  say  k»aelb,  the  imaginary  part  of 
k will  enter  the  wave  solutions  as  an  exponential  factor 
causing  the  amplitudes  of  the  waves  to  be  either  Increasing 
or  decreasing  in  magnitude,  depending  on  the  sign  of  b.  This 
attenuation  is  proportions!  to  the  axial  coordinate  z and, 
superficially,  implies  that  the  soureeless,  nondissipative 
medium  (in  whicli  the  wave  is  supposed)  is  either  adding  to 
or  absorbing  the  energy  transported  by  the  waves.  The  energy 
balance,  however,  is  not  affected,  as  Pilant  ^4]]  ^ibd  Oliver 
[67]  explain:  when  the  axial  wave  number  Is  complex,  the 

traveling  waves  of  given  frequency  and  wavelength  occur  in 
pairs,  traveling  In  opposed  directions.  This  phenomena 
creates  a standing  disturbance  with  zero  group  velocity. 
Consequently,  energy  transport,  which  is  proportional  to 
group  velocity,  is  lacking.  Standing  waves  with  complex 
axial  wave  number  occur  at  the  periphery  of  elastic  bodies 
and  have  been  observed  by  Oliver  [67],  Once  [68],  Zemanek 
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[69],  Shau  [70],  and  others. 

When  k is  a pure  Imaginary  number  the  propagating  com- 
ponent of  the  wave  solution,  that  Is,  exp  (ItRelkl- lut) , 
vanishes  and  is  replaced  by  exp  (lAit),  which  Is  characteris- 
tic of  spatially  attenuated  vibration  solutions. 

It  is  emphasized  that  when  k is  complex  the  bearing 
of  the  attendant  solutions  on  the  case  of  traveling  waves 
in  an  infinitely  long  cylinder  is  lost  since  waves  travel- 
ing in  a nondissipative  medium  must  vanish  if  the  amplitude 
contains  an  axial  attenuation  factor  of  the  form  exp  (-bz). 

An  oddity  in  the  anisotropic  theory,  lacking  in  isotro- 
py, is  the  occurrence  of  complex  radial  wave  numbers  for 
reel  values  of  u and  k.  In  this  instance,  also,  there  is  no 
energy  loss.  If  p and  q are  complex  for  real  k,  then  the 
wave  propagates  without  axial  attenuation.  The  attenuation 
is  confined  to  the  radial  direction  of  the  field. 

For  completeness  the  radial  wave  numbers  p,  q,  and  s, 
and  the  solutions  for  and  f will  also  be  presented  for  the 
cases  k complex  or  pure  imaginary. 

Let  k = a + lb.  Then 


If  bdO,  then  is  complex  and  if  keib,  s^  is  real  and 
positive  for  all  values  of  b.  Consequently,  the  representa- 
tion for  f(r)  is 


(4:5) 


(4:6) 
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(b) 


(4:7) 


The  occurrence  of  Bessel  functions  of  eoniplex  argument 
in  Che  solution  for  f(r)  may  be  vieued  in  the  same  light  as 
a sinusoidal  function  of  complex  argument  - for  the  analytl* 
cal  behavior  of  the  two  is  roughly  comparable.  In  Section  6 
of  Chapter  V the  effects  of  this  type  of  solution  will  be 
fully  disclosed.  For  the  moment  it  is  sufficient  to  indicate 
that  Che  contribution  to  cho  displacement  field  by  such  a so- 
lution will  be  of  an  oscillatory,  exponentially  decaying  na- 

To  continue  tlie  discussion  concerning  s and  f(r),  Che 
representation  for  f when  °°44>  bhaC  is,  the  turning 

point  value,  will  be  given.  Since  s‘  «=  0 for  this  value  of 
Che  Bessol  operator  reduces  to 


This  is  an  example  of  the  so  called  Euler,  or  equldlmen- 
sional,  differential  equation,  with  solutions 

f V"  ■,(  n*0.  (a) 

= j (q.g) 

I Aj  +Bjlor  if  n = o.  (b) 

It  is  instructive  to  obtain  this  representation  directly 
from  (4:4),  Recognising  that  f(r)  depends  parametrically  on 
the  value  of  k , f(r)  may  be  written  as 


0. 
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AjCWs"!  +Bjas1WV6r1, 

Uslns  the  amall  arsumenC  approjcimacions  [66j  for  J 
and  f(c)  becomes 

[ _Bj(wsU"Cn-iV.(v.%f"r~"  if  n+o, 

■"  (4:9) 

Comparing  (4:8)  and  (4:9),  the  follovlng  limits  must  hold 
M.«a  = A,(o) , n>o,  (a) 

Lim.t  - -a»  n>0,  (b)  (4:10) 

W-.®  lUtt’' 

Ui«H  f(*o.  (c) 

notice  Chat  (4;10c)  is  nonlinear  in  the  involvement  of  the 
radial  coordinate. 

Ihe  last  possibility  of  interest  occurs  if  k is  zero. 
From  (4:1)  it  follows  immediately  chat 

S(V*o)a  I.,:  » >1^  • (4:11) 

Since  Sj.  is  positive  without  equivocation,  €j  is  +1  and  equa- 
tion (4:4a)  is  the  appropriate  solution  for  f(r). 

Equations  (4:4),  (4:6),  and  (4:8)  comprise  a eoioplete 
listing  of  all  possible  solutions  for  f(r)  over  the  widest 
range  of  variation  of  r,  k,  and  t>, 

2.  Analvs is  of  the  Radial  Wave  Numbers  2 and  a,  and  (pCr) 

Attention  is  directed  now  to  equation  (3:31)  for  p^  and 
q^  . For  convenience,  this  equation  is  reuciccon  here 


4l'-(<#-'»)*M(ti-'"*)*(‘-*'«i‘r-  <•>  <‘-‘2) 


[a,b]  «iU  be  used  to  denote  a<  x<.  b,  a<  xs  b,  and  a4  x*  b, 


(4:13) 


is  in  [0,033).  4 c®"  "o*^  >’®  ®®>^o 

in  [C33,»)  since  all  three  terms  in  A are  strictly 


'’«L“  <a) 


and  (4:16) 

-.0  . (b) 

as  is  most  e.-»sily  seen  be  rewriting  - B as 


SI 


Purtherraore,  Hhen  is  in  Cc^ft>C33)>  - B again  can  not 

vanish.  This  is  seen  by  examining  the  form  of  A^  - B employed 
in  (4:lZa).  However,  when  is  in  [0,c^^),  A^  - B may 

pass  through  zero  (from  positive  values)  and  become  negative. 

The  qualitative  behavior  of  equation  (4:14)  is  described 
by  the  following  observations: 
a)  If 

both  (i  = 1,3)  are  complex. 


both  and  are  equal.  Their  magnitude  may  be 

obtained  by  substituting  (4:18a)  into  (4:14) 


Equation  (4:18b)  indicates  that  (^33*^44)^  must  be 
smaller  chan  (^13  + ^44)^  if  A*  -B  is  to  vanish  for  a positive 
value  of  ^w'/k^. 
c)  If 


r <,c  (^■V*<vy- 

" ** 

both  ^yk^lii  (it"  1,3)  are  real. 
(4:14)  is  less  than  the  term 


(4:19) 


If  the  radical  term  in 


(c,c,^-  (4:20) 

and  this  term  is  also  positive,  then  both  (1  = 1,3) 

ore  positive.  If  (4:20)  is  negative,  both  (i  = l,3) 
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are  negative.  Finally,  if  the  radical  terra  of  (d:14)  is 
greater  than  (4:20),  then  is  positive  and 

The  quantity  can  only  be  a positive,  real  number 

because  of  the  stipulation  that  k be  real  end  the  require- 
ment that  bj  is  positive.  Applying  this  criteria,  and  the  qual- 
itative behavior  of  (4:14)  outlined  above,  the  following  ob- 
servations concerning  the  qualitative  behavior  of  equation 
(4:12a)  hold: 

i)  If  Ci,4  + (<=i3+C44)'^/(=33-e44).  »<>  positive, 
real  value  of  will  cause  Che  radical  term  of  (4:12a)  to 

vanish,  and  consequently,  p’’  and  q’'  are  real  for  all  u and 
k in  E0.®)‘ 

ii)  If  Cii=C44+(Oj^3  + C44)*^/(c33-c^),  and  if 
is  equal  or  greater  chan  0,  the  radical  term  of  (4:12a)  will 
vanish  when  ^/k*  = j . For  this  ease  p*’  and  q^  are 

always  real  numbers  (positive,  or  negative)  for  all  sJ  and  k 
in  [0,»). 

ill)  If  C44  + (Ci3  + C44)*'/(e33  ”044)1  P°®* 

sible  for  the  radical  tern  of  (4:12a)  to  be  imaginary  for 
real,  positive  (u'/k^.  As  outlined  in  part  (c)  of  the 
previous  paragraph,  this  condition  prevails  whenever 
is  positive.  Moreover,  the  smallest  value  of  ^/k*'  for 
which  the  radical  term  of  (4:12a)  can  remain  Imaginary  is 
^'yk*'lj^^  or  soro,  depending  on  whether  is  positive  or 

negative,  respectively.  Consequently,  p‘‘  and  q^  ere  complex 


(conjugate)  numbers  when 


53 


The  above  discussion  has  been  illustrated  in  Figure  1 
b}*  the  curves  marked  and  These  curves  are  obtained 

b^r  inverting  (4:1&)  to  fom  the  following  inequality  in- 
volving Cj]_ 

tji  - 

Inequality  (4:21)  gives  Che  ceecrlcCions  on  c^^,  in  terms  of 
elastic  constants  and  ^/k'’ , such  chat  A*  - B in  (4:12)  will 
be  negative-  If,  in  (4:21),  ^uVk^'  is  varied  from  0 to  c^ 
curve  in  Figure  1 results  from  the  left  hand  side  of 
(4:21)  and  ^ from  Che  right  hand  side.  A horizontal  line 
drawn  through  any  value  of  Cj^j^  then  intersects  these  two 
curves  in  none,  one,  or  two  places.  These  intersections 
correspond  to  ^Vk’"!,^^  snd  ^<^/k^lj.i  , and  the  larger  (or  only) 
value  is  }“'/k'■^_|  . Notice  that  the  horizontal  line  also 
intersects  curve  which  is  equation  (4:25).  From  this 
it  is  to  be  noted  that  when  ?«Vk\^  is  in  [0,c^^)  it  lies 
intermediato  to  the  valties  given  by  (4:14).  As  previously 
noted,  wlien  is  greater  than  + ^^®33“®44^ 

there  are  no  intersections  with  curves  and  in 

consonance  with  the  requirement  that  k be  real. 

From  (4:21),  if  =0 

|A  t c„t  (4:22) 

and  if 

c = 


(4:23) 
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FIGURE  1. 


(1-1,2, 3) 


Che  radia 


(curve  (T))  and  Inequality 
) sbovlns  Che  llnlcetlons 
.1  wave  numbers  may  (the 
: be  complex. 
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Equation  (4:22)  sets  limits  on  such  that 
is  always  nesative.  Equation  (4:23)  sho»s  the  requirement 
on  Cj^j^  for  uhich  there  is  no  range  of  values  of  left- 

adjacent  to  such  that  p*-  and  are  not  complex. 

Before  returning  to  a description  of  the  solutions  for 
i^(r),  it  is  worth  while  to  discuss  the  ordering  of  the  two 

Because  the  elastic  constants  are  independent,  end  be- 
cause Cj^2  i&  not  limited  to  positive  values.  It  is  not  clear 
under  what  circumstances  one  or  the  other  ratio  la  the  larger. 
By  assuming  arbitrarily 


and  inquiring  into  the  truth  of  this  inequality,  the  follow- 
ing requirements  are  found  to  exist  between  Cj^3>  '^33>  ^/j/, 

(4:25) 

But,  the  positive  definiteness  of  the  strain  energy  density 
requires  c^>  0,  henco  the  quadratic  must  be  positive. 

With 

(a) 

and  (4:26) 

= -•^{4C,^-3C„’l  t •g\JCrt(BC,^t3C„Y  , (b) 

equation  (4:25)  yields  to  the  form 
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ThuG,  if  (4:24)  holds 


Hhen  the  two  branches  of  the  parabola  (4:25), 

given  by  and  join  and  (4:24)  holds  for  any  C33, 

if  Cj^3<  -9C33/8. 

Figure  2 illustrates  these  facts  about  (4:24).  The  co- 


i lines  have  been  taken  a 


=44/033  ana  C33/C33,  eno 

t33  by  assumption,  the  maxisum  value  of  the  ordi> 
nate  can  not  exceed  1.  The  requirement,  from  (4:18),  that 

be  positive  when  Cj^j^  is  given  by  (4:17)  is  indicated 
by  Che  curve  Isbled  For  material  constants  with  values 

to  the  right  of  this  curve,  is  positive. 

Turning  now  to  the  radial  wave  numbers  p^  and  q^,  it  is 
noted  (see  Fig.  1)  that  only  when  is  in  C=44>e33)  is 

it  possible  for  there  not  Co  be  a complex  range  for  p^  end 
. Transversely  Isotropic  materials  for  which  p®-  and  are 
always  real  (for  example,  magnesium)  will  be  called  less  an- 
isotropic in  analogy  with  isotropic  smcerials.  If  I® 


•44.C33. 


,)  end  thee 


are  complex 


n the  range  of  p^ 


, positive  , tbo  anisotropy  will  be 
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Within  the  shaded  region  the  reverse  of  (4:24)  holds. 
Curve  ® delimits  the  regions  of  material  constants 
such  that  to  the  right  of  the  curve  the  quantity  toZ/k'li,, 
of  equation  (4:18b)  Is  positive. 
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tormsd  more  anisotropic,  class  1.  If,  on  Che  other  hand, 

^ C0»C/|/|)i  there  are  always  complex  zones  for 
and  and  such  anisotropies  will  be  referred  to  as  more 
anisotropic,  class  II.  Class  1 more  anisotropic  materials 
exist  When  ever  the  material  constants  satisfy  the  follow- 
ing Inequality 


Class  II 
ever  the 


< c,  < c..,,  H. 

anisotropies  (always  more  aniso 
material  constants  satisfy 


(4:29) 

splc)  exist  when- 
(4:30) 


On  the  basis  of  this  description  any  material  isay  be 
classified  In  a manner  which  brings  to  mind  the  behavior  of 
the  radial  wave  numbers,  for  instance  beryllium  and  zinc  are, 
respectively,  eleae  I and  class  II  more  anisotropic  materials. 

The  classification  has  further  utility,  as  will  now  be 
discussed.  Consider  class  1 materials  end  recall  that  Che 
turning  point  numbers  are  chose  values  of  %‘^/k  for  which  a 
radial  wave  number  vanishes.  Referring  to  (4:12)  it  is  seen 
that  when  ^/k  P*  9* 


(4:31) 


(h) 


Thus  =Cgj  is  a turning  point  number  for  q. 

For  in  the  interval  »°33^  - B term  is 
elways  greater  than  A‘  (-B  is  positive),  hence  q‘  is  nega- 
tive throughout  the  interval,  and  p^  renmlna  positive. 
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Sines  materials  of  elass  1 are  under  consideration,  it  is 
apparent  that  uhen  ^u>/k  <=  c^  clio  A term  is  negative  and 
(a) 

(4:32) 

= 0-  (b) 

Thus  {w/k  =c^  is  a turning  point  for  p. 

As  decreases  Into  Co,c^)  both  p^  and  become 

negative  since  A Is  negative  and  now  A^  - B Is  less  than  A^. 

If  the  material  is  less  anisotropic,  that  is,  if 


p*'  and  q^  remain  negative  throughout  the  interval  C0,Oj^^); 
otherwise  p*"  and  q^  pass  from  negative  real  numbers  to  com- 
plex conjugate  numbers  when  and  back  to  neg- 

ative real  numbers  when 

For  f»Vk*  contained  in  C?“Vk*^lu  ,“) , P and  q are  real 
until  ^/k’'  decreases  below  their  respective  turning  point 
numbers,  at  which  point  they  become  imaginary.  Uhen  is 

in  (^yk'lvs  P <J  negative  complex  conjugate 

numbers , that  is , 

q,  = -P*.  (4:33) 

where  the  asterisk  means  conjugate.  The  range  of  the  argu- 
ment of  complex  q is  "n/2  < arg(q)  ‘ 3^4,  and  for  complex  p is 
tr/4  < arg(p)s 'n/2.  Figure  3 illustrates  the  behavior  of  p,  q, 
and  pi  , qi"  for  class  1 materials  when  0. 

Consider  class  II  materials.  As  for  class  I materials, 
{u>yk’’  “Cjj  is  a turning  point  number  for  q,  resulting  in 
equations  (4:31)  for  p*  and  q*  . When  ^•^k’'  is  in 
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for  class  I more  anisotropic  materials. 
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the  A term  of  (4:12)  reroains  positive  and  A*  - B greater  than 
A*.  Thus,  is  positive,  negative.  When  “®44> 

q again  vanishes  and 

(4:34) 

CL^o.  (b) 

As  ^u4'k‘  decreases  into  C0,e^^),  q^  becomes  positive  again 
since  now  f?  -B<A^  and  A>0.  Since  this  is  the  class  II 
ease,  there  is  always  a » 

both  and  q*  are  complex.  When  real 

part  of  p^  and  q^  vanishes  and,  as  decreases  below 

this  value,  reappears  negatively.  As  previously  discussed, 
the  complex  interval  extends  either  to  or  0.  Finally, 

when  ^u^^k*  5 ? 0,  p^  and  q*  again  become  negative  real 

numbers . 

For  ^/k^  contained  in  C^Vk^|^,“>),  p is  always  real 
end  q passes  from  real  to  imaginary  to  real  as  fsi^k^  passes 
through  the  turning  point  numbers  Cjj  and  c^.  When  ^u?^k*  is 
in  P «l  negative  complex  conju- 

gates, related  as  in  (4:33),  The  arguments  of  p and  q are 
bounded  such  that  0'arg(p)<1/2  and  H/I  * arg(q)  s H.  Figure  4 
shows  this  behavior  for  p,  q,  and  p^,  q^  for  class  II  ma- 
terials with  ?"'yk‘l^j^<  0.  Both  Figures  3 and  4 are  non-scale 
sketches,  included  strictly  for  illustrative  purposes. 

Employing  the  symbology  of  equation  (4:2),  with 
associated  with  p^  and  ^2  ^itb  q^,  and  the  extended  defi- 


nitions 
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FIGURE  4.  Radial  Wave  Kumbers  versus 

The  botcem  sketch  shows  the  ranges  of  p and  q relative  to 
each  other  and  relative  to  material  constants  and  ^u’/k’ 
for  class  II  more  anisotropic  maCerials. 
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it  .c„  , 

-.f  , 

if  fu'AX>rVk‘ 

il  fwVi.’^>  ^wVk' »o, 


i?  ^ ?-=Vk?l,,  , 

1?  fwVl>]^,?fuVk.*  >0. 

-,f 

if  > JcoVl?  > , 

;t 

If  »o, 


> c„, 

,-WL,>5.W^5-Vl‘L,. 


(a) 

(4:35) 

(b) 


(a) 


(4:36) 

(b) 


p “ , 


(a) 

(b) 


(4:37) 


II  ,.Vk'  1=  oom.i«.a  1«  [0.s.»«‘U]  » [f"Vl>n„.«)i 


(4:39) 


i,  (dn.*  1=  1.  (s«»i,‘L,,^f/k"L,). 


+ -V  AjZ„(VvUtkW„VV  (4:41) 

Here  Z^  and  are  regular  or  modified  Bessel  functions  of 
the  first  or  second  kind,  respectively,  depending  on  whether 
(i  = l,2)  is  +1  or  -1.  Notice  that  equation  (4:4)  may  be 


I ^ (a) 


IwCMB.lnryfc.j.c*,',^  AiI,Uy,r)  ■*  e.KoCW^O, 


f\zo;»=o, 


(4:42) 


( + 6Xa..’l/c„,cUl  - A.r'’  ^ B.r"'  , f^>o  • y,  .0, 
•PW-  (4:43) 

(k.(A,J.Ck*r')  + S.Yoa«r)ya.jK«)*A,-»5^V«r  ,n.o-,  ^=0. 


f equacions  (4:10)  v; 


SGrs  will  be  equal  when  jJ/lc"  = 

fi  (4:44)  holds,  the  dlffercntlcl  equation  for  degener> 

(4:45) 

a [71]  of  this  hyper-Bessel  equation  is 


is  only  applicable  it  0. 


(4:47) 
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’ l('  K,‘  (Ir.  • ')  * 1 

‘ <C  * T ‘“{ '<>'  • l‘iS>  (ir,  ’ '1  * - ‘«)1  • 


-.'-liMg. 

Equation  (4s41)  Is,  for  t 


1 (4:41).  (4:42),  (4:43),  and  ( 


R - + 4(«b0^‘  ± 2^  LCfl^  r>»A“.e  - iabDi.nsl  . 


D=  -c„Cii-ci  , 

If  neither  a nor  b is  zero,  p and  q are  complex.  If 
0,  (4:49)' reverts  to  (4:12),  and  if  a = 0,  (4:49)  gives 

('*  = 50) 


of  cp(r)  for  the  whole  range  of  k is  given,  simply,  as 
C^Crt.  k{A,H.„(pr),.B,W„lpriyfc«,cW)  ^ f'.^nCs.’-^-fB.WnCv),  (4:51) 


by  (4:49). 


CHAPTER  V 


1.  Preliminaries  of  the  Problem  and  Hondlmensionalization 

Consider  a right  circular  cylindrical  shell  of  infinite 
axial  length  (see  Figure  5).  Let  the  shell  consist  of  three 
layers,  each  transversely  isotropic  and  perfectly  bonded  to- 
gether at  their  mutual  interfaces.  Let  the  inner  material, 
bounded  by  the  radii  r = a,  r = b,  be  referred  to  as  material 
1,  and  all  equations,  elastic  constants,  scalar  potentials, 
etc.,  when  uricten  for  material  1 be  superscripted  as  ”(n". 
Similarly,  for  the  layers  bounded  by  r = b,  r = c,  and  r»>c, 
r=d,  let  the  superscripts  ''iii"  and  "tsl",  respectively,  be 

The  axial  uave  number  h and  frequency  u of  harmonic 
plane  waves  traveling  in  the  axial  direction  along  the  cylin- 
der will  be  the  same  for  all  layers  since  the  layers  are  per- 
fectly bonded.  The  eolution  of  the  anisotropic  equations  of 
motion,  prasented  in  Chapter  II,  is  valid  for  each  layer  of 
the  cylinder. 

Before  proceeding  with  the  development  of  the  frequency 
equation,  all  equations  to  be  used  will  first  be  nondimen- 
sionalizcd.  Executing  the  nondimenslonallzatlon  at  this 
stage  will  eliminate  the  necessity  of  developing  a subsequent 
notation  for  the  numerical  analysis. 
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FIGURE  5.  Cross-Section  o£  Cylinder 
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For  an  Infinite  flat  plate  of  thickness  h,  the  sheer 
velocity  of  the  first  antlsynmetrlc  thickness  shear  mode  Is 
fpf  , the  frequency  Is  F?j7/hjf  [72].  By  analogy  the  fre- 
quency given  by 


(5:1) 


will  be  used  to  nondljnanslonallzo  w.  The  tern  is 

the  velocity  of  equlvolumlnal  waves  traveling  In  the  direc- 
tion of  the  elastic  symmetry  axis  of  a transversely  Isotrop- 
ICi  infinite  body  [73]  with  the  mechanical  properties  of 
layer  2 of  the  cylinder.  Hence,  (5:1)  Is  the  first  mode 
thickness  shear  frequency  for  a plate  of  materiel  2 and 
thickness  of  the  middle  layer.  In  Che  limit,  ftirthermore , 
as  the  constants  pass  to  the  Isotropic  values,  (5:1)  reduces 
exactly  to  • Using  (5:1),  the  nondlmenslonal  fre- 

quency of  the  traveling  waves  Is  defined  as 


(5:2) 


defined  as 


(5:3) 


where  L Is  the  half  wavelength.  The  nondlmenslonal  variables 
(1  and  't  will  be  retained  throughout  the  analysis  as  the  de- 
pendent and  Independent  variables,  respectively. 

Denote  by  the  quantity 

‘ '>^'’1.  (5:6) 

and  by  the  ratio 


. (i  . 1,1, «. 


(5:5) 
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ote  by  V™  the  quantity 
co  ratio 


E=oA.  (5  = 8) 

n (5:2)  and  (5:3)  one  obtains 

is  seen  to  be  the  nondi.mens  Iona  1 phase  velocity  of 


'“4'  --i' 


a-i.i.,3).  (5:10) 


by 


.“.(Cpr . 


and  (4:39) 

(4:35)  and  (4:36),  are  given  now  as 
.I  ZVb‘’>‘>i. 

.t  Z™Vb‘-’A 


b Z.Vb‘’’^  > O , 

zVb‘«‘  >E‘^Vb‘«, 

zL‘;Vb‘'’‘bz.Vb‘'‘’^>,o, 


.1  (a) 
(5:13) 
.11  (b) 


;f  EVk!'"  >d‘^V4’, 
it  dV’/di’ 

it  E‘"Vb‘'’AEVt5‘’‘ 

it  Z.SVb‘'>SE-Vk.i‘’Ao, 

•,t  >a5:'V4''. 

;{  df/di"  >EVb"Ai, 

it  . 

z.^;Vb‘“  ^EVh^’^-X). 
a equations  (5:13)  and  (5:14),  the  r 


I 
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from  (6:14)  to  be 


. (5:13)  and  (5:14),  a‘” 


=li|g.-)4(¥>-vf(|''Tr-*|:(S.-oa-')i" 

d in  [0,Ob‘«]  or  [z‘f/b«S»):  or 


If  zVl^""  is 


(S:18) 

> (5:18)  to  the 


(5:19) 


(5:21) 


f = 2^,  . (S!22) 

t“ln-K'£„CV6“’n  + B‘”wT,C’!g‘'2\. 


uS  = -I.,|'J„l;|ir\C(d?.d?)!}=i?  .||"\sm<.e=«hS'!U),  (5:24) 

. C’  5 (A‘,V'“-  4"f.*  Ji“)  ]->-  ™<5>  -a'l. 
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In  (5:24)  the  formulation  op(r)  + c‘^3> +‘(’2“’ 

was  used  along  with 


2,  Derivation  of  the  Freauenev  Ea nation 


In  order  to  complete  the  solution  of  the  dieplacoment 
field  It  is  necessary  to  evaluate  the  integration  constants 
contained  in  the  displacement  potentials 
i^{?)  and  f(f).  The  boundary  and  interface  conditions  of  the 
cylinder  are  used  to  establish  relationships  between  the  in- 
tegration constants. 

For  a traction  free  cylinder  the  surface  stresses  must 
vanish  along  the  inner  and  outer  surfaces , that  is 
=0  at  f -a  , 


and 


(5:27) 


-C^P  *=0  at  ? *a,  tj  . r.e.i). 

The  interface,  or  continuity,  conditions  for  a cylinder 
with  perfectly  bonded  layers  are 


o„i  -x':;>‘xi] 

and 

In  these  equations  the 
der  layers  have  been  Indlce 
(5:22)  with  r replaced  by  i 
In  order  to  employ  the 


(5:28) 


cA  f -h  , ti- t 
al 

:ondimenslonal  radii  of  the  cylin- 
id  by  a bar,  and  ere  defined  in 


tions  the  stress  components 
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i (2:14)  forT‘"  (j=r,e,z),  the 


<-c^,^r[r.rv 


f ih4" j coi  ne  cosC^i  -«) , (a) 


+^L]'pr+ 


(5:29) 

(b) 


5La(?£-Ql)(c) 


(5:30) 


where  L = length  and  Fwforce,  Nondlrensionelizing  the  in- 
tegration constants  so  that  the  displacement  field  retains 
tho  dlmeceion  of  length  yields 


In  order  to  employ  equations  (5:23)  for  qf”  and  f'"  in 
equations  (5:24)  and  (5:29),  differentiation  formulas  for 
Zjj  and  Wjj  must  be  developed.  Consider  Initially,  and 
suppose  Che  argument  of  this  function  is  yf,  where  y may  be 
real,  imaginary,  or  complex.  Z may  represent  either  a reg- 
ular of  modified  Bessel  function  of  the  first  kind,  hence 

Since  tho  modified  Bessel  function  derives  from  the  reg- 
ular function  when  the  argument  is  imaginary,  it  follows  that 
the  following  differentiation  formula 

dSnCyf’)  = -tyZn„(y7)-<-Gj.^(-,r-)  (a)  (5:32) 

holds  for  either  understood  to  be 

the  modulus  if  the  argument  is  pure  imaginary,  and  if,  fur- 
ther, £ is  a sign  factor  of  the  type  previously  employed  for 
p'’,  s''’>  shd  s‘";  which  is  -1  only  when  the  argument  is  Imag- 
inary. In  a similar  fashion  it  follows  that 

^r,^7r)  c -y  W„,|C7£)  * 5 Wn(yf),  (b)  (5:32) 

since  the  differentiation  formulas  for  V„(t^) 


(5:23) 


(5:24) 


making  use  of  (5:31),  and  using  (5:32) 


♦ A?uS’  •Jr)l-4“?A'“' 

1 of  (4‘’»Jr)[-tA“’FW„„(Sf!«(l . nw„ci;j»o]t 

0?n«/.(l;£“'plj  c«  „s  =.5(l;i-tHl,  <■> 


- Aft-'f  t s“Pi..,A  1"')  t aH„15s".1]  - 


' '‘'"d|v4“  *"*)  * 

+ *?TF  . 
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* A?(4“tj,“)  [aj“  .?  , {e4»„,.-,,  - 

-fSU'V.V«.4f. 

-v(2a“\Cr,-.V 


(a 


*sm  n0cos(?7-ak), 
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- Aj’n>jf£nCV\"r)  -S.j’o^rW^^’ff'f-iUosnes.nC^l.fliV  (f)  (5:33) 

■'  (cent) 

Equations  (5:33)  are  applied  in  the  boundary  and  inter- 
face condition  equations  (5:27)  and  (5:28)  to  form  a system 
of  18  honogeneous  equations  in  the  18  integration  constants 
Aj\  Bj  (1,J  » 1,2,3).  In  order  that  nontrivial  solutions  for 
the  integration  constants  be  obtained,  it  is  necessary  (but 
not  always  sufficient)  that  the  coefficient  determinant  of 
the  Integration  constants  vanish,  resulting  in  the  equation 
dellCijUo,  (l,j«l,l,...,is).  (5:34) 

This  equation  is  the  so  called  frequency  or  dispersion  equa- 

tion. 

As  mentioned  above,  the  vanishing  of  (5:34)  is  not  al- 
ways sufficient  to  assure  nontrivial  solutions  of  the  inte- 
gration constants.  Equations  (5:33)  are  not  always  the  ap- 
propriate solutions  for  the  stress  and  displacement  equations 
as  was  revealed  in  Chapter  IV.  Thus,  when  attains  the 

value  of  a turning  point  number  for  p*  , q^"',  or  a”,  or  when 
gt/jjCoi  such  a value  that  one  of  the  three  pairs  of  p'*' 

and  g'”  are  equal,  equations  (5:33)  fail  to  represent  the  dis- 
placement and  stress  field  components.  To  this  extent,  equa- 
tion (5:34)  is  insufficient.  When  viewed  in  the  broader  light 
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of  Che  seneral  and  particular  soluclona,  (5:34)  is  the  only 
condition  ioposabla  to  obtain  nontrivial  solutions. 

In  ensuing  discxissions  It  will  be  necessary  to  under- 
stand the  construction  of  the  deterainant  (5:34)  In  sone  de- 
tail. The  rows  of  (5:34}  will  be  written  from  the  boundary 
and  interface  conditions  arranged  in  the  following  order; 


row  condition 
1 ‘ 0 

2 ‘ 0 F = 0. 


The  colujims  of  (5:34)  will  be  written  for  the  18  integra- 


tion constants  arranged  in  the  (logical)  sequence 


CO  lumn : 


123456789  10  11  12 

i“  »;■'  aJ  b';'  aJ  b“  a‘“  b"  4“  b“  aJ  b" 
nt):  13  14  15  16  17  18  (5:36) 

a‘“  A™  A»3’  Bf 


, Cy  AB.  1 


c„  ^ A v[Adl"nCa-.i  - 


-^1?# 

- .jf’ 

-jr^Ao-WAB-Ap”'  A bV'  f. .-arilwAlV""!, 
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n4'\dr*dr)[';r  a 

, = -4'’LnCn-0 

- d'"Ur>-0  - 6^'’\’-V'’'QVi]w„C^V'’«), 

p:37> 
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C,s  = r,W„l^6«^ 

Ca.'-'"?2n(V«‘'’6-), 

Cg,=  -(d^S4‘')nZ.„Cf/S“0, 

Ce,=  -Cd's’  + 4r)nV„C>r/i<''0, 

C«s"  - nE„Ct;lt‘'’0, 


c,,  --  ? s [j;'’4V'"  - * 4'’i*»i?/i"'s, 


®7,i»  '=8,1>  ®9.i  ^ 


"7.1.  "8.1.  "9.1  ' 
(1  = 1.2 


''17 

(i  = 1,2,  . 


i*  °1Z  i = 13,14,  . . . ,18)  =Cg  J^,  Cg  J^,  Cj  ^ 

. ,6),  respectively,  with  5 and  <>i  replaced  by 

i.  <=15.1  = 12>J’-Ci,i.  -Cj.i.  -C3.I 

« ,6),  respectively,  ulth  a and  replaced  by 


Ci5  ^ {1  = 13,14 18)  = Cj^ 

,6},  respectively,  ulth  a and  ( 


Cjg  ^ (1  = 13,14 18)=C3  Cj  j^,  Cj^  j 

,6),  respectively,  vlth  1 and  <0  replaced  by 


3,  The  fl-C  Plane 


Hie  plane  of  coordinates  defined  by  the  nondlmenslonal 
frequency  fl  as  ordinate  and  the  nondlmensional  wave  number  ^ 
as  abscissa  will  be  referred  to  interchangeably  as  the  n-t 
plane  or  the  frequency  spectrum.  Equation  (5:34),  being  a 
transcendental  relationship  between  frequency,  wave  number, 
and  number  of  circumferential  waves  n,  yields  an  infinite 
family  of  values  of  A for  each  chosen  value  of  ? and  n.  If 
the  values  chosen  for  ^ are  dense  enough,  there  results, 
rather  than  discreet  A-t  point  sets,  an  infinite  family  of 
curves  for  each  clioice  of  n.  Each  of  these  curves  is  a branch 
of  the  frequency  equation.  Each  branch  of  the  frequency  equa- 
tion originates  at  a specific  value  of  n (called  the  cut- 
off frequency)  when  f is  zero  and,  as  % increases,  passes 
across  quadrant  1 of  the  A-t  plane  (see  Figure  6).  The 
branches  are  frequently  given  names  which  are  indicative  of  the 
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character  of  the  dlaplacement  field  at  the  cut-off  fre- 
quency of  the  branch  (for  example,  flexure  mode,  breathing 
node,  etc.).  Although  the  displacement  field  will  vaiy  mark- 
edly as  is  increased  along  a particular  branch,  the  cut- 
off frequency  names  are  retained  for  purposes  of  identifi- 
cation and  convenience. 

The  coordinates n and  ^ have  been  chosen  to  depict  the 
behavior  of  the  frequency  equation  because  of  the  volume  of 
information  revealed  by  the  branch  curves  when  plotted  against 
this  choice  of  variables.  Kot  only  are  the  values  of  the 
frequencies  for  each  value  of  X apparent  at  a glance,  but 


it  is  apparent  that  the  phase  velocity  of  a traveling  wave 
of  given  wavelength  and  frequency  is  available  from  the  spec- 
trum as  the  slope  of  a line  from  the  origin  through  the  point 
with  coordinates  A and  S • Purthormore,  since  the  group  velo- 
city (the  velocity  with  which  energy  is  transported)  is  ob- 
tained as  the  derivative  of  the  frequency  with  respect  to  the 
wave  number,  it  follows  that 
dtO  _ ,,  dfl 


(5:39) 


Hence  the  slope  at  any  point  on  a branch  of  the  frequency 
equation  is  the  nondimcnslonalized  group  velocity  of  a wave 
with  the  frequency  and  wave  number  given  by  these  coordinates. 
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}1lBeorlcall7  the  branches  of  Che  frequent  equation  uere 
plotted  in  coordinates  of  phase  velocity  versus  wave  nuaber. 
Since,  at  the  cut-off  frequencies,  the  phase  velocity  is  in- 
finite, such  plots  could  not  reveal  the  value  of  the  cut-off 
frequencies • The  frequency  spectrum  circumvents  this  problem. 

From  the  point  of  view  of  analytical  (as  opposed  to  numer- 
ical) analysis,  there  are  ocher  advanCases  to  the  use  of  the 
fl-^  coordinates  for  analysis  of  the  frequency  equation.  Chap- 
ter IV  revealed  chat  the  turning  point  numbers,  the  commence- 
ment, and  the  termination  of  complex  domains  of  the  redial 
wave  numbers  (for  real*^)  are  all  readily  identified  in  terms 
of  values  of  the  ratio  fl/f.  Thus,  on  the  frequency  spectrum 
the  turning  point  numbers,  and  complex  zone  numbers  represent 
rays  of  specific  slope  extending  from  the  origin  of  coordinates 
across  quadrant  I of  the  n-\  plane.  There  ere  three  turning 
point  linos  and  two  complex  zone  lines  (for  more  anisotropic 
materials)  for  each  layer  of  the  composite  cylinder.  Any 
analysis  must  necessarily  be  concerned  with  the  ordering  of 
the  slopes  of  this  array  of  lines  since  the  nature  of  the  ele- 
ments comprising  the  frequency  equation  determinant  vary  as 
n/V  ranges  in  its  values  from  greater  chan  Co  less  than  their 

Since  the  different  degrees  of  transverse  isotropy  under 
consideration  may  or  may  not  produce  complex  radial  wave  num- 
bers, but  always  admit  the  Cvuming  point  values,  Che  ordering 
of  the  turning  point  lines  will  be  considered  initially.  The 
turning  point  line  for  s'"  has  the  same  slope  as  that  of 
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(of  the  same  layer)  or  the  lowest  of  q"’  if  g''’  does  not  van- 
ish, consequently,  there  are  at  best  six  distinct  rays  eina- 
nating  from  the  origin  of  the  frequency  spectrum  which  will 
be  turning  point  lines.  Those  rays  divide  the  frequency 
spectrun  into  seven  sectors  (Figure  6),  the  lowest  being 
numbered  1 and  the  others  numbered  counter-clockwise.  The 
slopes  of  turning  point  lines,  which  form  the  sector  bounda- 
ries, are  the  six  numbers  b“’, ’b‘"  (ic'1,2,3),  Moreover,  for 
a given  layer  ’i"  > b‘”  since  Cj j > c^.  Since  this  is  the  only 
restriction  in  the  ordering  of  the  turning  point  lines,  there 
are  90  distinct  ways  that  the  seven  sectors  can  be  bounded. 
Many  of  these  90  orderings  have  little  other  than  pedagogical 
significance. 

The  field  of  90  orderings  of  the  sector  lines  is  nar- 
rowed by  the  following  observations.  Light  weight  sandwich 
constructions  utilize  low  strength,  low  density  cores  which 
carry  transverse  shear  in  their  primary  role  of  separation 
and  stabilization  [74^  of  the  outer  layers.  The  outer  layers, 
or  facings,  on  the  other  hand,  are  high  strength,  more  dense 
materials,  Chu  [42,753  and  Yu  [41,76,773  have  used  such 
characteristics  for  the  three  Layers  in  their  approximate 
theories  for  light-weight  anisotropic  and  isotropic  sandwich 
shells.  Since  c^^  and  c^'^  ere  directly  related  to  the  axial 
Young's  modulus  and  shear  modulus,  respectively,  and  ace  also 
proportional  Co  Che  velocity  ratios  and  b'‘\  light  weight 
sandwich  shells  will  be  specified  to  be  chose  shells  for 
which 


a=i.M. 


:idcd 


(5:40) 

the  following  six 


This  requironenr  can  be  expan 
orderings  of  the  turning  point  line  slopos,  none  of  which  re- 
strict tho  theory  to  the  requirement  for  identical  facings 


fe”'  < < t,'’'  < 

t'*-’ < ■b"-'' < < b"’  < Ts'"  (a) 

and  ' light  freight  (5:41) 

Bnnrtwirh  Rhnilfi 

r b"'  < v’<  b"'<  V" 

I b‘'’  < b'”<  t‘'*<  ti”’  (b) 

[ < h^')  < 0="’  < -b"'. 

If,  on  the  other  hand,  the  core  layer  of  the  cylinder  is 
a high  strength,  high  density  material  relative  to  the  facing 
layers,  then  the  axial  Young's  modulus  and  shear  modulus  of 
the  core  will  be  larger  than  the  corresponding  moduli  of  the 
layers.  Such  core  layered  shells  will  be  restricted  by  the 
following  ordering  of  the  velocity  ratios 


The  above  two  listings  of  tiuming  point  line  slopes  do 
not  include  the  important  class  of  shells  for  which  the 


igth  of  the  individual  layers 


sparable.  Ag< 
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possibilities  ulll  be  listed  - based  on  the  restriction  that 
the  facing  layers  be  of  the  setae  material.  Thus 

b"’  c b‘”  , = 1^’', 

b"-'  < < b”’  < T,™. 

fa‘"  < b“*  < t'"  < 1>‘",  <a) 

b‘”<  b‘"  < t'"  < 

and  aenernl  shells  (5:43) 

b'’’  < < b‘"  < , 

k'«  < b"'  <'b"*  < -b'",  (b) 

b'’-''  a b'”  <t‘"  < V'-’. 

Inequalities  (5:43)  actually  bridge  the  gap  between  (5:41) 
and  (5:42),  that  Is,  the  first  of  (5:43a)  corresponds  to  all 
of  (5:41)  and  the  first  of  (5:43b)  corresponds  to  all  of 
(5:42)  for  the  case  of  Identical  facing  materials.  Similarly, 
If  all  three  materials  are  the  same,  the  last  of  (5:43a  and  b) 
correspond  to  the  single  requirement  for  the  turning  point 
line  slopes  of  a single  layered  cylinder  If  the  first  and 
third  Inequality  signs  are  changed  to  equality  signs.  All 
eighteen  orderings  listed  above  in  no  manner  affect  the  ge- 
ometry of  the  shell  layers,  but  do  represent  a cross  section 
of  useful  possibilities  of  layer  strength  end  density  ratios 
which,  if  properly  pursued,  will  reveal  the  trend  behavior  of 
the  90  possible  cases. 

Tiiming  now  to  the  slopes,  z"  and  z\_,,  it  was  estab- 
lished in  Chapter  IV  that  (a)  only  for  more  anisotropic  ma- 
terials will  these  lines  exist,  and  (b)  the  slopes  of  these 
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FIGURE  6.  A-\  Plane  with  Zone  Linea 

The  iI-<;  plane  for  a light  weight  sandwich  shell  depicted 
by  the  velocity  ratio  ordering 

< "U"  t < T5"  ‘ Tf’> 

of  equation  (5:41b).  The  outer  layer  is  assumed  to  be 
class  1 more  anisotropic,  with  b''*<  Zt?<  b‘®  and  z'J,  c b'*’. 
The  inner  layer  is  assumed  to  be  class  II  more  anisotro- 
pic, with  b"‘ and  0.  The  middle  layer  is  leas 

anisotropic.  The  broV.en  lines  are  complex  zone  lines. 
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lines  ere  elvays  less  than  the  smallest  turning  point  line 
slope  for  a given  material.  Restrictions,  relative  to 
and  the  three  constants  and  «ere  listed  in 

Chapter  IV  by  which  it  can  be  ascertained  whether  the  lowest 
complex  zone  line  has  zero  or  positive  slope.  Consideration 
of  these  facts  in  light  of  inequalities  (5:41),  (5:42),  and 
(5:43)  indicates  that  a complex  zone  for  an  outer  layer  of 
a light  weight  sandwich  shell  may  have  its  upper  boundary 
line  in  any  sector  from  1 thru  5,  tho  same  being  true  for  the 
core  layer  of  the  shell  class  described  in  (5:42).  Uithin 
tho  frameworlt  of  (5:41),  (5:42),  and  (5:43)  there  is  no  more 
specific  way  to  indicate  the  occurrence  of  the  complex  zones 

The  above  comments  indicate  the  need  for  the  proper  re- 
presentation of  the  when  tl/f  has  a turning  point  line,  or 
complex  zone  line  value.  In  the  next  two  sections  of  this 
chapter  these  formulations  will  be  obtained. 

4,  Turning  Point  Line  Expressions  for  the  Cj^j 

As  mentioned  above,  there  are,  in  general,  6 distinct 
values  of  C2/f  which  are  turning  point  numbers.  Only  the  wave 
number  vanishes  for  the  turning  point  number  which  is 
the  largest  of  the  two  for  a specific  layer,  and  for  class  1 
or  II  anisotropies,  respectively,  s*’’  and  g*"  or  s”  and  g]" 
vanish  for  the  lower  value,  b‘'\  Consequently,  recalling  the 
construction  of  the  frequency  equation  deteminent,  it  follows 
that  the  are  effected  in  groups  of  two  or  four  columns 
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when  Z attains  a turning  point  value.  Thus,  any  given  turn- 
ing point  number  only  affects  a minor  portion  of  the  fonnu- 
latlon  of  equation  (5:37). 

The  proper  form  of  the  can  be  obtained  by  use  of  the 
appropriate  solution  for  or  f*'*  given  in  Chapter  IV,  for  ex- 
ample, if  g”'  = 0 equation  (4:8)  can  be  used  in  the  appropriate 
portions  of  equations  (5:33)  to  rederive  the  for  columns 
17  and  18  in  (5:34).  This  method  is  tedious  since  the 
found  by  the  above  process  will  be  unpredietably  different 
depending  on  whether  n is,  or  is  not,  zero. 

The  awkwardness  of  this  situation  is  eliminated  by  utili- 
zation of  equations  (4:10),  which  indicate  the  alteration  in 
the  form  of  the  integration  constants  as  a radial  wave  number 
vanishes.  These  equations  can  be  used  in  applying  the  limit 
to  the  of  equations  (5:37)  to  obtain  the  proper  turning 
point  line  formulations.  Reference  to  (4:10a  and  b),  for 
n(6  0,  shows  tliat  no  order  of  differentiation  effects  the 
form  of  the  constants  A(0)  or  B(0).  However,  when  n = 0,  the 
fact  that  the  limit  of  a derivative  is  not  necessarily  the 
same  as  the  derivative  of  a limit  must  be  entertained.  Equa- 
tions (5:33)  show  that  either  ef'‘'or  f”’are  differentiated  at 
best  only  once  in  obtaining  the  Hence,  for  n<°0,  the 

Bessel  functions  of  order  1 in  equations  (5:37)  derive  from 
zero  order  Bessel  functions.  Thus,  taking  the  limit,  as  the 
radial  wave  number  vanishes,  of  a Wj^  Bessel  function  must  be 
treated  as  taking  the  limit  of  the  derivative  of  (4:10e) 


rather  than  that  of  (4:10b) 


The  difference 
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6,(kOi -*  (5:44) 

Equation  (5:44)  ia  the  proper  limit  for  the  Bessel 
cions  in  (5:35)  for  the  axisyraetrie  ease  n = 0. 

In  order  to  utilize  (5:44)  and  (4:10)  in  (5:35)  the  ap- 


a 2/t  as  in  (5:44),  are  factored  out  of 
ance  with  the  limits  of  (4:10)  and  (5:44) 


IS.V 

c„- 


■ 'id  ' 


Iidl'ei;'  .Ji'Vdd.o  - di"  J?  ", 

^ uf.jn  - Ji“  4"  V'  5‘(-|l,-  0 . “ ■ ” . 


(5:45) 
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C„-  -nvadrs", 


:45J 

(cont] 
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Ctj= 


Cft,= 

^-■lo,  .... 

Cgj  • -{d" +ct4'’)  rvfe", 
f-(d^'+di")n6"‘ 


Jn%  ^ 

“‘‘I-.,.-, 


C„  = C,.  =0  , 
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The  reniainlne  are  formed  from  the  ebove  by  use  of 
the  Instructions  Included  In  equations  (5:37). 

Equations  (5:45)  are  not  to  be  used  as  a Ilstlns  of  the 
general  elements  of  the  frequency  equation  as  are  equations 
(5:37),  Only  columns  from  the  frequent  equation  for  vhleh 
the  radial  vavo  number  is  zero  at  a specific  choice  of  Z 
ehould  be  composed  from  the  corresponding  column  elements 
listed  in  equation  (5:45), 

5,  Complex  Zone  Line  Expressions  for  the 

The  complex  zones  of  the  radial  wave  numbers  occur  only 
in  p"’ and  g“’,  hence  columns  (61  - 1)  and  (6i)  (i  = 1,2,3) 
of  the  frequency  equation  determinant  do  not  Involve  Bessel 
functions  of  complex  argument.  Equation  (4:45)  suffices  as  the 
solution  for  <4‘  when  2"  "9  • This  solution  may  be  written  as 

where  cp'”  is  the  same  as  used  previously  in  the  general  solu- 
tion but 

(5:46) 

Consequently,  only  columns  (61  - 3)  and  (61  - 2)  (ic  1,2,3) 
are  effected  by  the  occurrence  of  complex  zone  line  values 
for  Z,  The  now  forms  for  the  altered  C^j  will  be  listed  be- 
low, It  is  important  to  recall,  in  using  these  terms,  that 
since  the  three  layers  are  generally  dissimilar,  only  two 
columns  of  the  frequency  equation  determinant  will  be  ef- 
fected for  any  given  complex  zone  line  value 


of  Z. 


C,,^  24"t4'dOVs’''5^W„^;t^E",5^  2 Q^24\4%OC8*CCn-Q- 

- d“t' a>U*el" . 4' Jl- 51, 

c„  . 2 j;>  a 5(4"  ,“  a ^ 

c„-  24“an(4'’.<;‘)[l;«“5«,„(t.”5)-(ti,0W„„l55Mal], 

c„.  a;"a'sld“.“5“-d:'<.-. 

c,»  • ldS\j."H[-1s"tW.„Ci«l  Hd>0«„Ra"ll] , 
C„.-nUd“-.jr)£.„II«“5,l, 


C„.  dS(J».d«)W„„l,5"a, 

c„-  ft  Idl'd?’ »"•  -d;<tf^+d;']z„„(,4«i), 

c„.  fk[j“t»d“ 


8,  9.  15. 


by  R 


e =X;(R«*UO, 

^--tf(R»-lU\ 

d„<pf')  = Rj„c?-fwaj„('f-o , 

Y„(prt  = R , 

j.i^h-hi"j:(i!H, 

v.cutM-iin.’iff), 


3 complex  conjugate. 
Erooi  (5:37)  are  noted 


(5:48) 


<1, -a»  = d,R  e«  * U,l«»  = Rt,  *iir. , 

■),-/- = EFi->UF,, 

(5:49) 

RF,-llFi, 


equation 


in  columns  (6i  - 5), 
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(6i  - 4),  (6i  - 3),  and  (6i  - 2)  (1  = 1,2, 3)  when  and 
are  complex.  Study  of  Che  element  llering  in  (5:37)  showa 
that,  element  by  elereac,  column  (61  > n)  and  (61  • n + 2) 

(1  = 1, 2, 3;  n = 4,5)  are  complex  conjugate  pairs  with  every 
element  in  the  g‘'  colusina  modified  by  the  aame  factor 
(d^  + d^’}/^  relative  to  Che  columns.  Assuming  that  1 = 1, 
the  frequency  equation  may  be  written  as 

del  ICgl  =dellCi,C.j^trtCv,Bl=  0, 
where  S is  Che  last  14  columns  of  the  frequency  equation,  and 
C^j  (J  = l,2,3,4)  ore  the  four  columns  which  are  complex, 
the  can  be  written 


Substituting  these  expressions  into  dctlCj^jl  and  expand- 
ing according  to  the  addition  rules  for  deceminancs  yields 
ditlC^jl  = -4deilCi,  Ci^CijCuBl  = 0.  (5:51) 

Thus,  when  a given  sot  of  g‘"  and  are  complex,  Che  4 
columns  affected  by  the  transition  can  be  replaced  according 
to  (5:50)  and  (5:51)  by  the  real  and  Imaginary  parts  of  the 
column  elements  of  which  p‘’’  is  the  radial  wave  number.  The 
same  substitution  could  as  well  be  made  in  terms  of  the 
column  elements.  Using  (5:50)  and  (5:37)  the  elements  of 
the  frequency  equation  which  are  affected  by  complex  g'*’  and 
g"'  are  given  as 


1* 

in  ff  r I n <S:S2) 


tidW^n|i;a[R»,"UY„. 


' I F "Ua  ( R»s'n  J„„  * I «">  RJ^  J - n IJ  J ]|  ^ 


nIJ„l  J, 

- IF?’ [-? a (R»“iY.„  • I«"  R V.)  • a I Y J 1 ♦ (cont) 
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-«-  R F,"  [-■;  5 Cr^"”  IV„..  -V  X^S"’  R'Xn,,)  ^ n I Y„  1 1} . 

c„  - t{  -ft  f r.RJn  li  * 

+IY{-Y^  * W"  Rj„,  1 - n 1 , 

?|-’?fc[R<‘‘''RV  * !>«'/„  I*  g^5Z) 


Cs.=  -«^'?RJ„-U-?^J„, 

- '■  ‘ Lkp.'"  1 \ ■'  in'"  «'<-.]  - 

The  argument  of  all  Bessel  functions  in  (5;S2)  Is 
in  ro»s  1,  2,  3,  and  in  rows  4 through  9,  The  comple> 


CipCij  (l  = i.2....,is-,  j .i,i,7,8,re.,*).  (5:5 

E,}uations  (5:45),  (5:47),  and  (5:52)  taken  in  conjunc- 
tion with  (5:37)  provide  the  complete  specification  of  the 
elements  of  the  frequeno’  equation  (5:34).  Additionally, 
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frequency  equation  by  the  limit  given  la  (3:53).  Complying 
with  this  limit,  all  terms  of  the  elements  of  (5:34)  con- 
taining (dj’+  d^')  in  either  numerator  or  denominator  vanish. 


CHAPTER 


The  frequency  equation  obtained  in  Chapter  V may  be  aL- 
teted  quite  siinply  in  several  specific  ways  to  obtain  the 
frequency  equations  for  other  useful  cases.  In  this  chapter 
Che  general  frequency  equation  for  Che  three  layered,  tractior 
free  cylinder  will  be  simplified  for  the  cases  of  axisymmetcy 
(n>=0).  Infinite  axial  wave  length  (^  = 0),  and  the  combina- 
tion of  both  n = 0 and  ^ = 0,  The  various  steps  necessary  to 
obtain  the  frequency  equations  for  cylinders  with  clamped 
(rigidly  held)  surfaces,  and  cylinders  with  combined  condi- 
tions of  free  and  clamped  surfaces  will  be  derived.  The 
number  of  layers  in  Che  cylinder  will  be  generalized,  and  the 
conditions  to  be  met  in  case  the  inner  radius  passes  to  zero 
will  be  indicated.  Furthermore,  the  condition  of  thin  outer 
layers  will  be  considered. 

1.  Freauenev  Equation  for  Three  Layered.  Traction  Free. 


Although  equation  (5:34)  passes  over  to  the  isotropic 
case  automatically  when  dy»  d^’,  dj  d^'i  and  d^-»  d^',  the  ele- 
ments of  the  determinant  will  be  listed  below  due  to  the  fun- 
damental importance  of  the  Isotropic  elastic  theory.  This 
listing  is  accomplished  prior  to  tasks  set  by  the  opening 
paragraph  of  this  chapter  so  that  both  the  anisotropic  and 


104 


iotropic  equations 


dV’.;"— . . u“. 

df-V-drS.-'? 

J(d;“,df'i«;"«“«. 

d».;v»‘.d“'-^>i^=.iCd?’.dV’). 

UslJis  equations  (6:1)  and  (6:2)  in  (5:37)  results  in  the 
following  formulation  for  the  elements  of  the  frequency  equa- 


(6:3) 

C.i  = Zdi”  fl  V a W„..c?  -"a)  V (n(n-v)  - cV'  , 


c»  - iiait/"5W„fTs»ai  -■5>a'UV’e"*-oV»iwa], 

C„ . Znd“[-Ii;’ir«“ai„, ftp'll , tn-.«.n»“a], 

C,t  - ^ad^^'ra'>QVn„C■^d.‘"a^  d Cr>-OWnC?aaQl), 

■ Cddrl.»t*“>a!S„,Cl*“»l  -<"-1E„«s.“5l], 

a..  ■ df’L'r*'"o  Si  «f"  »l], 

C„-  (S.J. 

d.i  • I„(W'«1  * -e;\>»«Vli„«.“all, 

Cdi-  . l(a«a-0  ^a4‘'’»a*i«.(Te"all, 

c..»  drs««?V«-ill-t»”'aW„Ci;t»ol 

c„-  -Id£’0't[-y,»“aw„f';e"’5l  • aw„i'B5«all, 
c„-  -aVddi'’2„(ta»al, 
c^-  -Btndi'VnUo'aaX 
c„  ■ il-aV’t»“bE„„Cts“a  » az.n»-»)l, 
a..  ■ «„,<1  «""‘l  E“  all , 
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Ct,  • ^ nW„C'?/i‘''fc’)], 

C-,5  = 

C„  ' 

Cg^‘ 

C6^  ' 

Cs,=  (6:3)^ 

Cgs=  feV’';4‘'’b2.n*,C';s‘''V) 

Cb(.'  •• 

C„*  65'’'?^®“'’^  BH-nCfe‘"f), 

C,*= 

C.,s»  C^=  0. 

Hie  remainder  of  the  olcments  of  the  frcquon^'  equation 
sre  formed  from  equations  (6:3)  end  tlie  Instructions  of  equa- 
tion (5:37). 

The  turning  point  line  expressiona  of  Chapter  V are 


directly  applicable  to  the  isotropic  equations  v 
por  portions  of  (6:2)  are  cnployedi  The  scope  < 
plicabiliryi  hovrever,  is  limited  Co  that 


3 the  values  of  n for  which  (5:34)  c 


n t = 0,  that  Is,  equations  (4:11)  a 

a;”a"Sj,»e 

-rf’t')-  0, 


(6:5) 


Finally,  since  p“,  q%  end  sf  are  in  sector  7.  only  reg- 
ular  Bessel  functions  will  appear  in  the  frequency  equation. 
Using  (6:4)  and  (6:5)  in  (5:37)  results  in 

ti,  » = 0, 

c,  - , [zi':'  i„(q|fa), 

c,,.  i<''l-<’aJ„,(s"a)  *cn-ov<”5)], 

C.S-- s"a  V„.,ti,''’aV(«-i)VnCs?  oil, 

^ii  * “ °i 

i„.nd;n,VcJ.„Cce-al-<n-.U„(<'ul, 

c„  ■ -d"  S?a  j.„Ci;’al  -fi'tnic-0  - ,S''taVi] 

C.C  ■ -J?  i'cYcttS'cl  -iS'l«n-i7  -i“6Vil-(.lsV'«l. 

Cl,  ■ ’"■'.'Ffa'l, 

Cii' 

=C,c, 
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^•7,  ' 0, 

+>n  , 

C,*-- 

n J„Cs4"6), 

c^*  r\V„Cs5‘’B'), 

Cej=  -nJn(<;'’BX  (6:6) 

(cont) 

C^- 

Cg^  = fe  ^ 77 fe) , 

Cgfe- 

c».  • 

C.i‘ 

The  remaining  elementa  of  tho  determinant  are  formed 
from  the  above  terms  according  to  the  instructions  in  ela- 
tion (5:37). 

Using  the  elements  listed  in  (6:6),  the  frequency  equation 
decomposes  into  the  following  product  of  determinants  in  the 
limit  as 

Limli  deil  C^l  = 0,-Di  »0, 


(6:7) 


s"  s"  =" 


Ill 


(6:9) 


112 


Using  (6:9) 


( 5 ! 24 ) gives 


u‘;’.u'“  = o. 


(6tl0) 


s d"’  (1  = 1,2,3)  are  i 


lated  to  strains  only  through  tho  d“' Ci  = 1,2,3). 

On  the  other  hand,  if  Dj^-O,  then  a“’,  b“’ (I .=  1,2,3; 
j=2,3)  are  nonsero,  and  (5:23)  yields 


t’"  ■ Af  . 


(6:11) 


folloving  die- 


uf.o. 

Dj.O  l.  tl. 


C6:12) 
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Moreover,  alt  of  the  elastic  constants  except  and  are 
involved  in  the  determinant.  This  is  to  be  expected  since 
the  stress-strain  relations  in  Chapter  II  show  that  when 
ej’  = 0 the  elastic  constants  c‘j|j(=  ) and  Cj^Cet 

play  no  r&le  in  the  determination  of  stresses. 


==44°3  > 


3.  Frequency  Equation  for  Axlsynnat 


/eling  »av 


To  eliminate  the  angular  coordinate  from  the  dispiece- 
raent  equations  the  integer  n may  be  set  equal  to  zero.  When 
this  is  accomplished  certain  simplifications  accrue  to  the 
formulation  of  Che  frequency  equation.  Using  equation  (5:37), 
Che  elements  of  the  frequency  equation  deteminant  can  im- 
mediately be  written  as 

c,4*  0, 

(6:13) 

= 0. 

C,4= 

Ciq  - - 


H4 


c,,. 

c„  * -c4S4''U"?/s^-'t 
c,*  - -Cdf  ^ dl") 

tis  ' ^ic  --  0- 
Ce."  <^L'- 
c,j,» 

Cet-' 


c,..  ^tId;''<;v“  - 4" 

C^‘  -4" 


(5:37). 
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UW.l  D,'D,  = 0 
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With  regard  to  equation  (3:57)  the  displacement  field 
resulting  from  the  requirement  chat  Dj  = 0 is  given  as 


(6:17) 

BjW, r)]\;i"'sm(tl-ni). 


Thus,  the  frequency  equation  Dj  = 0 corresponds  to  the 
family  of  axiaym.'ieeric  modes  of  traveling  vaves  generated 
solely  by  Che  f'^(f)  (i'=l,2,3)  potentials.  The  displacement 
field  is  torsional  in  character  end  has  nodal  cylinders  for 
^ sufficiently  large.  Consideration  of  the  elements  in 
shows  a dependence  on  both  d^  and  dg  which  Suggests  the 
isechanisms  of  tangential  and  radial  shear  employed  to  resist 
torsion.  When  the  limit  is  taken  to  isotropic  conditions, 
d^-»  dg,  and  the  role  of  yi'''(cdg'')  is  clearly  revealed  to  de- 
rive from  die  two  separate  anlsocroplc  mechanisms. 

C.  = 0 is  the  frequency  equation  for  nontorsional  axl- 
ayrotieCric  wave  motion.  Only  the  integration  constants  Aj  and 
(is  1,2,3)  are  tero,  consequently,  equation  (5:23)  yields 


+ (6:18) 

f"\n‘o. 

Using  equation  (5:26),  the  displacement  field  Is  given 
^ (6:19) 
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The  displacemont  field  is  seen  to  be  comprised  of  redial 
and  axial  components  end  is  caused  separately  of  the  f^’^Cr) 
(ie  1,2,3)  potentials.  In  the  limit  to  isotropy  the  two  sets 
of  potentials  and  ^ (i  = l,2,3)  pass  over  to  equivolumi- 
nal  and  dilatational  potentials,  respectively,  and  Che  fre- 
quency equation  dependence  on  elastic  properties  can  be  re- 
lated CO  the  three  Poisson's  ratios  v“’ (i  «=  1,2,3).  In  Che 
anisotropic  formulation  all  of  the  elastic  moduli  are  pre- 
sent, indicating  an  involved  dependence  of  Che  modes  of  wave 
notion  on  the  elastic  parameters  of  the  layers  of  the  cylin- 
der. This  result  is  to  be  expected  on  the  basis  of  the 
transversely  Isotropic  stress-strain  relations. 

4.  Frequency  Equation  for  Axls'nrcrietrle  Haves  of  Infinite 

If  both  n and  ere  set  equal  to  zero  the  frequency 
equation  degenerates  to  a product  of  three  6X6  determinants , 
reflecting  the  fact  that  Che  three  sets  of  potentials  (fp 
ipj',  and  f^''  (t  = 1,2,3)  generate  uncoupled  motion  in  the  com- 
posite cylinder.  The  appropriate  fomrulaCion  of  the  elements 
of  the  frequency  equation  may  be  derived  from  (5:37),  or 
either  of  the  last  two  sections  may  be  used  as  a starting 
point.  Choosing  to  proceed  with  the  formulation  established 
in  Section  2,  Che  elements  are  given  as 


C,,'  = 


(6:20) 
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<:.4  ' 

Cm  ‘^<'0, 

c,.^  = - t 

!^i4-  - 

Cji ' -<^4"  i\‘'’.a  JiCf^doi. 

c,4=  -4"  R‘-‘’Q\«.p'"a), 

(6:20) 

c„--c,^=.o, 

C,J«  -‘(t 
^•>4’ 

C,4  ' C?4  » 0. 

Ce5-=  <"t  j,(ij'’6l, 

C*t=  sr5Y,<s.'-’bl, 

C»> 

c,^'  V„{p'"6), 


Dj^  = 0 is,  thus,  the  frequency  equation  fc 
gitudinal-shear  vibrations , and  the  dlsplec 
ends  only  on  (f)  (i  = 1,2,3). 

If  03  = 0 only  A3’,  B‘3'  (i  = l,2,3)  need  be  r 
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and 


U.'’'  : = 0 , 


(b) 

(c) 


(6:24) 


It  is  clear  that  limit  hence  Dj  “ 0 yields  the 

cut-off  frequencies  of  the  axlsymmetrlc  torsional  waves. 
Since  Dj = 0 produces  the  plane-strain  displacement  field 
(6:24),  it  follows  that  the  torsional  waves  approach  plane 
strain-shear  vibrations  as^---0. 

Finally,  with  Dj  = 0 only  Aj’  and  (i  = l,2,3)  need  be 
nonzero,  yielding 

Hi”  - , tpl'’ = M 


(b) 


(6:25) 


u‘"  ■»  (A"’j.^‘^f)  + Si\(<gr))c.osat  . (e> 

Cj = 0 is  the  frequency  equation  for  symmetric,  plane- 
strain  extenslonal  vibrations  in  the  radial  direction.  In  the 
limit  to  isotropy,  (^2^1  = 1,2,3}  becomes  the  dilatational  po- 
tential, hence  the  vibrations  represented  by  (6:25)  are  dila- 
tational in  the  Isotropic  state. 

It  is  to  be  noticed  that  the  symmetric  plane  strain  fre- 
quency equation  V>2  degenerates  to  Dj  as  n-»0,  uncoupling 
the  f‘‘  and  £{‘2  (i  = l,2,3)  potentials.  For  isotropy  this  corre- 
sponds to  a dlsassociation  of  dilatational  and  equivoluminal 
motion.  Furthermore,  as  ^■•0,  the  nontorsional  frequency 
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equation  degenerates  to  Dj  and  indicating  that  the 
axisyirmetric  nontorsional  traveling  waves  eminate  from  the 
uncoupled  radial-extens ional  and  longitudinal-shear  vibra- 
tion fields. 

Two  broad  distinctions  in  the  behavior  of  the  frequency 
equation  are  discemable.  For  the  ease  of  axisymmetty  there 
are  three  families  of  cut-off  frequencies  which  are  derivable 
by  solving  the  three  transcendental  equations  Dj^cO,  D^cO, 
and  As  is  increased  from  zero  the  radlal-extcnsional 

end  longitudinal-shear  vibration  fields  couple  and  two  fami- 
lies of  branches  of  the  frequency  equation  exist  correspond- 
ing to  the  solutions  of  Dj  = 0 and  that  is,  torsional 

two  families  of  cut-off  frequencies.  These  derive  from  Dj^=0 
and  1^2°^  correspond  to  longitudinal  thickness-shear  and 
plane-strain  vibration  fields,  respectively.  As  '<(  is  in- 
creased from  zero  these  two  fields  couple  resulting  In  tra- 
veling waves  for  which  all  components  of  the  displacement 
field  will  be  generally  in  evidence.  There  is  Just  one  fam- 
ily of  branches  of  the  frequency  equation  corresponding  to  the 
solutions  of  equation  (5:34). 

5.  Frequency  Fouation  for  One  and  T;?o  Layered  Shells 

It  is  readily  demonstrable  that  when  two  layers  of  the 
cylinder  are  composed  of  identical  materials  the  frequency 
equation  degenerates  to  a product  of  two  determinants  - one 
of  which  is  the  frequency  equation  for  a two  layered,  traction 


f e ^ S' 


Layer  material 


E: 


aeilCijl  ■ r.-F;  ^ a 


c.i  c„  c,^  c„  c„ 
c.^  c„  c„  c,,  c,. 


n F,  =0  is  the  f 
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isotropic  or  isotropic,  respectively,  if  equations  (5:37) 
or  (6:3)  are  used  for  the  elesicnts,  ?or  the  Isotropic 
case,  Armenakas  [29]  has  obtained  this  equation  in  different 
form.  The  two  frequency  equations  coalesce  if  Fj^  is  altered 
by  replacing  columns  I,  2,  7,  and  8,  respectively,  by  the 
difference  between  columns  1 and  5,  2 and  6,  7 and  11,  and 
6 and  12 > The  differences  encountered  between  the  two  for- 
mulations arises  from  the  initial  choice  of  potentials  used 
to  separate  the  equations  of  motion. 

Repeating  the  process  Just  used  on  dctlC^j  I , but  this 
time  operating  on  yields 

Limil  F,  = Y,.Fi-0,  (6:28) 

where  Y2  is  of  the  same  nature  as  and  F2  is  given  by 


C„  C,i  C„  C„  C,5  C,^ 

Cii  C,3  C„  C„  C„ 


‘'•VJ  ^'•■4 

The  equation  Fj = 0 is  the  frequency  equation  of  an  iso- 
tropic or  anisotropic  homogeneous  traction  free  cylinder, 
respectively,  if  the  elements  are  taken  from  either  (5:37) 
or  (6:3).  Gazis  [2l]  and  Mirsl^  [39],  respectively,  ob- 
tained this  equation  for  homogeneous  isotropic  and  trans- 
versely isotropic  cylinders  by  solving  the  one  layer  cylin- 
der problem  directly. 
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6.  Layered  Cylinders  with  Clamoed-Froe  and  Clamped- 

In  Chapter  V the  boundary  and  Interface  conditions 
utilized  to  obtain  the  frequency  equation  for  layered  cylin- 
ders with  free  surfaces  are  given  in  equations  (5:27)  and 
(5:28).  It  is  a simple  and  straight  forward  problem,  based 
on  the  substantial  formulation  of  equations  (S:33),  to  re- 
write (5:27)  and  (5:28)  to  reflect  any  other  boundary  and 
interface  conditions  exprossible  in  tenrs  of  displacements 
and  stresses,  and  thence  obtain  directly  the  appropriate 
frequency  equation  to  replace  (5:34),  It  is  emphasized 
that  no  additional  formulation  beyond  the  above  mentioned 
is  necessary. 

The  three  displacement  components  at  a clamped  surface 
must  vanish  in  contrast  to  the  traction  free  surface  re- 
quirement. Thus,  recalling  the  construction  of  the  fre- 
quency equation  (see  (5:35)),  it  follows  that  if  a two  lay- 
ered cylinder  with  clamped  outer  surface  were  considered  di- 
rectly, a 12X12  frequency  equation  determinant  would  result. 
The  determinant  would  be  composed  of  the  listed  boundary  and 
interface  conditions  at  f = E end  r = b,  and  the  clamped  con- 
dition at  f = c,  which  has  the  form 


(6:30) 


Similarly,  if  the  inner  surface  of  a two  layered  shell 
nant  arises 


is  clamped, 


12X12  determit 


for  Che  frequency 
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equation  uhlch  can  be  considered  as  composed  of  the  conditions 


and  the  rema Inins  equations  of  (5:35)  at  r = c and  f = 3. 

Finally,  for  a homogeneous  cylinder  clamped  at  both 
surfaces , the  frequency  equation  is  a 6X6  determinant  com- 
posed only  of  the  conditions  listed  in  (6:30)  and  (6:31). 

In  the  limit  as  the  elastic  rigidity  of  the  Inner  lay- 
er approaches  infinity,  the  frequency  equation  for  a three 
layered  shell  uith  traction  free  surfaces  approaches  that  of 
a two  layered  cylinder  with  inner  surface  clamped  and  outer 
surface  free.  Alternatively,  when  Che  rigidity  of  the  outer 
layer  approaches  infinity,  the  frequency  equation  will  ap- 
proach chat  of  a two  layered  shell  with  traction  free  inner 
surface  and  clamped  outer  surface.  However,  for  8x183^11- 
metrlc  (n  = 0)  motion,  in  the  case  of  infinite  axial  wave- 
length, the  resulting  equations  will  include  the  effect  of 
axial  or  rotational  rigid  body  motion  of  the  outer  or  inner 
rigid  layer.  In  order  to  complete  the  clamped  requirement, 
it  is  essential  to  secure  Che  rigid  surface  by  an  external 
independent  force.  This  is  comparable  to  supposing  that  the 
rigid  layer  has  infinite  mass.  Consequently,  in  Chose  cases 
for  which  both  n and  '5  are  zero,  it  will  be  required  that 
both  the  rigidities  and  the  density  of  Che  inner  or  outer 
layer  pass  to  infinity. 

To  form  the  frequency  equation  for  a two  layered 
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cylinder  with  clamped  timer  surface,  the  rigidity  of  the  in- 
ner layer  of  the  three  layered  cylinder  will  be  llraited  to 
infinity.  Since  a rigid  material  is  capable  of  no  elastic 
response,  there  is  no  loss  of  generality  imposed  by  assum- 
ing that  the  inner  layer  is  initially  isotropic.  The  ro- 
oalning  layers,  however,  are  assumed  to  be  anisotropic.  In 
the  limit,  as  the  stiffness  of  the  inner  layer  goes  to  in- 
finity, b*''  reduces  to 


whore  the  ratio  remains  finite.  Horeo 

tlty  t reduces  to 


Thus,  for  finite  frequeneiesA,  it  follows  from  equations 
(6:1)  that 


(6:32) 


The  first  six  columns  of  Che  18X18  frequency  determi- 
nant, composed  of  Che  elements  listed  in  (6:3),  become  in- 
finite as  dg^eo.  To  retain  finiteness  of  the  determinant 
when  the  limit  is  taken,  these  columns  are  each  divided  by 
dg\  Distinguishing  the  resulting  elements  of  Che  frequency 
equation  by  equation  (5:34)  degenerates  Into  Che  pro- 
duct of  two  determinants  X and  E.. 

Limit  dellCyl  = Jt-E.,  = o.  (6:33) 
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E 


C„  C„  C,„  C,„  C,„  0 
C„  C,,„  C,„  C^„  0 


‘'.6  ‘'J.®  ^1,11  ° 


0 0 0 0 
0 0 0 0 
0 0 0 0 


E^  is  the  frequency  equation  for  a two  layered  cylin- 
der with  clamped  inner  surface  (f  «=b)  and  traction  free  out- 
er surface  (r»d). 

X is  a nonzero  6X6  determinant  comprised  of  the  ele- 
ments of  the  first  six  rows  and  columns  of  detlc^jl.  Since 
0 does  not  appear  in  X,  this  determinant  is  not  a frequency 


equation. 

For  the  case  of  axisymmetric  traveling  waves , use  of 
equation  (6:32)  with  n = 0 causes  the  frequency  determinant 
to  decompose  into  the  product  of  the  following  three  deter- 
minants, as  d^'^co; 

Umit  delRijl  = n=o,  (6:3S) 
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Again  Xj^,  being  a function  of  ^ alone,  does  not  repre- 
sent a frequency  determinant.  Hie  equation  £2=0  is  the 
frequency  equation  for  torsional  waves  in  an  infinite  cylin- 
der with  clawped  inner  surface  (f  = B)  and  traction  free 
outer  surface  (r=3).  The  equation  Ej  = 0 1s  the  frequency 
equation  for  nontorsional  traveling  waves  in  a clanped-free 
cylinder.  From  the  element  listing  of  equation  (5:37)  it 
can  be  seen  that  £2  depends  only  on  the  potentials  f’^’Cr) 

(1  = 2,3)  and  Ej  only  on  (pj(r)  (i  = 2,3;  J =1,2).  Thus,  the 
limit  in  equation  (6:35)  causes  an  uncoupling  of  tho  combined 
effects  of  the  potentials. 

It  is  readily  apparent  that  if  instead  of  the  inner  lay- 
er, the  outer  layer  is  limited  to  infinite  rigidity,  the  fre- 
quency equation  will  again  degenerate  in  a similar  manner. 
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the  outer  (isotropic)  layer 


. by  equation  (6:32)  with  1 


numbers  are 


(6:37) 


Limit  - CO. 


.to  Cetlc,^l  vithn.Oand 


U^Wetlc,h 


E4 


c„  C,,  C„  C,,  C,, 

C.c  c„  c„  c„  c,, 

C„  c„  C„  C„  C,,  0^ 

c*.  C„  c„  c„  c„ 

c„  c„ 

Cc,  C,,  c„  c„  c„ 

C„  c,,  c„  c„  c„  c„ 

C„  c,,  C„  C„  Cg, 

c,.  C„  c,,  c,.  c„  c,, 

0 0 0 0 0 0 


0 0 0 0 0 0 
0 0 0 0 0 0 


Cm  C«  C.,  C^„ 

c«  c„  C„  c,„ 

Cm  C,,  C„  C,..  C,,„ 

C.,  c,,  c„  c,„.  c,,„ 

C„  c„  c„  c,  „ c.„ 

C>,  c^  c„  c.,„  c,„, 

C4:  C,„  C,.  C„„ 

C,,,  c,,,,  C,,„. 

C,i,,  C„  , C„,  C^  .,  Cn„ 


5 ofA.  E^  = 0 is  the 


f f f f f f f 
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frequency  equation  of  a two  layered  shell  with  free  inner 
surface  (f  ca)  and  clamped  outer  surface  (f  = c}. 


When  n^O,  the  lijRit  gives 


The  equations  and  Eg  ■=  0 are,  respectively,  the 

frequency  equations  for  free  (raa)  - clamped  (r  = c)  shells 
for  the  cases  of  torsional  and  nontorsional  traveling  waves. 
As  before,  is  not  a frequency  equation. 

If  both  the  inner  and  outer  layers  of  the  shell  are 
treated  together  as  they  have  been  previously  on  an  indi- 
vidual basis  the  resulting  frequency  equation  will  be  for  a 
homogeneous  transversely  isotropic  cylinder  with  clamped 
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<r=B)  • 


clamped  (f  cc)  boundary  conditions,  limit 


Limtl  I | 
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C„, 
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C,i, 

C>v,, 

C.vi 

Ctt.ii 

Finally,  if  the  frequency  equation  (6:12)  is  written 
for  the  axisymmetric  case,  the  resulting  expression  is 

Limtl  detlCiil  o^O.  (6:44) 


(a) 


(6:45) 


(b) 


Eg  = 0 and  E^=0  are,  respectively,  the  torsional  end 
nontorsionsl  frequency  equations  for  a clamped-clamped  homo- 
geneous cylinder. 

In  the  limit  when  "V  approaches  0 the  axisyimietric  (n  = 0) 
frequency  equations  E2,  Eg,  Ej,  Eg,  Eg,  and  Eg  reduce  to 
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those  of  appropriately  clamped  two  or  one  layered  shells. 

It  should  be  pointed  out  that  in  the  limit  as  dg~<°  and/or 
dg-«o>,  the  general  frequency  equation  (6:21)  corresponding 
to  Che  axlsyrranetric  vibrations  of  three  layered  traction 
free  shells  does  not  reduce  Co  the  frequency  equation  of 
the  two  layered  clamped  shell.  As  was  indicated  previous- 
ly, these  equations  Include  Che  effect  of  rigid  body  mo- 
tion of  Che  outer  or  inner  rigid  layer.  To  obtain  the  ap- 
propriate frequency  equations  for  clamped  cuo  and  one  lay- 
ered shells  the  additional  limit  must  be  taken  as  and/or 
approach  infinity. 

7.  Freouenev  Eoiiations  for  Cylinders  with  Clamocd-Free 
and  Clamped -Clainoed  Boundary  Cond Itions  When 't  =n  » 0 

Consideration  of  equations  (6:23c)  and  (6:24c)  reveals, 
as  intimated  in  Section  6,  that  the  sxisymmecric  cut-off 
modes  of  vibration  can  allow  rigid  body  motion  of  a layer  of 
the  cylinder  if  the  rigidities  of  that  layer  are  indeed  in- 
finite, To  obtain  the  appropriate  frequency  equations  for 
clemped-free  and  clamped-clamped  boundary  conditions  the 
frequency  equations  of  Section  4 can  be  treated  in  the 
limit  as  first  dg  approaches  infinity  and  subsequently 
where  i is  1 or  3,  or  both,  depending  on  the  choice  of 
boundary  conditions  sought.  It  should  be  recognized  that, 
regarding  Dj  which  is  for  radial  extenslonal  modes,  it  is 
not  necessary  to  limit  ^'‘ot  the  rigid  layer(s)  to  infinity. 
That  a rigid  layer  can  not 


This  follows  from  the  fact 
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Assuming  thst  Che  Inner  (isotropic)  layer  is  limited 
first,  the  radial  wave  numbers,  for  finite  frequency  and 
density  become 


Limit  is’  , »0.  (6:46 

Applying  (6:46)  to  D^,  Dj,  and  Dj,  and  subsequently 
string  pass  to  infinity  reduces  the  frequency  equation 
5 the  following  product 

Limit  J Limit  ^5,- D,  I Djj  = E.  „ 'E-i  Eji=  0,  (6:47 


The  elements  of  the  above  frequency  equations  are 
listed  in  equation  (6:20).  Each  of  the  three  equations 
represents  the  frequency  equation  of  a two  layered  cylin- 
der clamped  at  r = b and  free  at  f = d.  The  appropriate 
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displacement  field  equations  for  E^2>  ^32  listed 

in  equations  (6:23),  (6:24),  and  (6:25). 

If  the  limiting  process  discussed  above  is  executed  on 
the  outer  layer,  the  follouing  product  of  frequency  equa- 
tions is  obtained. 


f the  three  equations 
two  layered  shell  clomped  at  r 


3 Che  frequency  cqtia- 
3nd  free  at 

rc=a.  The  elements  of  those  determinants  are  listed  in  c- 
quation  (6:20).  The  appropriate  displacement  field  equa- 
tions for  Sg|_i  and  Egj  ate  listed  in  equations  (6:23), 
(6:24),  and  (6:35)  with  (1  = 1,2). 

Finally,  enacting  the  above  limits  in  concert  produces 


135 


LitnK|  Limit  ^D,  D,  j DjJ  I E-»i  • ' E-,i  , (l'  i,«,  (6:51) 


These  equations,  representing  the  most  complete  reduc- 
tion in  complexity  of  the  general  frequency  equation  (5:34), 
are  the  frequency  equations  for  homogeneous  cylinders 
clamped  at  r = 5 and  f = c.  Again  equations  (6:23),  (6:24), 
and  (6:25),  respectively,  apply  to  Eg,  and  vith 

i restricted  to  the  single  value  of  2.  baltrukonis,  et  ^ 
^27,78]]  has  treated  many  clamped-free  cases  for  homogeneous 
cylinders  and  rods. 


FreouencY  Eouati.on  for  Cvlii 


An  important  problem  in  modem  technology  Is  Che  design 
and  analysis  of  composite  shells  with  thin  facings  (or  outer 
layers).  Such  construction  is,  indeed,  typical  of  light 
weight  sandwich  shells ■ 

It  is  well  knovm  that,  in  the  strict  sense  of  the 
phrase,  Bessel  functions  do  not  possess  addition  formulas 
[66j,  that  is,  it  is  impossible  to  write  Zjj{X  + x)  as  an  al- 
gebraic function  of  Zj^(X)  and  Zj^fx).  There  are,  however, 
generalized  formulas  which  pass  under  the  name  of  addl- 
cion  formulas  and  which,  in  the  case  Chat  (for  specific- 
ness) x«X,  can  be  used  to  provide  some  simplification  of 
Che  frequency  equation.  It  must  be  borne  in  mind  that  the 
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following  f onr(ulation&  are  a departure  from  the  exact  Bolu- 
tlons  presented  earlier  in  this  chapter  and  in  Chapter  V. 

From  a theorem  due  to  deunann  and  subsequently  gene> 
ralized  by  Graf  [66],  the  following  addition  formulas  for 
regular  and  modified  Bessel  functions  of  the  first  and  sec- 
ond kind  may  be  written; 

(a) 


= Ej->r 


(b) 


and  the  numbers  x and  X are  generally  complex, u 
trary,  and  the  angle  if'  is  defined  by 

X-Titos6  = 6 to^'f , isint^asin^. 


d = X--J  , t/<  *0, 
and  when  e = 

Using  (6:54)  in  (6:53),  with  x«X  a 
the  integer  n,  produces 

= (,  t^]z„(X)  (X], 


(a) 

(6:54) 

(b) 

3 replaced  by 


(a)  (6:55) 
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and 

= (l  ^ t (b)  (G:53j 

In  obtaining  (6:55),  recursion  fonKilos  (previously 
used)  and  the  additional  facts 
J„(X)  =(-irO„(Xl, 

1^(X)=  1„(X), 

wore  employed  in  conjunction  with  the  assumptions  that  x 
is  BUtficiently  smaller  than  both  X and  1 so  that  J„(»:) 
and  7„(*)  (n  = 2,3,.,.)  are  negligably  small  relative  to 
Jj^(x)  and  Ij^(x);  and  Jq(>;),  Jj^(x),  IqCx),  and  Ij^Cx)  may  be 
replaced  by  their  equivalent  small  argument  asymptotic 


(see 


If  the  three  layered  cylinder  has  than 
Figure  5)  then 


(6:56) 


The  approxir.’ations  (6:55),  in  light  of  (6:56),  are  va- 
lid in  the  low  frequency,  long  wavelength  region  of  the  fre- 
quency spectrum  and  co;isoquently  the  appropriate  choice  of 
Bessel  functions  are  initially  and  Hovjever,  since 

the  approx SjTiation  does  not  lose  validity  as  becomes  large, 
and  since  the  turning  point  lines  for  the  facings  arc  steep 
(see  equation  (5:61)),  the  goneral  formilations  of  (6:55) 
will  be  preserved  so  that  sectors  other  than  7 are  properly 
described. 


Equations  (5:37)  must  be  altered 


reflect  (6:55), 


Obviously  only 


c„^  • y^")z..C’5i'»'il-  . 
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These  results  apply  to  shells  ulth  conparatlvely  large 
relative  thiclcnesses , for  the  assumptions  made  do  not  limit 
the  thickness  of  the  core  layer.  In  the  extreme  case  that 
and  vanish,  equations  (6:59)  vanish  and  (6:5S)  be.* 
comes  exactly  the  frequency  equation  for  torsional  axisyia- 
metric  waves  in  a homogeneous,  traction  free  cylinder. 

Similar  treatment  of  the  frequency  equation  D.  = 0 will 
yield  corresponding  approximate  frequency  equations  for 
thin  layered  shells.  However,  for  this  case,  the  results 
are  rather  involved  and  consequently  will  not  be  discussed 
here. 

9.  N-lavered  Rods  and  Cylindrical  Shells 

To  generalize  tho  theory  of  Chapter  V to  describe 
multi-layered  shells,  it  is  only  necessary  to  Initially  sup- 
pose the  cylinder  to  be  composed  of  N layers,  and  to  write 
additional  interface  conditions  at  each  now  layer  interface. 
The  free  (or  clan5>ed)  conditions  at  the  inner  and  outer  sur- 
faces remain  as  previously  written.  Thus  for  an  N-layered 
shell  there  will  correspond  a 6HX6N  frequency  determinant. 

Any  attempt  to  manipulate  the  given  frequency  equation 
(5:34)  In  a manner  that  will  give  way  to  the  proper  descrip- 
tion for  a three  layered  rod  meets  with  only  partial  success. 
It  is  much  more  straight  forward  to  review  the  boundary  con- 
ditions and  to  observe  that  when  the  inner  layer  attains  the 
geometry  of  a solid  rod,  the  Inner  surface  bounda 
tlons  are  superfluous.  Furthermore,  since  the  Be 


sry  cc 
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functions  of  the  second  kind  approach  infinity  as  the  argu- 
ment approaches  zero,  the  inner  layer  solutions  must  be 
modified  by  setting  bV’C1c1,2,3}  to  zero.  This  assures 
that  a properly  obtained  solution  for  the  displacement  vlll 
be  finite  at  the  origin.  Thus , In  urltlng  the  frequency 
equation  for  a three  layered  rod  it  Is  necessary  to  doleto  the 
first  three  rows  end  the  second,  fourth,  and  sixth  coliimns 
of  the  cylinder  frequency  equation  previously  obtained.  The 
extension  of  the  rod  frequency  equation  to  N-layers  is  ob- 
vious and  the  corresponding  frequency  determinant  Hill  be 
(6M-3)X(6N-3). 

The  frequency  equations  generalized  as  Indicated  above 
vlll  degenerate  to  lesser  order  determinants  analogously  to 
the  degeneracies  previously  described  for  the  cases  of  n =0, 

? = 0,  and  n ■=■?  =0. 


CHAPTER  VII 


In  Chapter  VI  it  was  shown  that  all  branches  of  the 
frequency  equation  initiate,  when  V with  some  zero  or 
nonzero  cut-off  frequency.  The  axisymnetric  cut-off  fre- 
quencies are  determinable  by  numerical  analysis  of  the 
three  transcendental  equations  given  in  (6:32)  and  the  non- 
axisymmetric  cut-off  frequencies  follow  from  the  transcen- 
dental relations  listed  in  equation  (6:17).  These  five 
equations  are  complicated  Implicit  relationships  in  O with 
parametric  dependence  on  the  elastic  constants  and  Che  shell 
geometry.  Ko  analytical  aCtacU  will  be  made  on  these  equa- 
tions since  simplifications  can  only  be  expected  sftien  the 
Bessel  functions  are  representable  in  asymptotic  series  for 
either  very  large  or  very  small  fl.  Since  there  is  interest 
in  the  branches  of  the  frequency  equation  with  intermediate 
cut-off  frequencies,  a numerical  analysis  will  be  applied 
to  obtain  all  cut-off  frequencies  sought. 

When  the  axial  wave  number  is  not  zero  the  branches  of 
the  frequency  equation  may  bo  ascertained  by  solution  of  the 
two  equations  listed  in  (6:26)  for  axisymmetry , or  equation 
(5:34)  for  non-axlssTHmctry.  These  equations  involve!!, 
and  n implicitly,  with  parametric  dependence  on  the  elastic 
constants  and  shell  geometry.  Again  a numerical  solution 
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wilt  be  employed. 

If  the  slope  of  a straight  line  emanating  from  the  ori- 
gin of  the  0.-'?  plane  is  designated  by  m, 

(7:1) 

and  this  expression  for  Cl  is  used  in  the  frequency  equa- 
tion, then 

=0  (7:2) 

represents  the  requirement  which  values  of  ? must  satisfy 
for  each  m and  n in  order  for  waves  to  propagate  in  the  cy- 
linder. From  (7:1)  the  value  of  fl.  corresponding  to  each  value 
of  can  be  found.  If  the  slope  parameter  m is  varied  in  suf- 
ficiently small  increments,  the  resulting  CL-f  point  sets  be- 
come sufficiently  dense  on  the  frequency  spectrum  to  yield 
the  locus,  for  each  n,  of  Che  branches  of  Che  frequency  equa- 
tion. 

When  the  axial  wave  number  is  zero,  the  phase  velocity 
is  infinite.  As  ^ is  increased  from  zero  the  phase  velocity 
decreases  and,  on  physical  grounds,  should  approach  a finite, 
nonzero  limit  when  ^ approaches  infinity.  This  follows  from 
the  fact  chat  when  the  wavelength  is  sufficiently  small,  the 
medium  in  which  the  wave  is  traveling  takes  on  the  physical 
proportions  of  an  infinite  continuum  with  regard  to  the  mag- 
nitude of  the  wavelength.  Hindlin  and  Pao  [20]  have  shown 
for  homogeneous  rods,  and  Armenikas  and  Keck  [30]  for  two 
layered  rods,  that  all  branches  of  the  frequency  equation 
approach  asymptotically  to  specific  values  of  m as  ^ ap- 
proaches infinity.  The  limiting  values  of  m are  directly 
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related  to  phase  velocities  of  waves  propasating  in  infinite 
and  layered  eemi-inf inite  mediai  For  this  range  of  values 
of  Che  phase  end  group  velocities  of  the  waves  ere 

equal,  that  is,  the  waves  are  nondispersed . 

Along  a ray  of  slope  m the  equation  (7:2)  will  be  sat- 
isfied by  an  infinite  set  of  discreet  values  of  a con- 
tinuous set  of  values  for';,  or  no  ; at  all  except  ; ^0.  If 
the  values  of  ■?  which  satisfy  (7:2)  are  discreet,  then  the 
branches  of  the  frequency  equation  cross  the  ray  (7:1)  at  the 
values  of  t found.  If  the  values  of  ; are  continuous,  the 
branch(s)  run  along  the  ray.  If  there  are  no  values  of 
for  which  (7:2)  is  satisfied,  then  no  branch  of  the  frequency 
equation  can  progress  to  the  slope  m of  tho  ray  under  con- 
sideration. This  observation  forms  the  basis  of  the  asymp- 
totic analysis  of  the  frequency  equation  for  The  speci- 

fic values  of  m for  which  the  set  of  values  of  ¥ changes  from 
discreet  to  continuous  will  be  found,  and  it  will  bo  shown 
that  for  m less  than  these  specific  values,  the  left  aide  of 
equation  (7:2)  can  not  vanish. 

Initiation  of  an  analysis  for  large  axial  wave  number 
Is  beset  with  a confusing  array  of  conditions  related  to 
turning  point  lines  and  complex  zones.  Even  the  restric- 
tions of  equations  (5:4L),  (5:42),  and  (5:43)  do  not  alle- 
viate to  any  great  extent  tlie  dlveralty  of  conditions  which 
must  bo  embraced  by  the  analysis.  The  frequency  equation 
elements  change  from  regular  Bessel  functions  with  real  ar- 
guments to  modified  Bessel  functions  with  real  arguments 
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bade  to  resular  Bessel  functions  with  complex  arguments  ac- 
cording to  the  range  equations  (5:13),  (5:14),  and  (5:20). 
The  formulation  of  the  frequency  equation  is  the  same  in 
sector  7 for  all  of  the  anisotropies  under  consideration. 
Furthermore,  in  sector  1 the  formulation  Is  the  same  for 
all  less  anisotropic  and  isotropic  materials,  Hoxrever,  if 
a layer  of  the  cylinder  is  more  anisotropic  (either  class 
1 or  II)  there  will  be  a complex  zone  for  the  radial  wave 
numbers  p and  q of  this  layer  and  the  frequency  equation  is , 
consequently,  formulated  in  a different  manner. 

The  higher  sectors  are  of  little  Interest  In  the  analy- 
sis for  T-*®,  since  it  will  be  shown  that  all  brandies  of 
the  frequency  equation  progress  to  the  lowest  sectors.  It 
is  convenient  to  proceed  from  sector  1 to  higher  sectors , 
showing  first  the  conditions  under  which  branches  can  exist 
in  the  leanest  sectors  and  subsequently  that  the  branches  do 
Indeed  cross  the  higher  sectors  into  the  lowest  when  these 
conditions  are  met.  Hence,  the  analysis  will  be  conducted 
for  cylinders  with  all  layers  composed  of  less  anisotropic 
materials.  Under  this  restriction  the  sector  1 representa- 
tion of  the  frequency  equation  will  be  the  same  for  the  18 
cases  cited  in  equations  (5:41),  (5:42),  and  (5:43). 

1.  The  Froouenev  Eciuation  in  Sector  1 for  f --oe 

The  nine  wave  numbers  and  e“'  are  ell  imaginary 

and,  consequently,  the  elements  of  the  frequency  equation 


contain  only  modified  Bessel  functions.  Since  sect: 


BC  Of  b'‘’(i  = l,2,3)  (Figure  6),  ary  ray  used  in  this  sec- 
Lujii  (7:3) 


s (7:3)  and  (7:4)  i 


1 frequency  equation 
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■ ifSi-i"'"'"'  ■ '"f.- 


C„  » a-’B 


C,,«  (a‘'’.ari/5“’Be.’''’"'' 


0. 

Cs.-Cn' 


(7:5) 

(cent) 


Cfc  = 6'"Be''^“' 


‘ ' 0. 
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The  remaining  elements  of  the  frequency  equation  may  be 
formed  from  the  above  by  use  of  the  instructions  in  (5:37), 

It  is  significant  to  notice  that  equations  (7:5)  and 
(6:13)  have  exactly  the  same  nonzero  elements.  Thus,  it  fol- 
lows that  the  above  elements  will  yield  a product  of  two  de- 
terminants as  in  the  case  for  n>=0.  Further,  since  n does 
not  occur  in  equations  (7:5),  the  fortheomins  results  will 
be  valid  for  all  n. 

Using  (7:5)  in  (5:34)  gives 

L.mlt  = O,- D.  = O,  (7:6) 

where  Dj  and  are  arranged  from  the  elements  listed  in 
(6:25). 

Considering  first  the  possibility  that  = 0 yields 

, 0.  (7:7) 

The  first  term  of  (7:7),  incorporating  (5:20),  gives 

In  sector  1 this  equality  can  never  be  realized.  The  sarao 
Is  true  for  the  second  factor  of  (7:7).  Consequently,  re- 
gardless of  the  type  of  shell  under  consideration,  and  whe- 
ther n is  zero  or  not,  no  branch  of  the  frequency  equation 
for  torsional  motion  (or  coming  strictly  from  the  potential 
set  f“’(i»  1,2,3))  ever  penetrates  into  sector  1. 

Turning  to  D^,  using  (6:15b)  and  (7:5),  and  retaining 
only  the  largest  terms,  degenerates  to 
0,=  Dt’-D'""o5’"D4’'  =0, 


(7:8) 
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C^,  C„,  C,^„  C,,,^ 


The  determlnanta  Dg  (1°1,3)  represent  the  sane  fre- 
quency equation  but  at  the  free  surfaces  f aa  and  red,  re- 
spectively. Similarly,  the  two  determinants  Dy“  (iel,3) 
are  the  sane  except  that  Dy  is  for  r>b  and  Dj  for  fee. 
These  equations  will  be  discussed  separately. 


A.  The  Equation  =0  (1  = 1,3) 

For  isotropic  conditions  = 0 becomes 

pin.  I (7,j 

I *•;«“  -■?'»)  I 

This  is  the  frequency  equation  for  Rayleigh  surface 
waves  [47,79]. 

If  D^eO  is  satisfied,  only  the  integration  constants 
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and  B2' are  nonzero. 


Furthermore,  since 


-7K„.,CtF) 


the  displacement  field  equations  (S:24),  for  contain 

only  u“  and  u"’  components.  This  is  directly  comparable  to 
the  case  for  n = 0 and  was  to  be  expected  since  the  as3'mtot- 
ic  behavior  in  sector  1 is  independent  of  n. 

It  is  easy  to  shou  that  the  displacement  field  decays 
from  tho  f vg  free  surface  of  the  cylinder  in  the  case  of 
Dg  = 0,  or  from  the  free  surface  r = 3 if  Dg  = 0.  The  nature 
of  the  asymptotic  expansion  and  the  relative  size  of  ( is 
obtained  by  recalling  that  the  Rayleigh  surface  uaves  sat- 
isfy the  boundaiy  equations  and  equations  of  motion  uritten 
for  a single  materiol  half  space.  Thus  t must  be  sufficient- 
ly large  so  that  in  the  radial  direction  the  wave  motion  de- 
cays completely  within  the  thl^cness  of  the  inner  or  outer 
layer.  The  displacement  field  will  approximate  the  surface 
wave  phenomena  for  finite,  large  ; - being  a more  satisfac- 
tory approximation  for  a set  value  of  f in  thiclcer  layers 

Expanding  (7:10)  gives 
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In  this  form  it  is  clear  that  the  value  of  the  real 
roots  of  the  Ka/leigh  equation  are  dependent  on  Poisson's 
ratio.  It  can  be  shown  that 


as  v"' varies  from  0 to  %.  It  follows 
slope  bounded  by  the  limits 


lies  in  sector  1 when  1,  the  i-th  {i  ; 
cylinder  has  a displacement  field  which 
Rayleigh  surface  wave,  and  an  undisperse 
frequency  equation  lies  on  the  locus 


(a)  (7:12) 


(b)  (7:12) 
,3)  layer  of  the 


branch  of  the 


Notice  that  it  is  not  always  possible  for  this  locus 
to  be  in  sector  1.  Indeed,  if  0.871:'^ (1  s i ,3)  is  greater 
than  the  sector  L boundary  slope,  ma’l*'’  (i  = l,3)  can  never 
lie  in  sector  1.  This  restriction  docs  not  preclude  the 
existence  of  Rayleigh  waves  in  the  inner  or  outer  layers, 

It  merely  indicates  that  the  mechanism  involved  in  the  ap- 
pearance of  surface  wavos  is  not  so  simple  when  the  shell 
1s  fabricated  of  materials  placing  it  in  Che  sandwich  shell 
class  (5:41), 

Since  in  the  isotropic  case  - and  especially  for  core 
layered  shells  - it  is  possible  for  Dg'=0  to  have  a soluclor 
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in  Beceor  1,  ic  is  reasonable  to  expect  chat  there  will  also 
be  an  analogous  solution  for  at  least  soics  range  of  aniso> 
tropic  materials!  In  particular,  for  less  anisotropic  ma- 
terials, ic  is  reasonable  to  expect  a solution.  In  ever? 
term  of  the  frequency  equation  the  variables  (1  and  t appear 
only  as  the  quotient  Z,  hence.  If  a solution  exists  In  sec- 
tor 1,  regardless  of  the  degree  of  anisotropy,  ic  necessar- 
ily represents  an  undispersed  mode  of  wave  motion. 

Very  little  can  be  discovered  qualitatively  for  on  equa- 
tion even  as  simple  as  (7:10a)  \dien  all  the  terms  are  pro- 
perly represented  from  equation  (7:5).  Since  the  radial  wave 
numbers  appear  as  squares  and  first  powers,  little  can  be  ac- 
complished without  extensive  rationalization  • a process 
which  by  its  intrinsic  character  destroys  the  very  informa- 
tion boing  sought;  sector  1 values  of  m which  cause  (7:10s) 
to  vanish.  Numerical  analysis  will  have  to  be  employed  for 
the  anisotropic  case. 

B.  The  Equation  (1  = 1,3) 

As  was  done  with  the  Rayleigh  surface  wave  equation,  the 
discussion  will  commence  with  an  investigation  of  Che  isotrop- 
ic case.  To  avoid  needless  notation  the  superscript  (1)  will 
be  replaced  by  (1)  and  it  will  be  understood  that  similar  re- 
sults hold  if  (1)  is  replaced  by  (3).  The  equation  at  either 
interface  r=B  or  r = c is  identical,  the  only  difference 
being  in  the  set  of  integration  constants  which  must  be  non- 
trivial. 
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For  Isotropy  becomes 


(7:13) 


Equation  (7:13)  is  knovm  as  the  Stoneloy  vrave  equation. 
Stoneley  introduced  it  in  1924  [8l]  and  indicated,  rather 
summarily,  that  roots  exist  for  a restricted  range  of  values 
of  the  material  constants.  Since  then  Seholte  [82-87], 

Koppe  [88],  and  Sezawa  and  Kanai  [89-92]  have  discussed  it 
extensively.  A discussion  of  the  extent  of  variation  of  the 
material  constants  for  which  (7:13)  may  possess  a solution  in 
sector  1 will  be  pnsented  now  based  on  the  above  references 
and  the  author's  work  [30]. 

For  convenience,  the  following  parameter  will  be  defined 
p = (7:14) 

Ho  ambiguity  will  result  from  omitting  a superscript  on  p 
since  in  either  Dy”  or  Dy'”  only  one  of  the  dg’  (1  e 1,3) 
appears. 

To  proceed  with  the  discussion  a coordinate  plane  with 
ordinate  axis  p end  abscissa  b''’^  will  be  employed.  Roorrang- 
Ing  (7:13)  to  reflect  definition  (7:14)  yields 
ap^‘"  -p(i.0“’*)  2^"'  (i. -■■>)' 

, ip."'  (..s'"') 

Cl'"’. 


(7:15) 
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The  intent  will  be  to  plot,  on  a p,-b'®  plane,  the  bound- 
ary curves  which  determino  the  region  within  which  (7:15)  can 
be  satisfied  for  some  choice  of  m in  sector  1 of  the  Sl-5 
plane. 

If  b'''  is  the  turning  point  line  with  least  slope,  then 
the  slope  of  any  ray  considered  must  be  less  than  b'’.  But, 
b'“=i  by  definition.  If,  on  the  other  hand,  b'''  is  less  than 
b'*^,  then  no  ray  may  be  considered  which  has  slope  greater 
than  b'  . In  summary,  one  notes  that  the  perrnissablc  range 
of  values  for  the  ray  slope  m is 

0<mM  :{  (a) 


(7:16) 


C<m<b''’  (b) 

The  positive  axis  encompasses  all  possible  variations  of 
inequalities  (7: 16) . 

Expanding  (7:15),  yields 


This  form  of  the  Stoneley  frequency  equation  was  first  pre- 
sented by  Keck  and  Armenakas  [30]  and  is  conducive  to  in- 
vestigation on  the  p-b"'*  plane.  Expanding  the  «s''' and 
yields 
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t (l  -J  i-mVb™  J l-.-.V-i,<-‘  (7!17J 


e parametric  depend 


n ji  =0,  (7:17)  i 


= 0^  -m"’  ,Cp.  = o1  . 


(7:18) 
, to  v“'. 


thereyfl  = m??  is  stationary,  that  is,  tt=m‘“f  al«ay 
below  the  ray  Q =f  by  an  amount  bracketted  by  the  lit 
given  in  (7:12).  Thus,  for  p = 0,  only  if  (7:1 
holds,  and  this  value  of  m is  insector  1 of  the  fl-? 
only  if  b‘’’b  ra‘".  When  =m‘"  the  locus  of  the  bran 


1 given  by  n = ms't.  On  theji-b'’”  plane  the  1: 


d'-  = 0.  Which  is  satis: 

When,=».Dr  = 0 


s to  D'-^  =0,  V 


. that  m“b->  is  always  less  than  b“’  by  an  a 


(7:19) 
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equation  n = n'^b"'?  , lies  in  sector  1 of  Che  <1-T  plane.  On  the 
ji-o''  plana,  at  p the  line  parallel  to  the  b‘' ‘ axis  may 
be  thousht  of  as  the  locus  for  Che  solution  of  d‘^”  = 0,  which 
Is  satisfied  by  m ■=  mg!#’ . This  line  starts  at  b'”-  =0  and  con- 
tinues past  b'’"'  = 1 to  b'”^  = l/nia*"  which  is  the  upper  limit  of 
b'’'^  for  which  the  locus  remains  in  sector  1 of  cheG-^ 

plane.  At  this  upper  extent  the  lowest  turning  point  line 
and  the  lowest  branch  of  Che  frequency  equation  coincide. 

The  solutions  n = m‘,  (for  p = 0)  and  ni=in“b‘'’  (for  fi*'®) 
are  the  lowest  values  of  m in  sector  I for  which  Che  equa- 
tion =0  can  be  satisfied.  For  b"fl,  the  maxirauin  value 
of  ni  which  lies  in  sector  1 is  just  less  than  b‘'‘  itself. 

Thus,  if  m = b‘'’  is  substituted  into  (7:17)  the  resulting  equa- 
tion, When  plotted  on  the  p-b"'*  plane,  will  bound  the  region 
in  which  d'^’"=0  is  saClsfiable  by  m in  sector  1.  Performing 
this  substitution,  and  noting  that  1>*’ei,  equation  (7:17) 
gives 

J.-  \ <,  p.[3.(2.-kf"0  t 

t ^ (7:20) 

When  b“”’«=l,  equation  (7:20)  reduces  to 

p.»  V 1 = o. 

Thus,  the  quadratic  (7:20)  possesses  a double  root  p.»l  at 
b''**^  =1.  As  b"'*'  decreases,  the  two  root  values  of  p.  separate. 
When  b"’‘t=m®‘  the  lower  root  of  (7:20)  is  ji  = 0.  For  b"^  in?'' *• 
this  root  becomes  negative,  vihlch  is  unacceptable  on  physical 
grounds.  The  upper  root  alone  is  significant  for  b^'’‘  less 
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than  As  b'’'  -*  0,  (7:20)  becomes 

K (i  - 1 0 /TTvv^  -•»(!-  ib“VT."-’‘‘)]- 

which  Is  a quadratic  expression  in  jl  with  positive  root 


. 2./'^"^  ^ JT- 


Suimarizing  the  above  discussion,  Che  region  on  Che  {'•b'^' 
plane  for  which  Dj'’  = 0,  for  any  value  of  m in  sector  1 of 
the  ft-';  plane,  is  bounded  by  the  curves:  equation  (7:20), 
fl  = 0,  and  =1.  This  region  is  valid  only  when  the  low- 
est turning  point  line  is  Zsb  , which,  in  turn,  la  Che  low- 
est only  for  b * b =1.  In  the  region  where  Dj  =0,  the 
value  of  m (for  any  permlssable  choice  of  (1  and  b""- ) is 
greater  than  if  ji«  1 and  greater  than  m“6'’  if  ji  » 1. 

Figure  7 illustrates  these  results. 

The  discussion  for  the  alternate  case,  b*'*^*  1,  goes 
rapidly.  The  b'"‘'  axis  is  a boundary  curve  since  along  this 
axis  (flcO)  is  always  satisfied  by  njfwhich,  in  turn, 

always  lies  in  sector  1 when  b*'^  is  the  lowest  turning  point 
line.  When  the  line  parallel  to  the  b‘”‘  axis  extend- 

ing from  b*^^  = 1 to  b'"^  = l/m'**'  is  an  upper  boundary  since 
is  satisfied  if  ni  = iii"b'",  and  this  choice  of  m lies 
in  sector  1 of  the  Cl-t  plane  so  long  as  b“’*  i 1/m”*.  The 
curve  b = 1 is  a third  boundary  of  the  region  being  sought. 

The  remaining  boundary  of  the  region  in  question  is  ob- 
tained from  equation  (7:17)  by  using  the  substitution  m^cl. 
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This  follows  fran  the  fact  chat  so  long  as  mil,  equation 
(7:17),  if  satisfied,  is  satisfied  by  a value  of  m in  sec- 
tor 1 of  then-^  plane.  Using  this  value  of  m,  and  Che  fact 
that  1?’=1,  (7:17)  yields 


Equation  (7:22)  possesses  a double  root,  p = l,  at 
= 1.  When  b"”'  equals  l/mg'  one  root  of  (7:22)  is  jx-'m  , 
the  other  is  less  than  1.  For  b " > l/m^*,  only  the  lower 
root  retains  its  significance.  When  b“’*-»«  equation  (7:22) 


-i/b'"'  -dCi-  t ■. 

+ Cl-i/ib'->n  -4(1-  i/zb''*K>-'/a'U'’‘)]t 

which  gives  the  asynptote  for  (7:22)  when  b''’*—<o: 


(7:23) 


The  description  of  the  partition  of  the  jk-b'*''^  plane  is 
complete  and  Figure  7 shows  the  region  of  variables  for  which 
= 0 can  be  satisfied  by  a ray  with  slope  m lying  in  sec- 
tor 1. 

Two  contnents  concerning  Figure  7 are  in  order.  First, 
the  point  p = l,  1^'*  = I corresponds  to  the  case  for  single 
nuiterial.  Equation  (7:17)  vanishes  identically,  but  in  es- 
sence would  never  have  appeared  had  the  two  layers  been  1- 
dentical  and  also  perfectly  bonded.  Second,  the  line 
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b''"  =1  eorrssponds  to  Wieehert's  condition  for  Che  possibil- 
ity of  existenca  of  a solution  to  The  requirement 

b“"'  = 1 is  identical  to  This  is  the  condition 

vhich  Stoneley  [8l]]  iinplied  c?jsC  be  nearly  satisfied  for  a 
solution  of  the  determinant  to  exist.  Wiechert'a  condition 
vas  developed  as  a requirement  to  satisfy  certain  analytical 
models  for  earthquake  phenomena  in  layered  media.  From  Fig- 
ure 7 it  is  clear  that  conditions  very  divergent  from  the 
Wiechert  requirement  lead  to  a solution  of  d‘^  =0. 

Uhen  a degree  of  anisotropy  is  Introduced,  all  discus- 
sion must  be  based  on  equation  (7:9).  It  Is  again  reason- 
able to  suppose  that  lees  anisotropic  media  have  a material 
dependence  someuhat  akin  to  the  isotropic  case,  and  that  the 
anisotropic  version  of  D-j  possesses  a sector  1 solution. 

C.  Summary  of  Sector  i Limits 

For  sandwich  shells  both  b'”  and  Is”  are  greater  Chan  "b™. 
Considoratlon  of  normal  isotropic  aaterlals  shows  that  Ti'’  is 
roughly  1,5b”',  Thus,  the  Rayleigh  wave  loci  from  equations 
(7:9a  and  b)  can  not  oxist  In  sector  1.  Indeed,  if  Ts'’’  is  at 
least  l.lSlT'  then  m'Jb''’  (i  = l,3)  can  never  be  less  than  I?*®!, 
It  is,  however,  possible  chat  either  or  both  of  (7:9b  and  c) 
possess  a solution  which  lies  in  sector  1. 

Core  layered  sheila  present  the  opposite  condition. 

Both  b‘"  and  iJ*  are  smaller  than  !>*'.  They  are,  in  fact,  pre- 
sumed CO  be  less  than  0.86  sinceV'’  and  as  well,  are  less 
than  1.  In  this  cose  all  four  of  the  degenerate  cases  of 
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= 0 may  possess  sector  1 solutions  if  b*"  and  Ij”  are  not 

The  general  three  layered  shell,  vith  velocity  ratios 
characterized  by  inequalities  (5:43),  embraces  the  range  of 
values  on  the  plane  excluded  by  the  previously  dis- 

cussed shell  types.  The  Rayleigh-like  frequency  equations 
may  be  satisfied  by  sector  1 values  of  m if  the  value  of  b‘'' 
is  less  than  1 or  not  too  much  larger  than  1.  The  signifi- 
cant difference  in  the  behavior  of  the  frequency  equation 
is  in  the  fact  that  tho  Stonaloy-like  wave  equations  may 
be  satisfied,  for  a given  b‘"‘ , by  two  ranges  of  fi.  Figure 
7 shows  this.  As  a consequence  of  the  near  similarity  of 
Che  materials  of  the  shell  it  is  clear  chat  there  is  a 
much  greater  possibility  of  actually  realizing  the  Stoneley 
wave  loci  as  a part  of  the  asymptotic  behavior  of  the  fre- 

2.  The  Freouenev  Eouation  in  Sector  2 for 

Passage  from  sector  1 to  sector  2 is  a formal  process , 
enacted  by  simply  specifying  that  only  rays  of  slope  m lying 
in  sector  2 will  be  admitted  to  consideration.  The  frequency 
equation  trust  be  modified  to  a certain  extent.  A choice  must 
be  made  concerning  the  type  shell  under  consideration.  Equa- 
tions (5:41),  (5:42),  and  (5:43)  set  tho  bounds  for  sector 
2 for  ell  cases  under  consideration.  The  requireirent  to 
specify  tho  bounds  of  sector  2 is  of  importance  in  order  to 
decide  which  of  tho  g"  and  8^“’  (i  = 1,2,3)  change  from  imaginary 
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to  real  numbers. 

For  sandwich  shells  only  Bessel  functions  with  and 
s'“  in  their  arguments  are  affected  by  the  transition  to  sec- 
tor 2.  For  any  m > b"’,  the  Bessel  functions  involving  these 
two  radial  wave  numbers  are  always  regular  Bessel  functions. 

For  core  layered  shells  sector  2 may  be  bounded  by  the 
turning  point  lines 


(7:24) 


For  each  of  these  four  cases  only  the  radial  wave  numbers  g“’ 
and  s'”,  for  which  the  superscript  (i)  is  the  same  as  the  su- 
perscript of  the  lower  bound  of  the  range,  is  affected.  As 
long  as  m is  greater  than  the  lower  bound,  the  effected 
radial  wave  numbers  are  always  real. 

Finally,  for  general  three  layered  shells,  the  sector 
2 bounds  are 


The  remarks  following  (7:24)  hold  also  for  (7:25). 


(7:25) 


A.  Sandwich  Shells  - Sector  2 


Only  colirnn  elements  of  the  frequency  equation  determi- 
nant which  contain  Bessel  functions  in  p‘“  and  need  be 


(7:5). 


(7:27) 

V"'-  b'’'-?- . (b) 


and  12  (the  only  a 


••  «^58=  0.  <’’:28) 


C,,„=  - id‘'-'i«’‘-E=‘ 


f.. -4“].- r>. 
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C.,  = 

C.„=  -%™B  %>r><{‘«6,  (7:28) 

(cont) 

C.,.v'  6"’bcost?«B. 

c,, 

c,  - idf’ -4*’^;,  w?u»t“>6, 

c,,.  ■c,,„.o. 


n (7:28)  according  t 
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" rr(S,.-0-j4-m-(.-$.)(.-s;,)  ■ 

r 1 b s 


of  (7:29), 


will 
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repeatedly  pass  through  the  value  given  by  the  right  side. 
Thus,  on  any  given  ray  in  sector  2,  equation  (7:29)  possesses 
an  infinity  of  zeros,  each  demarking  the  passage  of  a branch 
of  the  frequency  equation  across  the  ray. 

Denoting  the  right  side  of  (7:29)  as^p,  chat  is,  Che 
principal  value  of  the  arctangent,  the  equation  can  be  «rit- 


When  J-.0 

(7:30) 

For  m > Che  branches  of  Dj  = 0 that  enter  sector  2 run 
CO  the  asymptote  m ■=  b . Thus,  all  of  the  torsional  modes  pro- 
pagate uith  the  speed  of  the  least  shear  wave  velocity  of  the 
three  materials  of  the  cylinder  when  the  wavelength  is  very 

Turning  now  to  D^cO,  the  folloving  equation  results 
when  only  the  dominant  terms  in  each  column  are  retained: 

0*"*’  DiOV’t)*  (7:31) 

where  Dg  (1»1,3)  arc  given  in  equations  (7:9)  and  Dg  is 
given  by  (7:32)  on  page  166, 

It  is  evident  that  the  passage  of  ^ into  the  real  do- 
main has  created  a displacement  field  in  the  center  layer  of 
the  sandwich  cylinder  which  is  not  entirely  a decay  field, 
hence  the  two  previously  obtained  Stoneley-like  wave  equations 
are  now  coupled.  Because  of  Che  trigonometric  torms  in  col- 
umns 7 and  6 (see  (7:28)),  Dg  = 0 possesses  an 


infinity  of 


i»l4»-J!-')-(4«.i“]* 


the  limit  «^'=0 
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t»05.  This  yields  the  locus 

a > (7:34) 

in  the  some  fsshion  ss  (7:30).  The  brsnches  of  both  and 
run  to  the  seme  locus. 

Reviewing  the  sandwich  cylinder  turning  point  line 
Inequalitlesi  inequalities  (5:41),  it  Is  possible  for 
(iol,3)  to  represent  nondispersed  mode  curves  In  sector  2. 

If  b“’— or  b™ — then  (1  = 1,3)  as  given  by  inequal* 

Ity  (7:12)  may  certainly  be  loss  than  and  V’  is  the  sec- 
tor 2 upper  boundary. 

The  very  notion  of  a nondisperslve  branch  of  the  fre- 
quency equation  lying  above  the  locus  of  an  infinite  family 
of  branch  curves  is  very  disagreeable.  It,  in  fact,  at 
once  confounds  the  physical  observation  that  the  Rayleigh 
waves  derive  from  the  lowest  brandies  of  the  frequency  equa- 
tion, and  the  requirement  that  branches  con  not  cross  one 
another. 

It  will  now  be  shown  that  the  Rayleigh  frequency  equa- 
tions in  the  product  (7:31)  may  represent  terraces  in  the 
high  order  branch  curves  of  the  frequency  equation  D^  = 0. 
Teinraees  are  described  as  zones  in  the  branch  curves  which 
have  the  approximate  slope  of  Rayleigh  wave  equation  loci 
when  (1  = 1,3)  and  thence  decrease  to  the  asymptot- 

ic limit  given  by  (7:34).  It  will  further  be  shown  that 
the  degenerate  frequency  equation  = 0 is  only  valid  when 
ID  is  approximately  equal  to  b'*’ and  '5'”i  . 

The  elements  In  D^,  fort>'l,  are  listed  in  equations 
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(7:5)  and  (7:28),  It  ia  to  be  noted  that  only  the  argu- 
ments of  the  exponential  or  trigonometric  terms  vary  with 
^ when  is  written  for  a particular  ray  with  slope  m in 
sector  2,  All  remaining  factors  ore  constant  and  bounded. 
Denoting  terms  which  are  constant  by  C,  terms  which  are 
exponentially  large  and  small  by  exp(x)  and  exp(-x) , re- 
spectively, and  trigonometric  terms  by  sin  and  cos, 

Slay  be  written  schematicaily  as  follows: 


Dg"  e"^^  0 0 0 0 


where  the  broken  line  encircles  the  elements  which  comprise 
Dg  and  the  ''terms"  outside  this  boundary  are  2X2  determi- 
nants which,  excepting  Dg’ (1  = 1,3)  indicate  only  the  order 
of  oiagnitude  of  each  such  determinant. 

Certainly,  when  Ixl— <s,  only  Dg'’,  Dg^  and  Dg  remain. 

The  underlined  elements  are  a sort  of  "soft"  coupling  be- 
tween the  three  diagonal  determinants.  If  D^,  given  above, 
is  expanded  for  the  purpose  of  grouping  order  of  magnitude 
terms,  and  only  terms  no  smaller  than  exp(-2x)  are 


retained, 
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that  Is  IX)  is  assumed  lorse  enough 
weakest  coupling,  becomes 

ft,  ‘0"ci.®[D,tZc*e-’sc»>  cos]  ■* 
■»  c‘  e"'*  Sirx  . 


to  exclude  all  Ixit  the 


(7:35) 


Any  value  of  ? and  m which  causes  to  vanish  locates 
a point  on  a branch  of  in  sector  2 of  the  frequency  spec- 
trum. Picking  such  a point  when  m is  epproxtniacely  iS*'  and 
Chen  varying  m,  so  chat  with  increasing  t,  is  always 

satisfied,  a branch  of  Che  frequency  equation  can  be  followed 
through  sector  2.  For  sandwich  shells  it  is  possible  Cltat 
mjb‘‘’  (i  = l,3)  both  lie  in  sector  2.  This  condition  is  as- 
sumed. It  is  asserted  that  as  f is  increased,  m decreases 
in  order  Co  obtain  points  on  the  same  branch  of  the  fre- 
quency equation.  As  m approaches  the  value  m^b"  (assuming 
mib‘'>  for  definiteness),  the  first  and  second  sums  of 

(7:35)  vanish  individually  due  to  tho  influence  of  D^’ and 
equation  (7:35)  reduces  to 


As  m is  further  decreased,  the  sum  JC®  sin  cos  e’^* 
varies  very  slowly  since  the  products  of  trigonometric 
terms  tend  to  counteract  each  other  and,  moreover,  e*^* 
makes  changes  in  magnitudes  nearly  imperceptible  when 
|x|»l.  Consequently,  as  m-m'Ji’'  (D^— 0),  equation  (7:35) 


reduces 
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Again  the  change  of  m as  ly  is  increased  is  possibly  very 
alo»  vhen  m is  nearly  equal  to  . As  is  increased 
still  further  finally  con  be  described  accurately  by 
Dq  alone,  since,  when  t-»»,  the  exponential  factors  vanish. 
Hence,  the  branch  runs  asymptotically  to  the  slope  iD  = b‘'’'. 

The  coupling  exhibited  in  (7:35)  can  never  be  ignored 
when  m is  not  very  close  to  b''^  since  it  is  possible  for 
either  of  the  last  two  sums  of  this  equation  to  play  the 
dominant  role  in  the  satisfaction  of  the  requirement  chat 
D^  = 0.  Thus,  the  weak  coupling  introduces  the  possibility 
of  terracing,  and  also  provides  the  conclusion  that  ° 0 
does  not  replace  the  requirement  D,;,  = 0 unless  m is  approx- 
imately b“'. 

B.  Core  Layered  Sholls  - Sector  2 

Inequalities  (7:24)  delineate  the  bounds  of  sector  2. 
Observe  that  a duality  in  the  four  conditions  exists.  The 
ensuing  discussion  will  illustrate  quite  clearly  that  only 
one  set,  for  example. 


need  be  considered  in  detail.  If  (7:36)  is  employed,  the 
sector  2 representation  of  the  frequent  equation  will  con- 
tain oscillatory  factors  in  the  elemancs  of  coluans  13,  14, 
17,  and  18.  If,  on  Che  other  hand,  the  second  two  of  (7:24) 
are  used,  Che  elements  of  columns  1,  2,  S,  and  6 will  be 
oscillatory.  Because  of  this  symnetry  of  the  frequency 
equation,  the  results  obtained  by  using  (7:36)  con  be  di- 
rectly reinterpreted  for  the  alternate  ease. 
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Equation  (7:5)  can  be  used  for  the  elements  of  the  fre- 
quency equation  that  do  not  change  when  moving  from  sector  1 
to  sector  2.  Using  the  requirement  of  (7:36),  the  elements 
of  columns  13,  14,  17,  and  18  are  obtained  from  (7:28)  as 
follows : 


placed  by  3. 

From  equations  (7:5)  and  (7:37)  it  is  apparent  that  the 
asymptotic  form  of  the  frequency  equation  again  degenerates 
to  the  product  The  determinant  yields 

The  first  factor  of  this  equation  has  no  solution  in  sector 
2,  and  the  second  factor  reduces  to 


For  isotropy,  equation  (7:38)  reduces  to  the  Love  wave 
frequency  equation.  In  the  usual  manner  the  branches  repre- 
sented by  Dj  = 0 in  sector  2 all  asymptote  to  the  sector 
boundary  given  as  m = t>'’ regardless  of  the  value  of  n: 


It  is  interesting  to  note  that  if  the  second  qualifica- 
tion of  (7:36)  3:8  employed,  then  in  sector  3 the  frequency 
equation  = 0 may  bo  written  as 


lb;  J = 13,14,1/,18J  = - Cj^  j 

,8,  . . . ,4j  j =7,8,11,12)  with  superscript  (2)  re- 


.1  = J ■=  13,W,L/,i»;  J 

: ia  \£l  17  1®^  w^qh  c re* 


(7:37) 


(7:38) 


(7:39) 


- 4“  i"  = 0. 


(7:40) 
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Both  the  inner  end  outer  layer  ere  seen  to  have  displacement 
fields  which  are  similar  to  the  Love  wave  field  and  the  core 
displays  a decay  displacement  field.  This  is  reasonable 
since  oao  may  picture  the  wave  front  in  tho  high  velocity 
coco  as  wheeling  or  tviming  to  the  outside  layers  and  thus 
concentrating  most  of  the  core  displacement  at  the  interfaces. 
= 0 takes  the  form 

Cl.'  • 0,  (7:41.) 

where  and  Dy'^'are  given  by  (7:9a  and  b),  and 


0 0 0 C,.,., 

0 0 C,5„  0 

As  concerns  Dg' and  only  wlicn  the  first  of  Condi- 

tions (7:36}  prevails  is  it  possible  that  the  effects  of 
these  determinants  are  not  felt  in  sector  2.  However,  if 
the  second  of  (7:36)  holds,  or  If  is  sufficiently  near 
I,*  in  magnitude,  Dg"  and  D^'"'’(if  the  materials  of  tho  inner 
layer  and  the  core  permit  the  existence  of  a Stoneley-like 
wave)  may  cause  terracing  of  the  branch  curves  in  sector 
2 at  the  values  of  m for  which  each  of  Dg  and  D^'“  vanish. 

The  frequency  equation  D^'g  = 0 reveals  the  locus  of 
all  branches  when  ra  is  sufficiently  close  to  b‘”.  This 
statement  relies  heavily  on  previous  arguments  and  is  sup- 
ported here  by  tho  fact  that  if  d'^q  is  expanded  the 
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resulting  terms  have  the  form  J ICCcos^^d  cob>'*c  + sin'^®d 
atnV'c)  or  * K(sln’{’’S  eos-f'’'c  - sinVc  coBi^”a).  Thus  each 
term  contains  either  sini^^{a-c)  or  cosl!’'(5  - c).  This  con- 
sequence forces  the  conclusion  that  for  any  ray  in  sector  2 
there  arc  an  infinite  number  of  values  of  ? for  which 
d'^qsO.  All  modes  must,  then,  asymptote  to  Che  sector  low- 
er bound  given  as 

-«=  . (7:43) 

since  thoy  do  not  progress  into  sector  1. 

For  Che  alternate  cases,  posed  by  the  last  two  of  (7:24), 
the  branches  of  0^  = 0 and  D^*=0  asymptote  to  the  limit 

C.  General  Shells  - Sector  2 

Recall  that  a great  deal  of  symmetry  was  effected  in 
writing  Che  sector  bounds  for  general  three  layered  shells. 

The  two  facings  are  assumed  materially  identical  and,  as  a 
result  of  this,  the  inequalities  (5:43)  result.  For  the 
pair  of  sector  2 bounds  given  as 


the  frequency  equation  in  sector  2 is  mathematically  equi- 
valent to  the  core  layered  cases  discussed  above  with  the 
exception  that  the  two  outer  layers,  being  identical,  cause 
Che  occurrence  of  equation  (7:40)  instead  of  (7:38)  for 
DjbQ,  Both  factors  of  (7:40)  represent  the  asymptotic  lo- 
cus m = b . The  frequency  equation  = 0 degenerates  into 


, both  of  which  represent  the  asymptote  m = b , 
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As  previously  shown,  the  above  forms  of  the  frequency  equa- 
tion are  valid  only  for  m very  near  its  asymptotic  limit 
value  and  f » 1. 

For  the  cases 


b'”  < n « V’ , 


the  sector  2 frequency  equation  is  mathematically  similar  to 
Che  sector  2 frequency  equation  for  sandwich  shells.  This 
being  the  case,  D3  0 yields  Che  locus  m»b“’  and  D^eO 
yields  the  locus  tn  = b''  for  the  asymptotic  limit  of  all 
branches . 

3,  The  Fraouenev  Eciuatlon  for  Higher  Sectors 

Having  hy  this  time  become  rather  familiar  with  the  gen- 
eral structure  of  the  frequency  equation,  it  can  be  shown 
now  most  easily  chat  all  branches  of  the  frequency  equation 
(5:34)  pass  into  sector  2 on  the(l-'(  plane  and  only  a max- 
imum of  four  ever  pass  below  Che  sector  2 lower  bound.  The 
oscillatory  nature  of  the  regular  Bessel  functions  may  be 
generally  regarded  in  Che  same  manner  as  Chat  of  the  sin- 
usoidal functions.  Consequently,  in  sector  7 of  the  frequency 
spectrum,  where  all  Bessel  functions  are  regular,  any  chosen 
ray  must  Intersect  an  infinite  number  of  branches  excepting, 
possibly,  those  with  zero  cut-off  frequency.  In  particular, 
the  ray  mcco  intersects  all  branches  at  their  nonzero  cut- 
off frequencies.  When  m attains  the  slope  of  the  lower 
bound  of  sector  7,  two  columns  (at  least)  of  equation  (5:34) 
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cease  to  be  oscillatory.  If,  when  this  happens,  (5:34)  is 
expanded  It  is  clear  that  there  remain  oscillatory  factors 
in  each  term  of  the  expansion.  Further,  ? is  the  only 
variable  in  this  expansion.  Consequently,  an  infinite  set 
of  values  for  C may  be  found  which  depict  the  passsse  of 
the  branches  across  the  ray. 

Frosresslng  in  a similar  manner  from  sector  7 to  6, 
and  so  on  to  3,  at  each  sector  loi.-er  bound,  the  expansion 
of  equation  (5:34)  retains  oscillatory  factors  in  every 
term.  There  is  always  the  possibility  of  selecting  an  in« 
finite  number  of  values  for  ^ which  satisfy  the  requirement 
that  the  frequency  equation  vanish. 

In  sector  2 this  behavior  continues  to  hold  until  the 
slope  of  the  ray,  along  which  the  examination  is  being  car- 
ried out,  attains  the  value  of  the  sector  2 lower  bound. 

When  this  occurs,  sll  oscillatory  componenta  of  (5:34)  van- 
ish, ? no  longer  occurs  as  a free  variable  but  appears  im- 
plicitly in  ffl.  The  detailed  investigations  in  sector  2 
have  shown  that  all  branches  approach  but  do  not  cross  the 
sector  2 lower  bound.  This,  of  course,  is  with  exception 
to  the  few  nondlsperslve  branches  which  were  shown  explicit- 
ly to  exist  in  sector  1. 

With  these  facts  in  mind,  it  is  clear  that  all  branches 
progress  into  sector  2.  One  reasons  as  follows:  when  m‘«o, 

an  infinity  of  branches  may  be  found  which  cross  the  ray. 

As  ra  decreases  there  is  always  an  infinity  of  branches  which 
are  countable.  An  infinity  of  branches  have  been  shown  to 
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asymptote  to  the  sector  2 lower  bound  slope,  say  m'.  Sup- 
pose, however,  some  subgroup  of  brsiichos  of  (5:3i),  ema- 
nating from  nonzero  cut-off  frequencies,  do  not  pass  to  the 
asymptote  m'.  Then,  there  is  a slope  ra'  below  which  this 
subgroup  does  not  pass.  If  the  subgroup  of  branches  is  an 
interior  group  of  the'  infinite  family  of  branches  of  (5:43} 
then  only  a finite  group  of  branches  (those  below  the  sub- 
group} pass  to  the  limit  slope  m'.  This,  however,  is  not 
possible  since  it  was  shown  that  an  infinite  number  of 
branches  run  to  in'.  Hence,  there  can  be  no  such  subgroup 
and  all  branches  nust  pass  to  the  asymptotic  slope  m'. 

4.  Suir-nnrv  of  the  Asymptotic  Behavior  of  the  Frecuenev 
Equation  \>hen  1 

The  most  important  result  of  the  foregoing  analysis  is 
that,  with  the  exception  of  some  of  the  lower  modes,  all 
branches  of  the  frequency  equation  approach  the  slope  of 
the  lowest  turning  point  line  for  any  choice  of  materials 
of  the  layers  and  for  all  n»0.  In  other  words,  when  the 
wavelength  is  very  short,  the  traveling  waves  propagate 
with  a phase  velocity  equivalent  to  the  least  shear  wave 
velocity  of  the  materials  of  the  cylinder. 

It  has  been  shown  that,  for  the  case  of  isotropic  lay- 
ers, and  for  very  small  axial  wavelength,  the  phase  veloci- 
ties of  at  most  two  modes  can  approach  the  Rayleigh  wave 
speed  for  waves  traveling  in  semi-infinite  continue  composed 
of  Che  same  material  as  the  inner  or  outer  layers.  Moreover, 
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If  the  material  constants  compl7  ulth  requirements  similar 
to  those  set  forth  by  Piguro  7,  Stoneley  waves  nay  exist. 
Then  it  is  possible  that  the  phase  velocities  of  at  most 
two  additional  modes  can  approach  asymptotically  the  wave 
speeds  for  Stoneley  waves  traveling  at  the  interface  be- 
tween semi-infinite  bodies  composed  of  the  materials  of 
layers  one  and  too  or  three  and  too,  vdien  Bota  Bene: 

pure  Itayleigh  or  Stoneley  waves  can  exist  only  in  sector  1 
of  thea-t  plane. 

Three  distinct  eases  may  be  considered. 

a)  The  upper  boundary  of  sector  1 is  the  turning  point 
line  with  slope  b"’.  For  this  case  there  is  always  a Ray- 
leigh wave  at  the  surface  of  Che  inner  layer,  traveling 
with  the  velocity  If  rajb'”  is  greater  than  t!",  no 

Rayleigh  wave  exists  in  the  outer  layer  of  the  shell  ex- 
cept as  Che  successively  higher  modes  terrace  on  the  Ray- 
leigh slope  In  passing  to  the  asymptotic  limit  t:”.  If 
b‘”  is  less  than  b"’  then  there  will  be  a Rayleigh  wave 
at  the  surface  of  the  outer  layer  for  the  case  The 

lowest  branch  of  the  frequency  equation  will  approach  as- 
ymptotically the  line  with  least  slope  (1  = 1 or  3) 

and  the  second  branch  will  approach  the  remaining  slope. 

If  the  material  parameters  ore  such  that  Stoneley  waves 
can  exist  at  the  interface  r = b,  Che  phase  velocity  of  the 
second  or  third  mode  will  approach  that  of  the  Stoneley 
wave  velocity  depending  on  the  ordering  of  this  velocity 
relative  to  the  two  Rayleigh  wave  velocities.  If  Che 
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materials  o[  layers  two  and  three  will  allow  the  presence 
of  Stoneley  waves,  then  the  second,  third,  or  fourth  mode 
will  approach  the  Stoneley  wave  velocity  for  the  materials 
at  the  interface  r cc.  This,  in  turn,  depends  on  the  re- 
lative magnitude  of  the  Stoneley  velocity  at  this  inter- 
face relative  to  the  previously  mentioned  three  velocities. 
If,  on  the  other  hand,  m^'b”  is  greater  b'" , there  can  still 
be  a Stoneley  wave  in  sector  1 for  the  r = c interface  of 
the  shell  since  the  Stoneley  wave  velocity  is  always  brac- 
ketted  by  the  least  Rayleigh  wave  velocity  and  the  least 
shear  velocity  of  either  material  involved.  Since,  It  is 
possible  for  mi’lS*’  to  be  less  than  b“\  It  is  possible  for 
a Stoneley  wave  to  exist  at  r = e. 

b)  Tho  upper  boundary  of  sector  i is  the  turning 
point  line  with  slope  b‘*.  In  this  case  the  observations 
of  paragraph  (a),  above,  hold  if  the  roles  of  layers  one 
and  three  are  interchanged. 

c)  The  upper  boundary  of  sector  1 is  the  turning 

point  lino  with  slope  b“\  When  this  is  the  case  there 
will  be  no  modes  which  approach  tho  Rayleigh  wave  apeeds 
if  both  (i  = l,3)  are  greater  than  b“’.  The  lowest 

two  modes  may,  hotrever,  approach  the  Stoneley  wave  velo- 
cities at  f - b and  r ° c,  providing  the  material  constants 
of  the  layers  allow  this  possibility  (see  Figure  7).  If 
Stoneley  waves  exist  at  either  interface,  the  velocities 
will  necessarily  be  intermediate  to  and  b^’. 

If  is  less  than  b‘  , a mode  will  approach  the 
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Rayleigh  wave  velocity  for  a semi-inf Inlce  body  composed  of 
the  material  of  the  inner  layer.  The  mode  which  approsches 
this  velocity  will  be  the  first,  second  or  third  depending 
on  the  ordering  of  the  two  Stoneley  wave  velocities  (if 
they  exist)  and  the  Rayleigh  wave  velocity.  Similar  obser- 
vations hold  if  only  m^*’b‘”  is  less  than  b*'’’.  If  both  m^’b''’ 
(iel,3)  are  less  than  b‘”,  then  two  modes  will  asymptote 
to  the  Rayleigh  wave  velocities  of  the  inner  and  outer  lay- 
er materials.  In  this  case  the  lowest  four  modes  will  ap- 
proach the  four  velocities  noted  (provided  the  Stoneley 
waves  can  exist)  and  their  ordering  will  be  dictated  by  the 
ordering  of  magnitudes  of  the  four  nondlspersive  velocities. 


CHAPTER  VIII 


1 . The  Frequency  Equation  for  Infinite.  Three  Layered 

If  the  mean  radius  of  the  cylinder  is  made  infinitely 
large,  the  cylinder  passes  over  into  a three  layered  plate. 
If  this  limit  is  taken  on  the  frequency  equation  it  can  be 
shown,  also,  to  pass  over  to  the  frequency  equation  for  the 
plate.  To  show  that  this  is  indeed  the  case,  it  is  neces- 
sary to  use  the  asymptotic  large  argunient  formulas  for  the 
Bessel  functions,  expand  the  determinant,  end  collect  terms 
in  factors  which  reduce  to  cos(thickness  x wave  number)  and 
sin(thickness  xwave  number).  Hith  an  18X18  determinant  the 
elegance  of  such  a limiting  process  loses  most  of  its  luster 
and  a more  primitive  method  will  be  employed. 

Using  the  transformation 

d-i  = dr,  dn  s vde  , di  . di  , (8:1) 

the  requirement  that  r»«  inposes  that  d/dy*0,  eliminating 
all  dependence  on  y.  The  coordinate  system  In  the  cylinder 
is  replaced  by  that  system  shown  in  Figtire  8.  Although  the 
dimension,  a,  is  now  zero,  the  notation,  a,  will  be  re- 
tained (formally)  to  preserve  the  similarity  between  the 
forthcoming  and  previously  obtained  equations. 

The  displacement  field  reduces  to  the  conditions  for 
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'nils  choice  of  solutions  conforms  to  the  large  argument 
asymptotic  behavior  of  and  W^.  It  will  be  convenient  if 
appropriate  differentiation  formulas  for  (8:7}  can  be  pro- 
vided so  that  two  separate  developments  (for  £ j > 0,6j < 0)  of 
the  boundary  and  interface  conditions  do  not  have  to  be  made. 
These  are  found  to  be 

and  (8:8) 


with 


iin  7'"‘i 


(8:9) 


Koto  the  similarity  of  (8:8)  with  the  Bessel  function  dif- 
ferentiation formulas  (5:32). 

Using  equations  (8:7)  and  (8:8)  with  (8:6)  and  (5:31)  in 
(8:2),  (8:3),  and  (8:A)  yields  the  basic  set  of  equations 
from  which  the  frequency  equation  Is  formed.  As  pointed  out 


Umli  detUijl  = 


c„  = k'  4>“ , 4"  f »- 

0„-  i:,.-t>, 


■ '^"S“ 

c«. 

c„ . 4”  - 4"  -jy'l'! 

C_-  -if!  14"4V* 

c„-  C«-0, 


185 


C„. 

Ct,’  -(4''*4")'?/s‘'’><-,<?/^"b\ 
c„*c,,-o, 

Css»  6r’?i‘'’'3,(’;V''«'^,  (8:U| 

C„  . ^,k'u;v~-4’$..4'-l-}.(K"a, 

c„  - ^,  U"s“»‘“  ^4"*  • 

c„.  Y 


5 t 


Sin('^6Ml 


(a)  (8:12) 


0. 


IS6 


<b)  (8:12) 


Hie  equation  B.  =0,  under  similar  conditions,  degene- 
rstes  into  the  frequency  equations  for  symmetric  and  anti- 
synmetrlc  Rayleigh-Lame  waves  given  by 


where  the  4-1  is  for  symmetric  waves  and  -1  is  for  the  anti- 
syirmetric  case.  These  wave  equations  have  been  discussed 
fully  by  Hindlin  [72]  and  Meeker  and  Heitaler 

Equatiqns  (8:12)  are  quite  similar  to  the  elements 
which  compose  (6:15),  hence  the  behavior  of  (8:11)  and  (6:15) 
should  bo  comparable.  Furthermore,  when  ?■*•“,  equations 
(6:15)  pass  over  into  the  identical  form  as  (6:12)  which 
bears  out  the  conclusion  that  when  the  radial  wave  number 
is  very  large  the  wave  phenomana  in  the  cylinder  approach  Che 
plato  as  a limiting  case.  There  is,  consequently,  no  need 
to  discuss  the  asymptotic  behavior  of  (8:11). 


(8:13) 


CHAPTER  IX 


1.  Method  of  the  Hurr.Brloal  Analva Is 

A computer  program,  uritten  for  execution  on  an  IBM 
7044/7094  CCS  System,  and  using  the  modus  described  in  the 
opening  paragraphs  of  Chapter  VII,  was  used  for  obtaining 
a numerical  evaluation  of  the  frequency  equation  (5:34)  and 
its  degenerate  forms  (6:7),  (6:14),  and/or  (6:21). 

The  program,  consisting  of  e calling  program  and  30 
subroutines,  uae  devised  in  a manner  such  that  any  aspect 
of  the  problem  could  be  successfully  treated  regardless  of 
the  degree  of  anisotropy  dictated  by  the  material  constants. 
Previously  written  end  tested  subroutines  for  computing 
Bessel  functions  of  integer  order  and  real,  imaginary,  or 
complex  arguments;  and  subroutines  for  determinant  evalu* 
ation  end  matrix  equation  solution  were  incorporated  from 
sources  other  than  the  author's  pen. 

Use  of  the  "ray"  analysis  furnishes  a time  saving  faC' 
tor  in  the  problem  execution  since,  once  the  angle  of  the 
ray  is  specified,  and  the  sector  of  the  plane  in  which 
it  lies  identified,  there  is  no  requirement  to  reascertain 
this  information  on  each  successive  iteration  of  <T  or  '7  un- 
til a new  ray  is  specified.  Furthermore,  those  quantities 
depending  only  on  the  ratio  A/f  (for  instance,  ss'' , 
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ily  be  determined 


outset  of  computations 


for  each  new  ray. 

The  program,  in  normal  sequence,  executes  the  cut-off 
(;«0)  computations  - in  which  roots  of  the  determinants 
^1’  ^5  n^O,  or  and  for  n>0  are  found 

and  displacement  distributions  determined.  Subsequently, 
starting  at  a preset  minimum  value  of  n or  ^ , the  fre- 
quency equations  Dj  and  for  n = 0 or  detlCj^jl  for  n^  0 
for  finite  wavelength  (t*0)  are  examined  along  rays  ema- 
nating from  the  origin.  The  starting  angle  of  the  ray 
with  steepest  slope  may  be  any  prechosen  value.  For  the 
curves  presented  In  Section  2 of  this  chapter  a maximum 
angle  of  88  degrees  to  tho^-axis  was  chosen.  Each  suc- 
cessive ray  is  declined  at  2 degree  increments  until  the 
angle  60  degrees  is  obtained.  Below  this  angle  successive 
radiais  one  degree  apart  are  examined  until  a radial  of 
sufficiently  shallow  angle  is  examined  such  that  the  fre- 
quency equation  has  no  zeros.  The  density  of  radiais  ex- 
amined is  controlled  by  input  data. 

Along  each  radial  the  zeros  of  the  determinants  listed 
above  are  detected  by  sign  change  information  of  the  deter- 
minants, When  a sign  change  of  the  determinant  is  noted 
for  two  consecutive  values  of  the  variables,  the  increment 
is  halved.  A test  is  executed  to  determine  which  half  in- 
terval the  root  lies  in  and  this  interval  Is  subsequently 
tenthed  until  a sign  change  in  the  determinant  value  re- 
occurs . Coordinates  of  the  zeros  thus 
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deemed  to  be  miduay  between  these  two  smallest  increments 
in  A or  over  which  the  sign  change  occurred.  Since  the 
maxiimuD  size  of  the  nornial  increment  is  controlled  by  in- 
put data,  the  veracity  of  the  roots  so  located  can  be  con- 
trolled to  any  degree  of  accuracy  sought.  For  the  figures 
presented  herein,  the  accuracy  of  a and  ? is  controlled  to 
within  t 0.0025  - a quantity  deemed  sufficiently  small  for 
plotting  accuracy. 

The  following  scheme  was  employed  to  regulate  the  in- 
crement length  along  a ray  • a distance  which  is  always 
greater  than  the  increment  set  off  along  the  coordinate 
axes.  The  variable,  either  ft  or^i  which  is  chosen  as  the 
independent  variable  along  a given  ray  is  determined  by  the 
inclination  of  the  ray  to  the  "^-axis . For  rays  of  incli- 
nation greater  chan  60  degrees,  Cl  is  the  independent  var- 
iable with  a normal  increment  of  &A  . For  raya  of  incli- 
nation 45  degrees  to  60  degrees,  fL  is  again  the  independ- 
ent variable  and  afl/2  the  normal  increment.  Similarly, 
for  rays  less  than  45  degrees,  tp  is  the  independent  var- 
iable and  the  Increment  is  = AU  for  rays  of  angle  loss 
than  30  degrees  or  Af/2  for  rays  between  30  degrees  and 
45  degrees. 

The  computation  along  a given  ray  is  terminated  when 
either  the  prescribed  number  of  roots  have  been  found,  or 
a preset  ray  length  has  been  obtained. 

Special  means  have  been  <n-itten  into  tho  program  so 


localized  high  density  variation  of  the  variable 
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eland  t and/or  the  ray  slope  may  be  obtained.  This  facil- 
ity is  required  whenever  tvo  branches  are  very  close  to- 
gether or  when  undispersed  inodes  (uhich  occur  when  '^>>1) 

Hie  displacements,  computed  only  for  the  cut-off  cases, 
have  been  found  in  the  following  manner.  For  n = 0 the 
(ivl,2,  or  3)  of  the  middle  layer,  appearing  in 
or  Dj,  has  been  arbitrarily  set  to  1.0.  For  the  ease  n>0 
the  integration  constants  A^'  and  a‘^’  were  chosen  as  1.0 
for  Che  plane-strain  and  the  longitudinal  shear  displace- 
ment distribution,  respectively.  The  distributions  of  dis- 
placement, once  found,  are  normalized  relative  to  the  larg- 
est amplitude  of  Che  particular  distribution  so  that  in 
each  instance  comparisons  can  be  more  easily  made  from  one 
case  CO  another.  To  further  facilitate  comparison  of  sim- 
ilar modes  of  vibration,  the  displacements  for  each  dif- 
ferent case  of  geometry  and/or  materials  is  displayed  with 
the  corresponding  mode  for  a homogeneous  cylindrical  shell 
composed  of  the  same  material  as  the  core  layer  end  total 
thictaess  of  the  composite  shell. 

The  program  was  checked  by  comparison  of  output  with 
Gazis'  C^lD  results  for  the  case  n = l,  and  comparison  with 
Hirsky's  C'^0]  results  for  the  case  n = 0.  Although  both 
test  cases  were  for  homogeneous  shells,  it  is  easily  recog- 
nized that  algebraic  or  sign  errors,  if  present  in  the 
three  layer  shell  equations,  would  eliminate  consistent  and 
accurate  results.  To  eliminate  the  remote 
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self  rectifying  errors,  a third  test  case  for  two  layered 
composite  shells  was  also  employed.  This  case  was  furnished 
by  Dr.  A.  E,  Armenakas  by  private  arrangement. 

The  program  described  above  was  viritten  exclusively  in 
double  precision  and  left  less  than  40  unused  locations  in 
the  7094  core  memory.  For  this  reason  it  was  necessary  to 
write  a second  program  for  the  express  purpose  of  deter- 
mining the  displacement  and  stress  distributions  in  the  com- 
posite cylinder  for  the  case  of  traveling  waves  (?#0). 

This  was  accomplished  utilizing  many  of  the  subroutines  pre- 
viously employed  - modified  suitably  to  treat  only  the  non- 
honogeneous  set  of  equations  which  evolve  for  computation 
of  the  Integration  constants  based  on  the  preknowlcdge  of 
the  roots  of  the  frequency  equation.  Due  to  computer  time 
limitations,  this  program  has  been  debugged,  but  not  run 
for  output  in  the  dissertation. 

The  data  governing  the  eases  presented  in  Section  2 of 
this  chapter  are  given  below.  The  first  data  case  is  for 
Poisson's  ratloi’  = 0.3  and  all  layers  of  the  same  material. 
This  esse  reflects  the  conditions  for  a homogeneous  shell 
and  Is  included  for  Che  purpose  of  comparison  with  Che 
three  layer  eases.  In  both  three  layer  cases  the  inner 
and  outer  shells  ore  treated  as  of  the  same  thickness  and 
material  composition,  and  all  three  layers  are  isotropic 
.with  i5"’=0.3  (i  = l,2,3).  Moreover,  results  have  been  ob- 
tained only  for  the  axisymmetric  case  (n  = 0) . 
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1 = 3.5, 
‘^1.5, 


dj  = l.S, 
Case  II 


For  each  set  of  mecerlol  constants  noted  above  the  fol- 
lowing two  sots  of  geometric  constants  were  employed: 

H/R  = 1.0  and  H/R  = 0.L,  where 


defines  the  total  thickness. 

The  first  geometry  corresponds  to  the  case  of  a thick 
shell  and  the  second  to  chat  of  a thin  shell.  The  chin 
shell  ease  has  a total  thickness  to  radius  ratio  (chat  is, 
H/R  = 0.1)  which  is  within  the  range  of  validity  of  the  thin 
shell  bending  theories  [95,96],  Furthermore,  it  is  expected 
that  this  thin  shell  case  may  fall  within  Che  range  of 
validity  of  the  sandwich  theory  developed  by  Yu  [41,76,77]. 

Each  combination  of  material  constants  versus  H/R  ra- 
tio was  executed  for  the  single  choice  of  thin  outer  layers, 
as  depicted  by  the  ratio  A" =0.1.  Consequently,  a total 
of  six  frequency  spectra  were  computed,  each  showing  the 
first  ten  branches  of  the  exisyimetrie  frequency  equation. 
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All  of  the  frequency  spectra  have  been  plotted  to  a 
value  of  7 =5.0.  This  value  of  ( ucs  chosen  since  It  repre- 
sents 0 value  of  the  half  Havelength  L which  Is  considerably 
smaller  than  any  dimension  of  the  shell.  For  Instance,  In 
terms  of  the  outer  diameter  to  half  wavelength  ratio,  that 
la,  d/L,  the  abscissa  valuo  of  '^  = 5.0  corresponds  to 
d/L  = 62.9  and  d/L  = 9.0,  respectively,  for  the  two  geometries 
H/R  = 0.1  and  1.0,  when  h'Vl^*'  =0.1.  The  maxliraim  value, 

9 = 5.0,  Is  sufficiently  large  so  that  the  outset  of  the  as- 
ymptotic effects  discussed  in  Chapter  VII  can  be  observed. 
This  is  reasonable  since  when  5 it  is  evident  that  the 
shell  thlclcnesses  are  beglimlng  to  assume  very  large  propor- 
tions relative  to  the  half  wavelength  of  the  axially  travel- 
ing waves • 

2.  Pis cuss ion  of  the  Freouenev  Spectra 

Since  the  frequency  spectra  are  plotted  over  such  large 
ranges  of  values  of  the  variables  Cl  and  'p,  it  has  been  found 
advantageous  to  plot  the  spectra  in  two  separate  figures  (for 
each  case).  The  first  of  any  such  set  of  figures  pertains  to 
the  long  wavelength  region  of  the  variables , and  the  second 
to  the  mid-wavelength,  or  almost  asymptotic  range  of  varia- 
bles. The  inserts  in  Figures  15,  17,  and  19  are  for  the  low- 
est branches  of  the  frequency  equation  near  the  origin,  end 
have  been  so  included  for  the  sake  of  clarity.  Small  por- 
tions of  the  ninth  or  higher  branches  have  been  excluded  for 
the  range  of  values  of  C-  near  2.5  duo  to 


the  method  of 
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presentation  o{  the  curves,  but  the  omissions  obscure  no 
unusual  behavior  of  the  frequency  equation.  Due  to  the 
manner  in  which  the  computer  program  counted  solutions  of 
the  frequency  equation,  fragments  of  the  eleventh  and 
twelfth  branches  were  obtained.  These  fragments  were  in- 
cluded to  complete  the  orientation  of  the  branches  in  the 
cut-off  region. 

Cases  I end  11  both  fall  in  the  class  of  general 
shells  described  by  inequalities  (5:63).  The  ordering  of 
the  turning  point  lines  is  such  that 


Both  case  I and  II  yield  the  same  values  for  the  b'’ 
and  '<3  (1  = 1,2,3)  and  Che  apectra,  therefore,  present  a 
means  of  comparing  the  separate  effects  of  different  stiff- 
nesses as  well  as  densities  when  these  types  of  variation 
each  lead  to  en  analytically  identical  form  for  the  fre- 

For  ease  of  identification  of  the  branches  of  the  fre- 
quency equation  the  following  notation  has  been  used  to  i- 
deiitify  the  branches.  Since  only  the  axisymmetrlc  condition 
is  under  study  all  branches  originate  with  cut-off  frequen- 
cies determined  by  solutions  of  D^,  and  D^.  When  T#0 
the  branches  are  determined  by  end  D^.  Consequently, 
all  branches  emanating  from  the  solutions  of  Dj('<;  = 0)  and 
continued  by  solutions  of  D3(ty0)  will  be  called  torsional 
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branches  and  demarked  '^^>^2’  * * * > according  to  the  order 
of  the  cut-off  frequent  of  the  torsional  branch  relative 
to  other,  lower  torsional  branches.  In  a similar  manner 
Che  branches  emanating  from  solutions  of  Dj^  and  and 
perpetuated  by  will  be  called  noncorsional  branches  and 
demarked  where  i indicates  the  order  of  the  branch 

relative  to  nontorsionsl  brandies  with  smaller  cut-off  fre- 
quencies, and  j is  1 if  the  branch  is  a longitudinal-shear 
mode  at  its  cut-off  frequency  or  2 if  radial  plane-strain. 

A study  of  the  six  spectra  reveals  the  following 
generalities : 

a}  For  7 < 1 the  branches  of  the  frequency  equation 
are  comparable  for  oil  eases  of  the  same  geometry,  regard- 
less of  Che  material  constants  of  the  layers,  and; 

b)  For  7 > I the  branches  of  the  frequency  equation 
are  comparable,  regardless  of  Che  geometry,  for  each  of  the 
three  cases  of  material  constants. 

These  observations  indicate  that  the  shell  geometry 
plays  a strong  rSle  in  the  determination  of  branches  in 
Che  long  wavelength  region  and  this  influence  decreases  as 
Che  wavelength  decreases.  Furthermore,  Che  role  of  Che 
material  constants  appears  more  pronounced  in  Che  short 
wavelength  region, 

Tuening  now  to  more  specific  observations,  Che  long 
wavelength  region  will  be  considered  initially.  It  is  clear 
from  Che  identification  of  the  branches  of  each  of  the  six 


spectra  what  Che  nature  of  the  modes  of  vlbratior 
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the  cut-off  frequencies  of  each  branch.  The  displacement 
distributions  of  the  first  six  modes  of  vibration  have  been 
illustrated  in  Figures  21,  22,  and  23  and  will  be  discussed 
in  Section  3. 

Comparison  of  the  three  layer  spectra  with  the  homo- 
geneous spectra  indicates  that  the  first  longitudinal  shesr 
mode  initiates  with  a phase  velocity  similar  to  the  phase 
velocity,  J E/p  , for  isotropic  homogeneous  shells.  It  is 
further  noted  that  this  branch  is  the  only  one  which  does 
not  initiate  in  the  highest  sector  of  the  fl-C  plane.  Since 
the  outer  materials  are  identical  there  are  only  five  sec- 
tors on  the  Cl-'5  plane  and  the  branch  always  Initiates 

in  sector  3.  This  requires  that  the  radiel  wave  munbers 
g“’  (i  = 1,2,3)  be  imaginary.  Consequently,  the  dllatatlonal 
potentials  tp‘2’  (i  = l,2,3)  contribute  exponentially  decaying 
components  to  both  the  wave  equation  and  the  displacement 
field. 

The  radial  extensicpnal  branch  ® pronounced 

minimum  in  each  of  the  six  spectra.  Superficially,  this 
trough  indicates  a negative  phase  velocity  (see  equation 
(5:39)).  However,  it  is  well  known  (see,  for  instance, 
[72])  that  the  above  conclusion  does  not  apply.  Since  the 
axial  v;ave  nufi^er  f appears  in  the  frequency  equation  only 
in  squared  form,  it  is  immediately  obvious  that  spectra 
similar  to  those  exhibited  in  Figures  9 through  20  could 
be  plotted  for  -f.  Furthermore,  it  has  been  shown 

ocal  minimum  of  a branch  (for  real  7)  of  the 
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frecjueney  equation,  another  branch  eiianates  which  ectisflea 
tlie  frequency  eqtiation  for  real  0.  and  complex  % • These 
new  branches  either  decrease  ChroUQh  complex  values  of  ^ 

CO  another  local  minimum  of  a branch  which  le  defined  for 
real  X only,  or  decrease  to  the  complex  X plane,  for  which 
fL“0.  In  this  manner  it  is  possible,  throueh  branches  of 
the  frequency  equation  satisfied  by  complex,  negative,  or 
positive  values  of  Xt  begin  with  Cl  = 0 and  define  patlis 
along  the  various  branches,  which,  for  continuous  monoton- 
Ically  increasing  fl,  always  have  zero  or  positive  group 
velocity.  Thus,  that  portion  of  any  other  branch, 

as  well)  with  negative  slope  msy  be  looked  on  as  indicating 
that  the  mirror  reflection,  about  the  A axis,  is  the  appro- 
priate segment  to  be  chosen  in  a complete  description  of 
the  mode  curves.  Furthermore,  it  necessarily  follows  chat 
the  branch  is  the  only  branch  of  the  frequency  equa- 

tion which  initiates  with  the  value  f = 0 when  A = 0.  All 
other  branches  originate  (at  A = 0)  with  some  complex  value 
of  T. 

The  cut-off  frequencies  of  all  modes  ere  raised  moder- 
ately with  increased  density  of  the  inner  layer,  and  only 
slightly  effected  by  increasing  the  stiffness  of  the  outer 
layers.  This  is  seen  from  comparison  of  Figures  11  and  13 
with  9,  or  17  and  19  with  15.  As  previously  remarked,  the 
geometry  of  the  shell  more  extensively  effects  the  cut-off 
frequencies  and  it  is  observed  that  as  the  H/R  ratio  de- 
creases for  a given  set  of  mate 
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and  modes  initiate  at  notably  smaller  cut-off  fre- 

quencies, vhen  compared  uith  H/Rsl.O. 

&cairinin8  next  the  reQion,  of  all  spectra,  for  '^>1, 
several  interesting  facts  arc  revealed.  Consider  first 
the  two  sets  of  figures  for  Che  homogeneous  cylinders,  Chat 
is  Figures  9,  10,  15,  and  16.  Two  important  results  are 
immediately  observable.  The  lowest  nontorsional  modes, 

HT^^  and  NTjj  both  run  to  the  Rayleigh  locus  a=0.927^. 

This  locus  is  cliaracterized  by  the  ray  inclined  at  an  an- 
gle of  42.83  degrees  to  the  ( axis  and  is  indicated  in  the 
figures  by  the  broken  line.  Notice  Chat  HT^^  approaches 
the  Rayleigh  slope  from  below  and  NT22  from  above.  The 
analysis  conducted  in  Chapter  VII,  when  suitably  altered 
for  homogeneous  cylinders , indicates  this  behavior  of  both 
NTjj^  and  RT22’  lowest  nontorsional  mode  with  nonzero 

cut-off  frequency  (NT22)  approaches  the  Rayleigh  locus  as 
5*09,  Prom  the  fact  that  KTj^^  never  lies  in  the  highest 
Eoctor  (sector  3 for  homogeneous  cylinders},  a similar 
analysis  easily  yields  the  same  conclusion  for  this  branch 
of  the  frequency  equation. 

The  second  fact  of  significance  with  regard  to  Che 
homogeneous  cylinder  spectra  is  the  Cviinning  of  Che  T.  and 
NT^^j^  j 'M^anches.  It  is  noticed  that  decreasing  H/R  ratio 
brings  these  curves  Into  contact  when  H/R  = 0.1.  Apparently, 
for  1!/R<0.1  Che  T^  branches  will  drop  below  the  ^ 

branches.  It  has  been  observed  Chat  Che  contact  point  be- 
tween the  branches  for  H/R  = 0,1,  and  Che  locales  of  closest 
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approach  for  H/R=1.0,  fall  along  the  rayn.  = f?'f.  This  Is 
the  required  relation  of  fl  to  '?  for  the  existence  of  Lame 
[51]  modes  in  plates.  However,  Annenakas  [25]  has  ob- 
served that  the  Lame  modes  con  not  exist  in  composite 
shells. 

The  spectra  for  three  layered  shells  (Figures  11 
through  14,  and  17  through  20)  show  that  no  mode  decreases 
below  the  lowest  turning  point  line  0.  = ';  (recall  that 

1).  This  must  be  the  case  since,  as  was  shown  in  Chap- 
ter VII,  Rayleigh  waves  eon  exist  only  if  b“’  (i  = l,3)  < b“’. 
Furthermore,  the  choice  of  material  constants  will  not  allow 
the  existence  of  Stonoley  interface  waves. 

Of  paramount  concom  is  the  effect  on  the  branches  as 
the  stiffnesses  of  the  outer  layer,  or  density  of  the  core 
is  increased.  Increasing  the  core  layer  density  retards 
the  rate  of  approach  of  all  branches  to  the  limit  11  = '7.  Ko- 
tlce,  as  well,  that  all  hint  of  twinning  of  the  branches  is 
removed.  The  extent  of  the  approach  of  the  upper  bronchos 
to  the  asymptotic  limit  is  not  greatly  effected  by  the  val- 
ue of  H/R. 

When  comparing  the  spectra  of  Figurets  13,  14,  or  19, 

20  with  the  corresponding  homogeneous  spectra,  it  is  at 
once  noticed  tliat  not  only  is  the  twinning  tendency  obviated, 
but  furthermore,  the  Tj^  branches  intersect,  from  below,  the 
^^i+2  j branches.  This  is  a pronounced  difference  In  the 
small  wavelength  behavior  of  the  mode  curves.  Comparison 
of  the  slopes  of  the  two  families  of  branches  indicates 
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that  increaBing  the  outer  layer  stiffnesses  results  In  re- 
ducing the  rate  of  approach  of  Che  T.  curves  to  the  limit 
and  has  little  effect  on  Che  rate  of  approach  of  the 
curves  - when  compared  with  the  homogeneous  spectra. 
Comparison  of  Figure  12  and  14  or  18  end  20  show  that  the 
branches  approach  their  asymptotic  limit  more 

rapidly  in  the  shells  with  stlffer  outer  layers,  and  the 
branches  more  slowly. 

In  surmary  It  Is  pointed  out  that  the  cut-off  fre- 
quencies are  strongly  effected  by  Che  ratio  K/R  - approach- 
ing the  appropriate  place  frequencies  as  H/R-*-0,  and  the 
solid  cylinder  frequencies  (Fochhammer  frequencies  for  the 
homogeneous  case)  as  H/R->2.  Furthermore,  as  Che  stiffness 
of  the  outer  layers  or  density  of  the  core.  Is  increased, 
Che  rate  of  approach  of  the  branches  to  the  limit  la 

retarded,  relative  to  the  homogeneous  case;  it  being 
slightly  more  pronounced  for  small  H/R  chan  for  largo. 

3.  Discussion  of  the  Vibration  Displacement  Distributions 

The  identification  of  the  branches  of  the  frequency 
spectra  in  Figure  9 through  20  indicate  immedlacel.y  the 
nature  of  the  displacements  when  ( = 0.  The  branch  curve 
does  not  represent  a vibratory  mode  when  f is  tero 
since  the  frequency  is  also  zero.  The  branch  Tq,  noted 
on  the  two  homogeneous  spectra,  is  for  the  classical 
first  torsional  mode  given  as  Ug » (Ar}cosni.  This 
mode  is  not  obtained  from  the  solution  given  in 
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Chapter  V for  composite  eyllndcrs. 

Since  (for  ^ = 0)  does  not  represent  a inode  of 

vibration,  NT22  the  only  mode,  for  composite  shells, 
which  does  not  possess  a nodal  cylinder.  This  is  a so 
called  ring  mode  [22]  and  has  been  termed  the  breathing 
mode.  Ail  longitudinal-shear  and  circumferential-shear 
modes  possess  one  or  more  nodal  cylinders  and  have  been 
characterized  as  thickness  modes  Moreover,  the  se- 

cond and  higher  radial  plane-strain  modes  of  vibration  are 
also  thickness  inodes. 

The  first  two  modes  of  each  of  the  three  types  men- 
tioned above  are  displayed  in  Figures  21,  22,  and  23.  The 
modes  have  been  plotted  in  the  manner  previously  discussed 
in  Section  1 of  this  chapter.  Only  the  two  cases,  for  each 
geometry,  governed  by  homogeneity  and  the  case  1 material  con- 
stants have  been  plotted,  for  reasons  of  clarity.  The  mode 
shapes  for  the  cose  11  material  constants  are  generally 
bounded  by  the  pairs  of  curves  given  in  each  figure. 

In  comparing  the  curves  in  each  figure,  of  H/R  = 0.1 
and  H/R^l.O,  it  is  immediately  apparent  that  the  mode  shape 
curves  show  almost  no  attenuation  in  amplitude  with  respect 
to  thickness  for  the  ratio  H/R  = 0.1,  This  is  directly  the 
result  of  the  proximity  to  plate  geometry  that  this  small 
value  of  II/R  represents.  Furthermore,  an  additional  con- 
sequency  of  this  approadi  to  plate  conditions  is  evinced 
by  the  symmetry  of  the  distributions  relative  to  the  center 
plane  of  the  shell  for  small  K/R. 
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Figure  22,  for  H/R  = 1.0,  cotitaijis  the  only  plot  of  a 
case  II  data  displacement  distribution.  The  case  I distri- 
bution for  this  mode  nearly  coincides  with  the  homogeneous 
case.  Consequently,  compared  to  the  homogeneous  case,  the 
relative  stiffnesses  of  the  layers  Influence  the  ring  mode 
displacement  distribution  more  greatly  than  the  relative 
dens It lea. 

Study  of  the  displacement  distributions  for  all  three 
eases  of  material  constants , and  for  each  case  of  geometry 
reveals  that  the  effects  of  greater  outer  layer  stiffnesses 
is  similar  to  the  effects  of  greater  core  density  with  re- 
gard to  the  thickness  mode  shapes  near  the  surfaces  of  the 
cylinder.  The  outer  layers  of  the  composite  shells  have 
distributions  with  less  curvature  than  for  the  outer  ex- 
tremities of  the  homogeneous  shells.  This,  in  the  case  of 
stiff  outer  layers,  follows  without  explanation.  For  the 
case  of  shells  with  more  dense  cores,  the  result  is  due  to 
the  greater  inertia  of  the  core  layer  which,  on  the  basis 
that  all  three  layers  have  the  same 
observed  effect. 


stiffness,  implies  the 
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PIGUP.fi  10.  Ftequeucy  Spactirum:  Homogan.  (cont), 
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FIGURE  11.  Frequency  Spectrum!  C*se  I,  H/R  = l,0 
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FIGURE  12.  Frequency  Spc< 


i:  Case  I (eont),  K/R  = 1.0 
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FIGURE  13,  Frequency  Speecrum:  Case  II,  H/R*=1,0 
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FIGURE  14.  Frequency  Spaotrumr  Cnse  II  (cent),  H/Rol.O 
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FIGOiaC  15.  Freqvonoy  SpactKiia:  Homogen.,  H/Rr0.1 
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FIGUKiC  16.  Froquojicy  Spoetrun:  Ho.moeen.  (tent),  H/R^O.] 
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FIGURE  17.  Frequency  Spectrum;  Case  I,  H/R  = 0.1 
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FIGURE  18,  Frequency  Spectrum:  Cass  1 (cent),  H/R  = 0.1 
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FIGURE  19.  Frequency  Spectrum:  Cass  II,  H/R»=0.1 
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FIGUiiE  20.  Frequenej'  Spoctrun:  Case  II  (cont),  H/R»0.1 
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FIGURE  22,  Radial  Plane-Strain  Vibration  Distributions 
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FIGURE  23. 


Circumferential-Shear  Vibration  Distributions 
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